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ABSTRACT 



This report covers work done on performance analysis of a 5 inch, 

54 caliber gun-launched guided projectile with solid fuel ramjet (SFRJ). 

A computer program (TRAJET) was developed. The program contains 
ramjet and trajectory analysis. The ramjet part considers conical 
shock wave losses, inlet boundary layer losses, normal shock losses, 
subsonic diffuser recovery, expansion into combustor losses, heat losses 
at the combustor and nozzle losses. 

A flat earth trajectory with drag and thrust was considered. The 
various drag coefficients which were considered are: cowl drag coefficient, 

skin drag coefficient, wing (or fin) wave drag coefficient and wing (or 
fin) friction drag coefficient. Base drag is assumed to be zero due 
to the jet from ramjet nozzle. 

It was found that the cowl drag coefficient has a major influence 
on the results. Therefore, a separate program (AERO) was developed to 
calculate this parameter. 

The 5"/54 solid fuel ramjet has a capability to produce fuel 
specific impulse in the order of 400-900 sec. depending mostly on the 
flight altitude. The thrust coefficient varies in the range of 0.3±0.1 
depending on the internal areas. 

A range in the order of 50 miles can be achieved with the ramjet 
operation compared to only 13 miles achieved by the conventional pro- 
jectile. At low-altitude launch, a range of over 18 miles can be reached 
with the ramjet version. Launches at high elevation angles can be 
useful in air-defense scenario. The ramjet propelled projectile reaches 
the ranges mentioned above at high Mach numbers (Mq> 1.8). It is, therefore. 
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clear that the ramjet concept provides significant improvement and has 
an Anti Ship Missile Defense (ASMD) capability. 
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1. INTRODUCTION 



1 . 1 Background 

This report covers work done on propulsion and flight me- 
chanics of a gun-launched guided projectile. Gun-launched guided pro- 
jectiles have been developed by Martia Marietta for the U.S. Army 
(Copperhead) [1,2] and more recently, also, for the U.S. Navy [3,4] 

(5"/54 Mark 46). The USN round has solid rocket propulsion. 

The addition of a liquid fuel ramjet (LFRJ) to the Navy's 
version, was examined in the past by Brown [5]. This report concentrates 
on the addition of a solid fuel ramjet (SFRJ), instead. 

It is believed that solid fuel ramjet has some potential 
advantages compared with the liquid fuel. Some of these advantages 
are: 

-Simple design 

-High reliability in operation 
-Low cost 

-Fuel control system not needed 

On the other hand, there are also a few disadvantages to SFRJ compared 
to LFRJ. Some of these are: 

-Difficult to control magnitude of thrust 
-Difficulties in achieving high combustion efficiencies 

In both cases, the addition of propulsion improves dramatically, 
the performance of the projectile by multiplication of range and enhanced 
maneuverability. Even more; to operate and produce thrust, the ramjet 
engine depends only on its forward motion at supersonic speeds and does 
not employ any moving parts. This fact, which is especially emphasized 
in SFRJ, leads to some advantages of the ramjet concept over the other 
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propulsion alternatives, at supersonic speed. On the other hand, the 
ramjet engine requires an auxiliary booster to accelerate it to its 
supersonic operating regime. The boost required causes some system 
difficulties. But, while solving these problems, the ramjet system 
becomes even more attractive for use with gun-launched guided pro- 
jectile, like the U.S. Navy 5"/54. 

A computer program was developed to analyze the performance 
of the SFRJ. The computer program was written for the IBM-370 
computer at the Naval Postgraduate School, Monterey, California. HTPB 
was selected as a fuel, but performance with any other fuel can be tested, 
using the same model. A flat earth trajectory with drag and thrust was 
considered. Using solid fuel, a thrust-equal -drag trajectory is more 
difficult to achieve with SFRJ. Therefore, most of the results given 
in this report eliminate this case, and the exact value for drag, at 
each point, was calculated. However, if desired, it is possible to 
change the air mass flow, in order to obtain thrust-equal -drag flight. 

The computer program can calculate this case also. 

1 .2 Basic Concept 

The ramjet engine consists of an air inlet, which serves as 
a diffuser, a combustion chamber, and an exhaust nozzle [6,7]. The 
diffuser admits air to the engine, which is mixed with fuel (solid or 
liquid) at the combustor. After the burning process, which adds heat 
to the flowing air within the system, the gases are transfered to the 
nozzle. The nozzle converts part of the thermal energy into kinetic 
energy to produce thrust. 

The areas inside the ramjet engine are usually divided into 
six stations, as illustrated in Figure 1.1. Station 0 defines the cross- 



2 



PORT AT ENTRANCE 
TO COMBUSTOR 




o 



Figure 1.1 Schematic View of a Solid Fuel Ramjet 



5-INCH GUIDED PROJECTILE 




Figure 1.2 U.S. Navy 5"/54 Semi Active Laser Guided Projectile (SALGP). 



section area of the stream tube captured by the inlet of the projectile. 
Stations 1 and 2 identify the diffuser. Station 1 itself is at the 
throat of the diffuser, but sub-station 1C is after the conical wave; 
sub-station 11 and SI are ahead of the normal shock wave, located at 
the throat, or at the actual place, respectively; sub-stations 12 and 
S2 are as above, but behind the normal shock wave. Stations 3 and 4 
refer to the entrance and to the exit of the fuel grain within the 
combustor, respectively. Note that both and increase with time 
as fuel burns. Station 5 and 6 belong to the nozzle's throat and to 
the exit of the nozzle, respectively. 

1 .3 Data 

1.3.1 Dimensions 

In order to be compatible with the Navy's 5 inch, 54 caliber, 

Mark 46 gun mount, as modified for gun launched guided projectiles 
(Figure 1.2), a set of requirements were adapted initially. These were: 

a. External shape of existing 5" guided SAL projectile 

b. Length - 60.97" 

c. Length of combustion chamber - 23" 

d. Total weight - 104.7 lb. 

e. Muzzle velocity - 2500 ft/sec. 

Typical values for the internal areas in the ramjet within the 
projectile are (units - sq. in.): 

^;o ^1 ^2 ^3 ^5 ^ 

19.3 5.2 2.6 4.3 8.2 7.5 13.1 

Refer to Figure 1.1 for definition of the areas. A^ is a reference area. 

For typical flight Mach number of: Mg = 3, the appropriate Mach numbers 

at the main stations are typically: 
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2.1 0.56 0.3 

1.3.2 Combustion Process 



0.1 



Some of the losses in the total pressure were taken to be 
constant. These are: 

a. Inlet boundary layer losses typically = 0.93). 

b. Subsonic diffuser recovery typically = 0.93). 

c. Nozzle losses (t^; typically = 0.96). 

These values are typical and are not expected to vary too much. 

All the other losses in total pressure, which are dominant to 
the projectile performance,were calculated. These are: 

a. Conical shock wave losses (tt^). 

b. Normal shock losses • 

c. Losses due to expansion into the combustion chamber (tt^). 

d. Heat losses in the combustion chamber (tt^). 

See section A 1.4.2 for definition of various tt . Combustion efficiency 
was taken constant: Hj = 0.90 

Air heat capacity ratio was also taken constant: y, = 1.4; however, the 

a 

value for the gas heat capacities ratio of the combustion products (Y.p) 
was calculated from thermodynamic data for Hydroxy Terminated Polybutadiene, 
HTPB, burned in air [8]. The stoichiometric chemical reaction of HTPB 
burning with oxygen is as follows: 



4 Cy3HjQ30j + 393 0^ ^ 292 CO^ + 206 H2O (1.3.1) 

1.3.3 Trajectory 

In the trajectory part of the program, the various drag coef- 
ficients were calculated. Those are: 

a. Cowl drag coefficient (Cq|^), [9-12]. 
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b. Skin drag coefficient (Cq^), [13-15]. 

c. Wing (or fin) wave drag coefficient [16 - 20]. 

d. Wing (or fin) friction drag coefficient (C^^^p), [13, 19]. 
Base drag [13, 21-23] is assumed to be negligible due to the jet from ramjet 
nozzle. 

The model, which was chosen to calculate the cowl (nose) drag 
coefficient, was based on a theoretical development done by T. H. Gawain 
[9]. The modified program (AERO) is listed in Appendix G. However, 
program AERO as it is, appears to be too long to be used directly in the 
main program (TRAJET). Therefore, best fit curves for calculated results 
from AERO were used in TRAJET. 

Skin drag coefficients were calculated for either laminar or 
turbulent flows. The same routine was also used to calculate wing or fin 
friction drag coefficients. To calculate the wing wave drag coefficients, 
a psuedo 3-dimensional model was developed. 

The program also has an option to calculate a trajectory of a 
projectile without propulsion. In this case, the drag coefficients which 
are calculated are; 

a. Nose drag coefficient (a different model than the above). 

b. Base drag coefficient 

c. Skin drag coefficient (as above). 

1.3.4 Air Defense Scenario 

In the air defense scenario, the program takes into account only 
cases in which the projectile exceeds a Mach number of at least 1.8. This 
value of minimum Mach number (XMg) can easily be changed. 

1.4 Results 

Each section of the program was developed, tested and run sep- 
arately. The ramjet part was first run without the trajectory part using 
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constant altitude (Typically - 10,000 ft). The same was done with the 
trajectory part using vacuum case, thrust-equal -drag flight, or constant 
thrust case. The final version was programmed to give the following 
optional printings: 

a. Loop on all possible values of Ag/A^ and A^/A^ and print 
summary tables only. 

b. Print detailed, time dependent, tables for any specific area 
ratio chosen: 

- Results from combustion process (file name: CMB D) 

- Results from trajectory process (file name: TRJ D) 

- Various drag coefficients (file name: DRG D) 

c. Detailed print of every step during the calculation, for 
checkup. 

d. Variation of the above: 

Detailed print of cases that were found not to be 
suitable: 

- Reasons only 

- Full detailed parameters 

- Loop on Mach numbers, also (output of subroutine 
CALCM) 
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2. TYPICAL RESULTS 



Figure 2.1 presents the dependence of the fuel specific impulse 
in sec.) of the ramjet on the projectile Mach number at various altitudes. 

It appears that the 5"/54 ramjet has a capability to produce fuel specific 
impulse in the order of 400 - 900 sec., depending mostly on the flight al- 
titude. The dependence of I^p on the flight Mach number is weak. 

Figures 2.2 and 2.3 present the dependence of the thrust coefficient 
(C^) on the internal area ratios Ag/A^ and Ag/A^. In figure 2.2, the change 
of with altitude and with Mach number is also presented. The thrust coef- 
ficient (C^) varies in the range of 0.3±0.1 while Aq/A^ changes from 0.25 to 
0.40 and Ag/A^ changes from 0.42 to 0.26. 

The correlation between the fuel specific impulse (l^p) and the thrust 
coefficient (C^) is presented in figures 2.4 and 2.5. In both figures, a Mach 
number of Mg = 3.0 was selected. In figure 2.4, the correlation was checked 
at various altitudes and at various Ag/Ag area ratios. In figure 2.5, various 
internal area ratios (Ag/A^, Ag/A^) are presented. 

More detailed results are presented in Appendix H. The dependence of the 
projectile performance on the other internal area ratios (A.|/Ag, A2/Ag, Ag/A^) 
was also checked. Some typical results are presented in that Appendix. 

An altitude vs range dependence for various elevation angles is presented 
in figure 2.6. A range of over 80 km can be achieved with the ramjet operation, 
compared to only slightly more than 20 km achieved by the conventional pro- 
jectile. A low-altitude launch (in this figure, an elevation angle of 15’ 
was selected) is also presented reaching a range of over 30 km with the ramjet 
operation. The high elevation angles are mostly used in air-defense scenario. 
The drag of the projectile when the ramjet is not operating, for example. 
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after burnout, was not determined. The computer program does not account 
for the drag increase due to the ramjet not operating. Consequently some 
of the trajectory curves in figure 2.6 are in error. However, for tra- 
jectories at low gun elevation, the ramjet burns all the way to splash. 

These trajectories are accurate. For all trajectories, the curves are 
accurate to the point of ramjet burnout. The trajectories of interest to 
air defense are accurately calculated. Trajectory of thrust-equal -drag 
(vacuum) case is shown for comparison. 

Results for air defense scenario are presented in figure 2.7. Only 
ranges where the projectile Mach number exceeds at least 1.8 where con- 
sidered. The area ratios where chosen as specified in the figure. The 
two cowl angles, shown in figure 3.1 were 20° and 9.5° respectively. 

The gun elevation angle was varied from 7° to 80°. The change of atmospheric 
conditions with the altitude was taken into account. In table 2.1, some 
typical results for "Surface-to-Surface Mission" are presented. 
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Figure 2.1 Solid Fuel Ramjet; Dependence of Fuel Specific Impulse (I 

on Flight Mach Number at Various Altitudes 

Conditions: An/A -n oc n ^ 

0//^r-0*25, Aj/Aq = 0.47. A2/Aq = 0.827 

A 3 /A =0.427, A(-/A^=0.28, A./A^=l, e = 45 ° t = 0 
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Figure 2.2 SFRJ: Dependence of Thrust Coefficient on 

Internal Area Ratio (A«/A„) at Various Altitudes 
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Figure 2.3 SFRJ: Dependence of Thrust Coefficient on 

Inlet and on Nozzle Area Ratios 
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Figure 2.4 SFRJ; Dependence of Fuel Specific Impulse (I 
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Figure 2.5 SFRJ; Dependence of Fuel Specific Impulse (I 
on Thrust Coefficient at Various Internal Area Ratios 
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20 1 performance boundary due to ramjet unstart 
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TABLE 2.1 



Solid Fuel Ramjet Propelled 5"/54 Projectile 
Surface-to-Surface Mission 
Ranges (km) vs Gun Elevation Angles 



Elevation 

Angle 


7 ° 


25° 


45° 


65° 


o 

O 

CO 


a. Ramjet^ 


15.6 


49.9 


80.3 


15.9 


5.6 


b. Projectile 
Without 
Propulsion 


9.2 


17.5 


20.2 


16.3 


7.6 



Note: 1. Area ratios as in Figure 2.7 
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3. Discussions 



Looking back at figure 2.1, the dependence of on altitude and on 
Mach number should be explained. We shall do that by using the equations 
described in Appendix A. 

From equations; (1.4.9a), (1.4.12), (1.4.16) 
together with equations: (1.1.1), (1.1.2), (1.1.3) 

one obtains: 





'sp = •'f ”3 


(3.1) 


where; 




(3.2) 


and : 


Cf = X 2 - 


(3.3) 


where; 


X. = k,Mg‘^[X,-l] 

^ ^ ' -0.4 


(3.4) 




X] = ^0) ^ 


(3.5) 


The parameters k^ 


to kg are functions of the various 


area ratios, the 



temperature of air (Tg), the heat capacity ratio of air (y^) and the 
perfect gas constant (R,). These parameters are assumed to be constants 

a 

in discussing the influence of the change in altitude and in Mach number 

on the value of I^p. The altitude dependence is mainly due to change in 

atmospheric pressure (Pg). In the conditions chosen for figure 2.1, the 

0 4 

dependence of I^p on pressure is approximately Pg ' at Mg =3. That means 

that, in the region mentioned, I pressure dependence is mostly due to 

sp 

change of Xg (equations 3.1 & 3.2). From the same equations, the Mach 

number dependence of I^p can also be explained. At high altitude, the 

change of Xg (Mg^‘^^) is very close to the dependence on Mach number 

(Mg ■ ) and therefore I^p is almost constant while changing Mg. On the 

1 35 

other hand, at sea level, the change in Xg(Mg ' ) is smaller than that of 

-2 1 4 

C^(Mg ) and therefore I^p changes with Mg as shown in figure 2.1. 
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Figures 2.2 - 2,5 present similar dependences of the ramjet performance, 
and can well be understood using the same equations. 

Testing these results, together with those presented in Appendix H, 
the design of the ramjet internal area ratios can be completed. The 
results are listed in Table 3.1. 
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Table 3.1: Ramjet Design 



Dimensions : 

External diameter = 5" 



Total length 


= 60.97" 






Total weight 


= 104.7 lb (47.5 kg) 






Area Ratios: 








0 r 


Aj/Ag 


A5/A, 


Ag/Ai 


0.25 


0.47 0.887 0.426 


0.26 


1 



Reference Area : 

= 19.3 sq. in. (124.5 cm^) 

Combustor 

Solid fuel: Hydroxy Terminated Polybutadiene (HTPB). 

Fuel weight: 3 kg 

Fuel density: 971.56 kg/m^ 

Fuel specific impulse (Igp): 400 - 900 sec. 

Booster 

Booster weight: 2 kg 

3 

Booster density: 1650 kg/m 

Booster specific impulse (Igp): 240 sec 

Performance 

Muzzle velocity: 762 m/sec 

Velocity after booster: 863 m/sec 

Thrust Coefficient (C^): 0.2 - 0.4 
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In figures 3,1 - 3.4, the designed ramjet concept is presented. The 
guidance and the control sections as well as the warhead were not redesigned. 
The location of the tailfin is described in figures 3.3 - 3.4 [White, 24]. 

The configuration of the solid fuel ramjet presented here, is in accordance 
with the design of the liquid fuel ramjet done previously by Brown [5]. 

The new design is also in agreement with the Navy’s requirement to be 
compatible with its 5"/54 Mark 46 gun mount, as modified for gun launched 
guided projectiles. 

Figure 2.6 presents a full-range comparison of the SFRJ 5"/54 pro- 
jectile performance with that of the conventional projectile. The 
improvement in oerformance of the ramjet-propelled, guided projectile is 
very significant. The extended range of the ramjet concept can be used in 
a "Surface-to-Surface Mission" (Table 2.1). The main improvement of the 
ramjet concept might be in an "Air-Defense Mission" (figure 2.7). A 
minimum Mach number of Mq = 1.8 was assumed for both the ramjet concept 
and the conventional projectile. At lower altitudes, the projectile is 
limited by burn-out of the fuel. At higher altitudes, the performance 
boundary is due to ramjet unstart. The conventional projectile is limited 
by the Mach number decay to values less than 1.8. It is self-evident 
that the ramjet concept provides significant improvement, and, therefore, 
has significant ASMD Capability. 
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ALL DIMENSIONS IN INCHES 




'^2 



Figure 3.1 Solid Fuel Ramjet Propelled 5"/54 Projectile: Design 



MAXIMUM OUTSIDE DIAMETER 



* 4.96 ► 

ALL DIMENSIONS IN INCHES 




0.076 



SECTION B-B 



Figure 3.2 Solid Fuel Ramjet Propelled 5"/54 
Projectile: Design, Section B-B 
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Figure 



4 



MAXIMUM OUTSIDE DIAMETER 
4.96 

ALL DIMENSIONS IN INCHES 



ALLOWABLE AREA 
FOR EXHAUST 
NOZZLE 




SECTION A-A 



gure 3.4 Aft Body Fin Design: Section A-A [24] 
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4. Conclusion 



Based on the computer model discussed above, the ramjet-propelled, 
guided projectile provides significant improvement over the conventional 
projectile. The fuel specific impulse of the SFRJ is in the order of 600 
±200 sec, depending mainly on the altitude of the projectile. The 
appropriate value of the rocket is only 300 sec. The thrust coefficient 
varies from 0.2 to 0.4, depending on the atmospheric conditions (altitude) 
and on the geometry of the projectile (internal areas). Therefore, the 
ramjet-propelled, guided projectile reaches a range of about 80 km compared 
to range of slightly more than 20 km in the conventional projectile. The 
improvement of the ramjet concept might be in both "Surface-to-Surface 
Mission" and in "Air-Defense Mission". It provides an ASMD weapon which 
is complementary to guided missiles. 

More work is required to design an optical system incorporated into 
the inlet. For satisfactory ramjet performance at the flight Mach numbers, 
the lens must be a conical shape. 
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Appendix A: SOLID FUEL RAMJET: EQUATIONS 



A] . Combustion 

A1 .1 Computation of fuel - regression rate, weight rate of 
burning fuel, and fuel - air ratio 
Define G as the weight flow rate of air per unit through the 
entrance port to the combustor; see figure 1.1. Hence G is given by; 




2 

and has dimensions of Ib/sec.in . The simple form of the regression 
rate of a solid fuel, f, is given by: 

t = aG" (1 .1 .2) 

Where a and n are empirically determined constants. The dimensions for 
f are in/sec. Knowledge of the temperature dependence of the regression 
rate will allow the use of a more accurate model instead of equation 
( 1 . 1 . 2 ). 

The weight of fuel burned per unit of time is as follows: 

^^0313 (1.1.3) 

3 

Where is the density of the fuel in Ib/in , D 3 is the inside dia- 
meter of the fuel grain. Note that D 3 increases as 'the fuel burns. 

The length of fuel grain is L 3 inches. For the case of HTPB, the value 
3 

of 0.0351 Ib/in was taken. In an actual solid fuel ramjet, t varies 
along the grain; t is largest in the region immediately downstream of 
the entrance port of the combustor. However, for the ramjet model 
developed here, the value of t is assumed to be constant along the grain. 
Consequently: 

I t 

D 3 = D 3 Q + 2 j ^ t^dt (1.1 .4) 
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The integral form for the change in grain internal diameter is used 
since f may vary with time. The initial grain inside diameter is D^q 
inches, and the inside area for fuel grain as a function of time can be 
wri tten as : 



A 



3 






nJTT 



'30 





(1.1.5) 



By definition, the fuel - air ratio is: 




The value for is obtained from weight flow through the inlet, and the 
value for is calculated using equation (1.1.3). The total mass flow 
through the nozzle is given by: 

+ f) (1.1.7) 

For HTPB burning in air, the stoichiometric value for f is 0.0728; high 
combustion efficiency is difficult when f is less than 0.025. 

A1 .2 Computation of Combustion Exit Condition 
Combustor exit conditions are specified by four quantities as 
follows: stagnation temperature, T.j.^, Oj^; stagnation pressure, Py^, 

psi; ratio of heat capacities, y^; and gas constant, Ry, in/Oj^. In 

2 

the computer program, the appropriate mks units are used, i.e. Oj^, kg/m , 
m/Oj^, respectively. To determine the exit conditions, two input quan- 
tities are needed. These are fuel - air ratio and stagnation temperature 
at the combustor inlet, T^q. Note that Ty^ = TyQ has been assumed. 

From the thermodynamic data for HTPB burning in air, one determines 
Ty 4 (th), R^. The symbol Ty 4 (th) is a theoretical temperature which 
results from 100% combustion efficiency. Introducing the definition of 
combustion efficiency yields: 
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( 1 . 2 . 1 ) 



'T4 






TO 



+ T 



TO 



As discussed previously, a constant value of rij equal to 0.9 has been 
assumed. 

The value for pj^ is calculated based on one-dimensional, choked 
nozzle flow. A certain value of Pj^ is required to force a certain weight 
flow, Wj, through the nozzle. Assume that remains fixed through the 
nozzle. Define a function of as: 



r = 




2 _ 

Jf + 1 



(y^ + l)/[2(y^ - 1)] 



( 1 . 2 . 2 ) 



Define a characteristic nozzle velocity, 
c*, m/sec: 

. >'3 ''f Tt4 
f 

where g is the acceleration of gravity and has value of 
The required value for Py^ is given by: 



(1.2.3) 
9.807 m/sec^. 



WyC* 

PT4 = Wf 



(1.2.4) 



The decrease of flight stagnation pressure, Pyg, by inlet and combustor 
losses must not be too large. If the inlet and combustor do not provide 
the required Py^, the inlet will unstart and will decrease. 
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AT .3 Computation of Nozzle Exit Conditions 

The relation between the area ratios and the Mach number is well 
known by the formula: 






(Yf + 1 )/2 ) ^ ^ - 1 )] 



1 + <f 



r, - 1 



(1.3.1) 



Ag, Ag are the areas at the throat and at the exit of the nozzle, re- 
spectively. Knowing the exit Mach number (Mg) can be calculated 
for any nozzle area ratio (Ag/Ag). This indirect calculation is done 
in subroutine CALCM, using Newton - Raphson's iteration routine. 

The total pressure at the exit of the nozzle (Pjg) is defined by: 

Pt6 " Pt 4 ^n (1.3.2) 

where the total pressure at the exit of the combustor (Pj^) was calculated 
previously (e.g. 1.2.4). 

Knowing the total pressure at the exit of the nozzle (Pjg) and the 
Mach number at this point (Mg), the exit pressure (pg) can be calculated: 

/ Y -1 \TYf/(Yf-D 

^6 " Pt6 2 ~ ^6 j (1.3.3) 

A1 . 4 Computation of Thrust and Thrust Coefficient 
Al.4.1 Thrust Coefficient (C^) 

The thrust of the engine is the net rate of change in 
momentum at a steady state condition, and is given by: 

F = PgAg + itlgUg - PgAg - iflUg - Pg(A^ - Ag ) + Og(A^-Ag) 

^ ^ 6^6 '^ 6*^6 ■ 0 - 4 . 1 ) 
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where Uq, Ug, itig, iTig are the velocities and mass flow at the inlet 
entrance and at the nozzle exit, respectively. 

From the continuity equation, the following relation arrives; 

tfiU = pU^A (1.4.2) 

Substituting for the density (p) from the perfect gas equation of state: 

(1.4.3) 

gives: rtiU = ^ U“A (1 .4.4) 

From the definition of Mach number and speed of sound: 



P " h 

itiU - pj- U A 



= M^a^ ; a^ = yRT (1 .4.5) 



Therefore; 



ifiU = gyM^yRTA 
ttiU = oM^yA 



(1.4.7) 



Substituting (1.4.7) into (1.4.1) gives: 

The thrust coefficient is defined ; 

'0^ 






where: 






(1.4.8) 

(1.4.9) 
(1 .4.10) 



Combining equation (1.4.9) and 1.4.10) gives: 

F 



f y 



4 PO^O V 



(1.4.11) 



Substituting for the thrust from equation (1.4.8) turns equation (1.4.11) 
i nto : 



^ ^^6^^r r PtS^Pq ^ 2x ' 
^ y^Mg^ LptS^Ps ® 



2A 

T 



0 (1.4.12) 
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Al.4.2 Pressure Losses 



The pressure ratios in the above formula: 

Pt6^PT0 

Pje/Pe 

can be substituted by a function of pressure losses across the ramjet: 

(1.4.13) 

Pq Pt4 Pt3 PJ2 PjO Pq. 



We define the pressure losses as follows: 

PT6 



T4 



= iT^ = Nozzle losses 



PT4 

- — = TT. = Rayleigh flow losses 
Pj3 ^ 

Pt3 

- — = = Combustor expansion losses 

PT2 ® 



Pj2 

’ Inlet losses = (conical wave loss )*(boundary layer loss)* 

Pjo ^ 



■( normal shock loss )*(subsonic diffuser recovery) = 



Therefore: 



Pq 



T6 _______ _ _ „ PtO 

" "^n ^h ^e ^0 ^NS Pq 



n(iri) Pj^ 

_ i = losses Jpq 



(1.4.14) 



Finally, substituting for total pressure ratios (pjq/Pq; pyg/pg) gives: 

-,Y,/(Y.-1) 

Y.-l 



PtS^Pq 

PtG'^Po 



n(TT.) 

_ i = losses -* 



1 + At41 ^ 



1 + 






2^ '6 



(1.4.15) 



This relation can be used in equation (1.4.12) which calculates the 
thrust coefficient of the system. 



A1.4.3 Computation of Thrust 

The thrust can easily be calculated from the thrust coefficient, 
using eq. (1.4.9): 

F = qQ (1.4.9a) 

The dimensions of F are Newtons (after multiplying eq. (1.4.9a) by 
the acceleration of gravity, g). The fuel specific impulse, I , in 
N/ kg/ sec, is defined by: 

Igp = F/W^ (1.4.16) 

The specific fuel consumption, SFC, in kg/hour/N is given by: 

SFC = 3600/Igp (1.4.17) 

Ramjet performance is specified in terms of the performance parameters 
C^, F, Igp and SFC. 
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A2. Check for Choked Nozzle 



The total pressure at the throat of the nozzle is given by: 

Pt5 ^ Pt 4 ^^n ^ 

Again, Pj^, is the total pressure at the exit of the combustor, and is 

calculated by eq. (1.2.4). is the nozzle loss. It follows that; 

/ 2 ( 2 . 2 ) 



Pc = D 



T5 



The pressure at the throat of the nozzle (p^) should be equal or greater 
than the atmosperic pressure (Pq)- 

Ps> Pq 

When inequality (2.3) is satisfied, the nozzle is choked. 

A3. Heat Losses at the Combustor 
A3. 1 Mach Number 

a. Conti nui ty 

At the combustion chamber, fuel is added and the continuity 
equation is: 

P 3 U 3 A 3 (1 + f/a) = P 4 U 4 A^ (3.1.1) 

where A 3 and A^ refer to the entrance and to the exit of the combustor, 
respectively. By assuming that A 3 = A^, the continuity equation 
(eq. 3.1.1) can be written as follows: 

^ (1 + f/a) = iji (3.1 .2) 

P 4 ^3 



Replacing the velocities (U 3 and U^) by the appropriate Mach numbers 
(eq. 1.4.5), turns eq. (3.1.3) into; 




^ f/a) 



^3 '^a ''a 



(3.1.3) 
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Where R , R^: are the gas constants of air and of the combustion products, 

d T 

respectively. Hence: 




(3.1.4) 



Applying the conservation law of momentum to the discussed 



problem: 

^3 *^3 *^3^ ^ P 4 P 4 ^ 4 ^ (3.1.5) 

From the definition of Mach number (eq. 1.4.5) and the perfect gas 
equation of state (eq. 1.4.3): 

P.j U.2 = Yi P^- M.2 (3.1.6) 

Substituting equation (3.1.6) in equation (3.1.5): 

^ ^ ^''a ^3^ (3.1 .7) 

P3 1 + Yf 



Again, from the perfect gas equation of state: 




Substitution of (3.1.8) into (3.1.7): 

\P3/ V4/ \1 + YfM^ 



(3.1.9) 



Substituting (eq. 3.1.9) turns equation (3.1.4) into: 




\1 + YfM^ 




1 

1 + f/a 



(3.1.10) 
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Replacing the temperatures and by the appropriate total temperatures : 

’X' *- 






Yi-1 2 



2 “i 

and assuming that the change in total temperature up to the entrance of 
the combustor, is negligible (Tjq = Tj^)^ one obtains: 



'M, 



'T4 



M / \ T 

oy Vto 



= /]_lVV 
V * 

c. Solution 



1 + 



Yf-1 






2 "4 



1 + 



_i — M 

3 



. A 



^^a'a 



1 



[1 + f/a 



(3.1.11) 



The Mach number at the exit of the combustor (^'' 4 ) can be 
solved knowing the conditions at the throat of the nozzle: 



Yf-1 



7tr-M4' 



1 + 



T-^4 



(3.1 .12) 



The computation is again indirect, using subroutine CALCM. Knowing 
M^, Equation (3.1.11) is used to compute M^. The solution is received 
by iteration. As first approximation, M^ can be solved from the 
following equation: 

^ " ( 3.1 .lla) 






where : 



B = 



1 + YrM.M /T^, \- (1 + f/a) 



TA 1 / 'T4 

M. } [ 



(3.1.13) 



y^R 



f'f 



consequentTy: -y 

+B-^B^-4Ya 
'^3N - 2— 



(3.1 .1/’-) 



The subscript M shows that the computation was made from the nozzle 
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direction. Mow, B is changed to be: 



B = B 



2'- 

1 + M32 






Y ^-1 



2 ^4 



(3.1 .15) 



This expression is consistent with equation (3.1.11). Substituting 
back into (eq. 3.1.14) gives an improved value for M^. This procedure 
can be repeated several times, but it was found that even after two 
interations, the value received for ^3N accurate enough, due to the 
small change in B resulting from (eq. 3.1.15). 

A3 . 2 . Total Pressure 
By definition: 



Pi 

Consequently: 






1 



?T 3 = ^ ) 


1 + -^3^ 


v<v> 


Pj 4 P4 j 




Yf/(Yf+l ) 



Substituting for (p^/p^) from equation (3.1.7) results: 



Pt3N " Pt4 



(3.2.1 ) 



(3.2.2) 



(3.2.3) 



Where the subscript N was defined previously. 
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A4. Computation of Maqh Number and of Total Pressure at the Various 



Stations of the Inlet 

A4.1 Initial Conditions 

In the previous sections (AT - A3) the pressure conditions at 
the combustor and at the nozzle region were calculated. Here, the 
pressure conditions at the inlet will be calculated independently. 
Knowing the total pressure conditions at the various stations of the 
inlet will allow the check of whether the inlet can supply the amount 
of air needed by the combustor. As will be seen afterwards, this check 
will also allow to specify the location of the normal shock wave at the 
ini et. 

Assuming that the static pressure (Pq), the static temperature 
(Tg) and the flight Mach number (Mg) are known from the trajectory part 
of the program, the total pressure and the total temperature can be 
calculated: 



Y -1 o Y,/(y.-1 ) 


PjO " Po ~2 ^0 ^ 


(4.1.1) 


Yo = ■^0 * 1 


(4.1.2) 


The weight flow through the inlet is given by: 


'^a " PO ''^O ^0 (% 


(4.1 .3) 



Usually, when flight Mach number is equal or greater than the inlet 
design Mach number, the value for Ag/Ag is unity. But for flight Mach 
number less than design Mach number, Ag/Ag becomes less than 1.0. A 
value of 0.9 was selected as a constant value for Ag/Ag. The additive 
drag due to Ag/Ag < 1 was ignored. 
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A4.2 Conical Shock Wave Loss 



In this section, the conical shock wave loss will be computed; 
the calculation results include the total pressure, Mach number and 
area behind the conical wave (Pji^j respectively). 

A4 . 2 . 1 Pressure 

The oressure coefficient can be defined as follows: 



C 



P 



PlC'°0 



(4.2.1) 



For a cone, the pressure coefficient (Cp) can approximately be for- 
mulated as: 



C 



P 




(4.2.2) 



where a is the cone half angle. The difference in pressure on surface 
and behind shock wave is ignored in this model. 

Knowing Cp, the pressure ratio (Pk^/Pq) can be calculated (4.2.1) 



Plf fa 2 

— = 1 + C ^ Mn 
Pn p 2 0 



(4.2.1a) 



The same pressure ratio, is also related to the Mach number, normal to 
the conical shock wave (M^^). Using (eq. 2.48a) in ref. [16], one can 
get: 

(4.2.3) 



.!Ia (V-’> 



2t. 



Knowing (p-jq/Pq) 'from equation (4.2.1a), M^ can be calculated from 



equation (4.2.3) ; 



M = 
n 



1 + 



Pq 



- 1 



2fa J 






(4.2.3a) 
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After computing the pressure ratio due to conical shock wave (equation 
4.2.1a) and the Mach number normal to the cone (equation 4.2.3a), one 
can use equation (2.54) in reference [16] to compute the total pressure 
ratio at the conical shock wave: 



TIC 

^TO 



[’• 



2y. 



] 






From equation (4.2.4) one obtains: 



[ (V’lT 1 

n 



V('^a-l) 



(4.2.4) 



PflC “ °T0 (4.2.4a) 

A4 . 2 . 2 Mach Number Downstream of Conical Shock Wave 
Using figures 4.1 and 4.2, the wave angle (0) can be defined 



as : 



6 



arcsin 




(4.2.5) 



On the other hand, the deflection angle (6) is defined from equation 
(4.10) in reference [16]: 



9 = arctan 



2cot (3) — 2 ~ 



(M ^-1) 
' n ' 



(4.2.6) 



Mg {y^+cos(23)}+2 - 
The Mach number behind the conical shock wave (M^g) niay, therefore, be 
obtained using equation (4.7) in reference [16]: 



'1C 



Y = ” ^ 0 



sin^(3-3 ) 






'^a-1 



(4.2.7) 



The relation between the deflection angle (9) and the wave angle (3) for 
various Mach numbers are shown in figure 4.3 (reproduced from figure 4.2 
in reference [16] ) . 
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Figure A4.3 Oblique Shock Solutions [16] 



A4.2.3 Area Ratio 



Using figure 4.2, the area ratio behind the conical shock 
wave can be obtained: 



^1 C _ 1 -tan(e ) 
Aq tan(S) 



CO s ( 9 ) 



:4.2.8) 



Equation (4.2.8) gives the area normal to the flow at the inlet lip. 
A4.3 Boundary Layer Loss 
A4 . 3 . 1 Mach Number 

The area ratio is related to the appropriate Mach number by 
by the formula: 








|(y3+1)/[2(y3-1)] 



(4.3.1) 



Where A.j is the area at the throat of the inlet. Similarly, 



^1 



((y+D/2 ) 

- ^ 11 ! 

1 + 2 ^1 



(y +1)/[2(y.-1)] 



(4.3.2) 



Where A.ji and M.j.| relates to the area and to the Mach number ahead 
of the normal shock wave. Dividing. these two equations gives: 
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Knowing equation (4.3.3) can be used to calculate 

indirectly, by subroutine CALM, which was mentioned previously. 
should be supersonic > 1) to prevent unstart conditions. 

A4.3.2 Pressure 

The total pressure in front of the normal shock wave is 

received from the connection: 



Pill ■ PtIC ^d' (4.3.4) 

Where , the boundary layer loss is assumed to be 0.93. 

A4.4 Normal Shock Loss 

The Mach number, behind the normal shock wave is defined as: 



12 






r 



"if * -^a-1 



(4.4.1) 



The total pressure behind the normal shock wave is defined as: 



Pt12 



1 

° PtII ' 


K'.„' 


) ■^a^l 


1 


ll^ 2 ^ 11 ' 






2y 



1 M 



3 M 2 

M, 






/(Y^-l) 






^a"l 



(4.4.2) 



A4. 5 Subsonic Diffuser Recovery 
Similar to equation (4.3.3) one can obtain: 



Ar 



Mi 2 \ 1 + 2 '-'2 



Y.-l 2 



M^ 



hz ”2 I 1 t Vi „ 2 



(4.5.1) 



2 '12 ; 

Knowing y^, Aj^2» ^ 2 ’ inlet 

can be computed using subroutine CALCM. The total pressure at this station 
is defined as: 

Pj2 ~ Pfl2 (4.5.2) 
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Where the subsonic diffuser recovery (irp") is assumed to be 0.93. 

A4 . 6 Expansion Loss 

A4 . 6 . 1 Mach Number 

On their way to the combustor, the gases coming from the inlet 
expand at station 3; a sudden change in area from ^2 to A^ occurs. The 
sudden change in area acts as a flameholder by creating a hot recirculation 
region. In this section, the loss in total pressure due to this expansion 
is calculated. 

From the continuity equation: 

P 2 U 2 A 2 = P 3 U 3 A 3 (1.4.2a) 

From perfect gas equation of state: 




and from the definition of Mach number and speed of sound: 
M = U/a; a = v/W (1.4.5) 

Equation (1.4.2a) turns, therefore, into the form: 



RT 2 ^2 ’'^a'^^2 ^2 " RT 3 ^3 ^3 (4.6.1) 

For sudden expansion of an incompressible fluid, the change in static 
pressure going from small area A 2 to large area A 2 is given by [16, 25]: 



P3 ~ P2 
^2 



2^23(1‘A23) 



(4.6.2) 



Where A 23 = A 2 /A 3 . Also, for the incompressible case, the stagnation 
pressure ratio is given by [16, 25]: 




q2/P2 

■ T+q^ 



2 



(1 - A23) 



(4.6.3) 



According to equation (4.6.3) as M 2 decreases Pj 3 approaches p-|- 2 - The 
model for sudden expansion with compressible flow is much more complicated. 
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It was assumed that static pressure is constant in the expansion. This 



assumption, which is reasonable for low values of 1 ^ 2 , is a conservative 
one, i.e, p_,/p_„ is lower. 

r '^a'\ 2" 

Substituting: 2 ~ ^ ^ — 2~^^2 



- 'a-^ 2' 

1 + 



and 



P 2 = P 3 , turns equation (4.6.1) to: 




(4.6.4) 



Solving equation (4.6.4) for M^: 



/ \J 1 +4a6-l 
31 ~ sj 2a 



(4.6.5) 



Where: 





(1 + aM2^) 



(4.6.6) 



The subscript I shows that the computation was made from the inlet 
direction. 



A4 . 6 . 2 Pressure 

By assuming again that P 2 = p^, 

Pt3 ^ 

Pf 2 PJ2/P2 

and therefore: 



(4.6.7) 



Pt3I " Pt2 









(4.6.8) 



Where subscript I is as defined previously. 
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A4.7 Location of Normal Shock Wave 



The main problem in computing the Mach numbers and the total pres- 
sures at the inlet arises from the fact that the exact location of the 
normal shock wave is not known, and must be found. The way of solving 
this problem is as follows: 

First, solve for two extreme conditions by assuming that the normal 
shock wave is located at the throat and at the exit of the inlet, respec- 
tively. After knowing the lower and the upper values for M^j, P-[- 3 j» 
iteration can be made to find the exact location of the normal shock 
wave. The criteria for this iteration is matching of values for M^, 
p.^.^ from both the inlet and the nozzle directions, i.e.: 

^31 " Pt3I " PT3N (4-7.1) 
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APPENDIX B: TRAJECTORY EQUATIONS 



B1 . Atmospheric Functions 

Best fit curves were calculated for basic atmospheric functions, 
pressure, density, temperature and viscosity of air as a function of 
the flight altitude. The basic formula which was used for this process 
is as follows: 

F = A exp(-B X 10"® h^) 0 •'' ) 

Where h is the altitude in meters, and A,B,C are numerical parameters. 

The appropriate atmospheric functions are as follows: 





Pg = 1.03322 X lo'^ exp (-59.148 X 10’® h^ 


(1.2) 




Pg = 1.224845 exp(-29.0144 X 10’^ h 
Tg = 288.16 exp(-13.232 X 10'^ h "I 


(1.3) 


When: 


h = 0 - 11,000m 


(1.4) 




Tg = 217.24°, when h = 11 ,000-32,000m 

Pg = 1.793 X 10'^ exp(-45.1374 X 10’^ ^^0.8984^ 


(1.4a) 


When: 


h = 0-11 ,000m 


(1.5) 




Ug = 1.41724 X 10"^, when h > 11,000 m 


(1.5a) 


2 

In these formulae, the pressure (pg) has dimensions of kg/m , the 


2 

density (pg) is given in kg/m , the temperature (Tg) is 


given in 



2 

and the viscosity (pg) is given in kg/(m-sec) (or: N-sec/m ). 

B2. Drag 

B2. 1 Cowl Drag Coefficient 

2.1.1 It was found that the cowl drag coefficient has a strong 
influence on the results. Therefore, a new model for this cowl drag 
coefficient was developed. The model was based on a theoretical 
development done previously by Prof. T. H. Gawain [9]. The main 
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difference between this development and the classical theory, is that 
the boundary conditions are applied at the body surface rather than 
along the axis. 



2.1.2 The model, which originally was developed for simple cases 
(cones, etc.), was modified to fit the shape of the projectile, illustrated 
in Figure 1.1. 

2.1.3 The modified program is listed in Appendix G. In the combined 
program (TRAJET) an interpolation procedure was used as a subroutine in 
order to simplify the calculation process. 

B2.2 Base Drag 

After checking the influence of the nozzle exit area (Ag) on the 
performance, it was decided to allow Ag to reach the maximum value 
possible (A^), in order to reduce base drag. In this case, the base drag 
is negligible. 

B2.3 Skin Drag Coefficient 

The Reynolds number is well known to be: 

_ pq^o^ 

■ li 

The transition Reynolds number of 
Re* = 2*10^ 



(2.3.1) 



is usually taken as criterion for transition between laminar flow 
(lower values) and turbulent flow (higher values). The incompressible 
laminar skin friction coefficient is related to the Reynolds number 
as follows: _ 



‘^DS,L " ^ •328/>'Rej^ 



(2.3.2) 



Where the subscripts DS.L stand for skin drag coefficient, and laminar 
flow, respectively. On the other hand, the turbulent skin friction 
coefficient (C^g j) is calculated indirectly from the formula: 
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(2.3.3) 



The computation is done in subroutine CALDC, which works in a similar 
way to subroutine CALCM which has been described earlier concerning the 
calculation of various Mach numbers. 

B2.4 Wing and Fin Drag Coefficients 

The wings and the fins of the projectile, also contribute to 
drag. Basically, each of these drag coefficients contains two parts: 

-Wing/ fin wave drag. 

-Wing/fin- frtction drag. 

2.4.1 A psuedo 3-dimensional model was chosen to simulate the wave 
coefficient. The basic formulae used in this calculation are: 






Y.-l 






Y,+l o -Y;,/(Ya+l) 
^TO ^ ^ ^ T" ^0 ^ 



The effective span could be taken as (see fig. B2.1) 
b' 



h ^ 

b - 2 



Substituting: tany = l/c , sinu = 1/M gives: 



b' = b - 



2 ^\r2T 

Mq -1 



The drag coefficient would, therefore, be: 

2 / '^01 ^'^TO - PQ2'^Pt0 



tb' 



'DWW ,^^^2 \ p7/p 



0^ 'TO 



Pq/Pto 



(2.4.2) 



(2.4.3) 



(2.4.4) 



(2.4.5) 



The way this simulation uses the above equations could easily be 
understood when looking at the formulae together with the flow chart of 
the appropriate subroutine (Appendix C). Interference drag is ignored. 
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Figure B2.1 Schematic View of an Wing/Fin 
b=span, c=chord, t=thickness 
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2.4.2 The friction drag coefficient was calculated using the 
existing model for skin drag coefficient. (See Section B2.3) 

B2.5 Calculation of Drag 

Define the dynamic pressure (9) as follows: 

Q ~ h Pq^o (2.5.1 ) 

7 2 

When q is in units of kg/(m.sec ), (or: M/m ). Also, define the following 

geometrical units: 

A = ttR^; S = 2ttRL (2.5.2) 

r r 

Where R,L are the radius and the length of the projectile, respectively. 
Similiarly, for the wings or the fins: 

= nbc ■ (2.5.3) 

Where n is the total number of wings/fins (a value of 8 was taken for n), 
and b,c are the span and the chord of the wing/fin (see fig. B2.1) 
Consequently, the drag (D) is given by: 

° ^ j% ‘'DM ^p *'DS ^WW ^‘'DWW ‘'DWF^ [ (2.5.4) 

B2.6 Drag Coefficient of a Conventional Projectile Without 
Propul Sion 

The program has an option to calculate also the trajectory of a 
projectile without a propulsion. The projectile is a conventional round. 
The formulae used to claculate the drag coefficients in this case are: 

2.6.1 Nose Drag 

Cpj^ = (0.083 + 0.096/Mg^)(a/10)^-^^ 

Where a is the cone half angle. 

2.6.2 Base Drag 

Cgg = (0.6837 - 0.3165 M + 0.0525 M^)(2/tt) 



54 



2.6.3 Skin Drag 

Skin drag is calculated as discussed in section B2.3. 

2.6.4 Drag 

D = Drag = q (Cp^ + C^r) + C 



P OS I 



B3. Booster 



The projectile has an initial muzzle velocity of 2500 ft/sec. Part 
of the combustor volume can be used as a booster to accelerate the 
projectile even more so that starting the ramjet will be easier. 



Define exhaust velocity (0^) as: 



U 



I, 



(3.1) 



e ^sp,B y 

Where I n is the specific impulse of the booster's fuel (in sec) and g 

Sp 5 D 

is the acceleration of gravity (in m/sec). 

From Newton's law [6, p. 323] 



dU 

"B "e ■ '"'p - dt 



= (m„ - rtiot) 



(3.2) 



Where m and rtio are the mass of the projectile and the mass flow of 
P D 

the booster respectively. 

dt 



dU = itin — . . 

B e m„ - itiot 

P B 



(3.3) 



'Consequently: 



m ~ rfinX 

AU = U(t) - U(0) = -Uln -2- 



m. 



(3.4) 



Where x is the booster burn time. Hence, 

m 



iU = U 



® '"p - <^B 



(3.5) 



AU is the change in initial velocity due to the booster, where m„ is the 

D 

mass of the booster. 
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B4 . Dynamics 

The flat earth trajectory with drag and thrust is well known. 
The differential equations of motion are: 



-g + (F-0) sin9/m^ 



= 

dt2 
d^x 

nir, — p = (F-D) COS0 



(4.1) 

(4.2) 



where y is the altitude, t is the time, g is the acceleration of 
gravity, F is the thrust, 0 is the drag, 0 is the elevation angle, 
mp is the projectile mass. 

For numerical solution of equations (4.1) and (4.2) a finite difference 
form can be used as follows: 



Xj+2 Ar/mp + 2Xj^^ - Xj 

yj+2 = C-g + C"-D) sin0/mp] At^ + 2y^.^., -y^. 
The values for j = 1 are from initial conditions, i.e. 

= 0, Yj =0 

X 2 = Uq * COS0 * At, Y 2 = LIq * sin0 * At 



0 = arctan 



'.i>2 - y.i>i 









(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 



Having calculated x^^^ yj+2 (®g* 4.3, and 4.4), one obtained new 
values for trajectory parameters using (eq. 4.7). Also thrust, drag 
and projectile mass are updated for j+2. 
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Appendix C : Flow Chart of the Computer Program (TRAJET ) 



C 1. Main Program 




ATM 

( 100 ) 



BOOS 



INIT 

(155) 
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PRILOO 

(CMB) 
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C2 . Command Subroutines 




W 



W 
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62 





I 
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C 3. 



Individual Subroutines 
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68 
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INTER 
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notes : 



1 . 



7 




compute : 






M (new)= 




T 



or 



1 = known mach number 

^ works in a similiar way, 



number , 
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respectively . 
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C3.8 INLET 
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DIFFUS 




Return 



EXPAN 



, A^ 



compute : 
M 

^31 



compute ; 
p 

T3I 



Return 
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C3.12 



TRAJ 



INTER 



compute : 
^DN 
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APPENDIX 0: 



PROGRAM TRAJET: LISTING 
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>ooon 



Fuei TRAJET 



FORTRAN A NAVAL POSTGRADUATE SCHOOL 



cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

PROGRAM TRAJET C 

c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 



THIS IS A PROGRAM FOR SOLID FUEL RAMJ ET^-TRA JECTOR Y 



CC 



CCMHON/GfO/AR6F,AO»Al, a 2 »A30,A 
f f AIA0,A2A0«A6A5, A3AR 
CC«M0N/A IR/G A, GAl ,GA2,GA 12 #GA3 
CC^'Mrf4/ FCFL/ RHCF , FIAT, A, N 
CC /I C S S/ P I Cl . P I 02 # P I a I , P I B 

COMMON /dLR/W A, UF ,'«T,F,ROQT ,CST 
CCMMCM/r,f I/RF.GF.GFl,(;f2,GFl2» 
C0PMCN/NCZ/'<6,PT6,P6 
common /OFS/CP. TFPUST .1 SP.S‘=C 
Cav^ON/ChC/?T5 .P‘5 
CC^mon/ I M/p tC.TTO 
CQPMON/CCN/PTIC.MIC .AlC, ALFA 
CCMMf.N/THFn/Ml l.PTll ,M2l .PT21. 
CG^Hf’N/NS/fn 2. PT12.M22.PT2 2.HS 
CO^MnN/OIF/M2.PT2.P2 
COMMON /CHF/M3N ,M3l . PT3 .PT3M.MA 
CCMMCN/TRA/P I . A. OR AG .RPP . LPP 
.IL.ILO.Y.YO. TGB.Xl.Yl.X2.T2.X3 
CCmmOM/OPG/TET P.CCN.CFT.COWW.C 
RFAL MO ./»IC, -Ml l.M 12 ,M2 .M3N.M31 
RFAL L 3 .N . ISP. M2l .M22.MSI ,MS2 
RFAL mua.LPR 

call oatta 

IFU 6 .GT.ARFFIG 0 TO 120 
Y* 0 . 

CALL ATM 

IFI IRAM.EO.OJGO TO 100 

ITRA=-l 

THF UST* 0 . 

Ot-LU^C. 
r,r TO 140 
100 thpust=i. 

CALI BDCS 
140 U=tO^OFLU 

MO=U/SQRT(GA*GRAV*R*TOI 

ROCT^O. 

061 T= 0 . 

WPRO-WPR 
ACARsAO/AREF 
A 5 AR = A 5 /APEF 
A 6 A 5 -A 6 /A 5 
ACC=AU 

aiao*ai/ao / 

A 2 A 0 =A 2 /AO 
IPRIM= I 

IFUTRA.GT.O) go to 155 
PRINT 80 C 

800 FCPMATI in, 30 X, ‘SOLID FUEL RAM 
, 5 X ,» AO/AP* , 6 X, *A 5 /AR*, 6 X. *TOB* 
, 8 X, » X< MAX ) ». 5 X, » Y(HAXJ » . 5 X. *TE 
THIS ^‘APT IS FOR AO , A 5 -C HANGES 
IFN 0 * 4 l 

IF( IRAH.eO.n I£N 0«26 
00 III I A 5 * 26 . IFNO 
IFI U 5 .EC.DCn TC llO 
A 5 *( IA 5 * U*O.Ol*AREF 
110 A 6 =A 6 A 5 »A 5 

IF(A 6 .GT.AREF)G 0 to 120 
AO=AOO 

OC 114 IAO= 26 ,IENO 

iLro=o 

Y* 0 , 

(FIl R AM.GT. I ) IPAMal 
IFI IRAM. 60 .UGD TO 145 
THPUSTsl . 



3.A5.A6. L3.A0AR, A5AR, ILOO.I PR 
,RHOA, TO.PO.UO.MO.GRAV 
2, PIN 

AR, PT4, TT4 

r,F3 



HSl ,PTSl 
2.PTS2 

.AS 

,WPR,U.WB.DELU.OELT,TT,TETAt IRAM. 

, Y3,wPRB,R,T0Ra. ITRA. TETAOO. XOB. YOB 
OWF .APR.SPR .SHW.O .XMO 
.M4.M5 «M6 



JET C TRAI eCTORY» ,//43X. ‘SUMMARY* ,// 
,ax. »xoB*. ax. • yob* . ax. • tof* , 

TP* ,7X.*T6TA*./I 
ONLY 



TRAOOOlO 
TRA00020 
TRA00030 
TRA00040 
TRA00050 
TRA00060 
TRA00070 
TRAO 3080 
TRA0009O 
TRAOOlOO 
TRAO JllO 
TR ^03120 
TRAO J130 
TRA00140 
TRAO 0150 
TRAO 01 60 
TRAOOl 70 
^RAO 0180 
TRAOOloO 
TRAO 0200 
TRAO J2I0 
TRA00220 
TRA00230 
TRAO 3240 
TRA0:250 
TRAOOPOO 
TRA00270 
TR AO 02 80 
TRAO 3200 
TRAO 3300 
TRAO 3 no 
TRA00320 
TRA0C330 
TPA00340 
TRA00350 
TRA00360 
TRA00370 
TPA00380 
TRA00300 
TRA00400 
TRA00410 
TRA00420 
TRA00430 
TRAO 0440 
TPA00450 
TRA00460 
TOA00470 
TR A00480 
TR A00400 
TRAO 0500 
TR A00510 
TRA00520 
TRA00530 
TPAO 3540 
TRA00550 
TRA0C560 
TRAO 0570 
TRAO 05 30 
TRA00550 
TRA00600 
TRA00610 
TPA0 3620 
TRA00630 
TRA00640 
•“RAODdSO 
TRAO 3660 
TR A00670 
TRA00680 
TRA00600 
TRA00700 
TRA'3 0 710 
TRAOC720 
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FlUE: TRAJET FORTRAN A NAVAL POSTGRAOUAH SCHOOL 



CALL ATM 

call boos 

145 U-UO^nCLL 

HC=U/SOR 7 (GA*GRAV*R*TO ) 
PCCT=0, 

oelt=o. 

WPRsWPR 0 
IT0A=-1 
TI =0. 

TCR9=0. 

XC8=0. 

tc°=o. 

IF( UO .EC.UGO TO 113 
AC = ( lAO- l)*0 •Ol*APEF 
Al-Al AO*AO 
42 = A2 AO 
113 CCNTInUF 

ACAP = AC/4PEF 
A5AR-A 5 / AP EF 

155 IFl Ipam.cq, DGC JQ 150 
CALL IN IT 
CAU «MP^ 

IFmnr.GT.l )GC TO 120 

IF ( ^hPilS T.LT •O.Ol JGO TO 150 

CALL NOZZ 

IF( ILfJ'KCT.l JGO TO 120 
CALL CHpKC 

IF( nno.GT.l )GC TO 120 
CALL hFAT 

IFULCU.GT.I JGC to 120 
CALL I'^'LFT 

IFdL^O.GT.l )G0 TO 120 
CALL rxPAN 
CALL CHECK 

IF( I LCQ.LE. DGO TO 130 
THW US Tr(J, 

TO=C^TI 

xnn = x3 
vrB = Y3 
ILUO^C 
GO to 150 
130 CALL CORVAL 

IF( ILUO.GT.l JGO TO 120 
CALL PFSLL 

IF( ILno.GT.l )G0 TO 120 
150 CALL TPAJ 

IF( ILHO.GT.I )G0 to 120 
I F( IPP.lT.O) go to 152 
IFtTH- UST.LE .0.01 )G0 TO 150 
IFUTOa.lT.O J GO TO 15 5 



152 



CALL PP ILCOI IPRIN) 

IF( ThkuST.LE .0.01 )G0 TO 150 



GO TO 155 
120 CCNTTMUE 

IFHTPA.GT.O) go to 157 
114 CONTIfiUF 
III CONTINUE 

157 PRINT 2 2 2,ILnO,AOAR,A5AP,AlAO*A2AO*A3AR. 

, TETP.T*=TaOO, [PR, ITR A.X.MO*? IDl. P!02, P IN 
222 format ( lhl,24( •♦*). 13. 24( •<'M*/5X,» INPUT OAT A; * , /2X . 5FL0. 3 . 
/2X.2R6.1.2I3.4F6.3) 

STOP 

ENC 

SLERruTUF INI T 

CCyMf'N/GEO/AREF.AO. Al, A2 ,A3O.A3.A5,A6.L3.AOAR,A5AR,IL0O.IPR 

CCmmpn/a IP/GA.GAI.GA2. GA12.GA3.RH0A. T 0 .P 0 . UO. M 0 i GR A V 

CCPMLN/FLrL/RHCF.ET AT, A.N 

CCVMPN/lCSS/PICI ,PI02,PIRl.P!R2.PIN 

CG^Mf N/HLR/W A, .*F , F , root ,CST AR, PT4 .TT4 

CC‘^‘^rN/GFT/PF,GF.GF I.GP2.GF12.GF3 

COy^LN /NCZ/M6. PTh ,P6 



TRA00730 
TRA00740 
TPA00750 
TPA00740 
TP \00770 
TP A00780 
TRA00700 
TRA00600 
TR 400R 10 
TR AO 0H20 
^P AOOblO 
TP A0ua40 
TPAOOrtSO 
TR AO GH50 
to AO 3H70 
TR AO OrtRO 
'^R AO 08NO 
TR400900 
TO A009I 0 
TP AO 3020 



TPA00930 
TRA0C^40 
TRAO 0950 
TP A00960 
TR AO O') 70 
r Q A009HO 
TP AO 0990 
TRAO 1010 
TRAO lOlO 
TR AOIO20 
TRAO 1030 
TRAO 1040 
TRAO 1050 
TRAO 1060 
TRAO 10 70 
"^P AO lORO 
TPAO 1090 
TRAO I 100 
TPAOlllO 
TRAO 1120 
TRA0U30 
TRAO 1140 
TRAO I 150 
TRAO 1160 
TRAO 1170 
TRAO 1190 
TPAO 1190 
■^P AO 1200 
TRAO 1210 
TP AO 1270 
TRAO 1230 
TR A 1 1240 
TRAO 1250 
TRAO 1260 
TPAO 1270 
TP AO 1290 
TP AO 1290 
TPAO I 300 
TR AOl 310 
IRAQ I 320 



TRAO 1330 
TRAO 1 140 
T°AO I 350 
TRAO 1360 
TPAO 1370 
TPAO 1390 
TRA01390 
TR AO 1400 
TP AO 1410 
TR AO 1420 
TPAO 1430 
TRAO 1440 
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FILES TRAJET FORTRAN A NAVAL POSTGRADUATE SCHOOL 



COHMCN/RES/CF,THeuST,lSP.SFC ! IRAQ 1450 

CCV*-CN/CH»^/0T5 *P 5 TRA01460 

CGP^TN/ I M /P TO.TTO TPA0I47O 

CCN40N/CCM/P TlC.?^lC,AlCfALFA TRAO I4rt0 

COyMC w/ThRP/Hil,pTU,.M2UPT2l* '^Sl ,PTSl TRA0 14P0 

CC'^MPN/NS/'^l 2, PT I2,M22,PT2 2.HS2,PTS2 TP AO 1500 

/n IF /M 2 , PT2 ,P2 ^RAOISIO 

CCf^MO‘i/CPE/M JN ,M3 I , PT3 .PT3M, Ml* ,AS IRAQ 1520 

RFal mo VI 1. Vl2 .M2 ,M3^i,,M3I ,M6 TR A01530 

REAL l 3 .N. ISP, V2i.M22.HSl ,MS2 TRA01540 

GAl=:(r,A f I. )/2, TRA01550 

GA2=(r,A - U )/ 2 . TRA01560 

GA3=GA/ ( CA-l , i TRA01570 

GA12 = G \ 1/GA2/2 . TRAO i530 

PTC = PC« ( 1. t-G A2 *H0**2 )*4GA3 TR I 5<)0 

TTC = TO* I l.+GA2^M0»*2 ) TPA01600 

QCTURN TRAO 1610 

FNO TWA01620 

TRAO 1630 

SLrfRC’JT INF 3U«N TP401640 

CO'-^ON/GFO/A REF, AO. Al. A2 . A 30 . A 3. A 5 . A 6. L 3 , AOAR , A5 AR . ! LOO. I PR TPAO 16§0 

r.CvM(;N/ A IP/G A,GAl,r,A2.GA 12.GA3 .RHOA. TO .PO . UO t HO . GRA V TRAO 1660 

CO^MFN/FL'^L/OHrF ,FT AT. A.N TP AO 16 70 

Cf ‘'“'IN /I CSS/PI Cl .Pf 02, PI3l,PIR2,PIN TP AO 16 30 

CC“Hf‘i/eLR/W A, ,WT , F,R00T,CST AP.PT4.TT4 TP A 016^^0 

CCmhpm/gfi/RF.GF ,GF I ,GF2 .GF12.GF3 TRAO 1700 

CC'VrN/NCZ/H5. FT6.P6 TPA0 1710 

CGVMf.^i/FFS/CF, thrust , ISP .SFC TRAO I 720 

CC““n\/CHG/P ^5 .PS TPA01730 

ccv«f/:/ T M /p TO.T TO TRA01Z40 

CC“Mf- /CCN/P Tl C.vir , Air, ALFA ' TPA0U5O 

CCMMCN/ TH?f’/Ml l,PTll.H2l,PT2l,HSl.PTSl TRAO I 7oO 

CCVHr'N/NS/Ml 2. Pi' 12. H22 .PT2 2, HS2.PTS2 TRAO 17 70 

COMMCM /D IF/H 2. PT2 ,P2 TPAOl7eO 

CnMHn^/CHP/H'»N .V3 I , PT3 ,PT3M,H4 .AS tra01700 

CU“MON/TFA/P I , M'j A. 0 R AG. R PR. LPR.hPR.U.WB.OELU.OELT.T I, TFT A, IRAM, TRAO 1.300 

.IL.ILO.Y.YO. TO0.Xl.Yl,X2.Y2,X3 , Y3 . 3P P3 . R , TORB. I TRA , T ET AO 0 . X OB . Y OB TRAO 1810 
REAL MQ.VIC. Ml l,yi2.H2.H<N.M3l ,M4,H5 .H6 TRA01B20 

PEAL L3.N.ISP. V2l,H22,HSl.HS2.HUA,LPR TRA01830 

A2l = A3 TRAO 1840 

A3SPI/4, ♦(Sy P T (4, /P I *A 30 H-2.*c nnT*Tl )**2 IRAQ 1850 

CCC IN CALCULATpiG WPR IS ASSUMED THAT THE ORIGINAL WPR EXCLUDE WB TPA01360 
;,ppsv.PR-( (A3-A31 )*L3*RHnF ) TFA01870 

IF<A3.L E.AREF ) GC TO 22 TRA'JIRRO- 

thruST = 0, TPA018‘J0 

TOFflsTI TPA01900 

XflP-X3 TPA01910 

YCH=Y3 TRA01920 

FETUFN TPA01930 

22 -A=PMnAML0*4C*.9 TRA01940 

CCC WA=PHOA*L0*AO* ( AC/AO I. TAXING i )».9 TRAOl^SO 

«AA3=«A/A3 TRA01960 

V,AA3*WA A 3*0. 02 54**2/0. 45 35 9 TRAO 1970 

CCC Gl^fS WAA3 IN L8/IN2 TRA019P0 

R0rT = A*wAA3»*N ■'■PAOIOOO 

CCC RCOT=0.06*(WAA3 LB/IN2I**0.6 IN/SEC TRA02000 

PnCT=F00T«0. 0254 TPA02010 

WF=3HnF*P0nT*SCRT(4.*3.l4l59*A 3I*L3 TP AO ^020 

CCC \^F = PHr)F#FOOT*P I*C3*L3.WHEN 03= SQRTf 4/PI* A3) TRA02030 

F*»»F/«A TRA02040 

V.T = wA^wF TRA02050 

CALL TEP“0(T TO .F . TT4,GF.RF ) TPA02060 

IF UL'^n.GT.l )PETURN TRA03070 

TT4 = FTA T* ( TT4- TTC )<-TT0 TPA02080 

GF l=(GF ♦ 1. )/ 2. TRA020a0 

GF2=<GF-1.)/2. TRA02100 

GF12='',F 1/GF2/2 . TRAO 21 10 

GF3 = Gc/< <:f-1 . ) TRAQ2120 

GFUfi = S )P TtGP )*n ./GFU**GF 12 TRAO 2 I 30 

CSTAF = SUPT (GRAV*3 F<«TT4 I/GFUN TP AO 2 1 40 

PT^,»v.T*CSTAR/(GRAV* A5) TPAJ2150 

CCC GIVES PTA IN XG/M2 TRA02160 
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C 



c 

ccc 

c 

c 



c 



c 



c 

ccc 

c 



PfTljrN 

ENO 



SUPPnuTlKE TEQMOiTTOPfFAP.TT^P.GFf^rRFPj'"' ' 

CnHMnN/GFn/ARfF.AO, Al,A2 • ' 30 # A ^ ♦ A 5 . A 6# L3 , AO AP . A5 AP ♦ ( LOO, I PQ 
CALL !NTFR(l ♦TTOr ,FA 3, TTAr^ I 
call JNTFP12 ,TTOF,pa P tGFPl 
CALL rNTFP(3#TT0P,>=AP.RFP) 

<Fiif nn.G^.i )iLnr ^2 

RETURN 

FNO 



SUnP<^MT (NE f NTFP MOATA ,ZP,XP •YP) 

rOHM''.N/r,fn/APFF, A'l, A I.A2 t A30,A 3, A 5 * A 6 , L 3 * A 0 A» , A54?-,(Lnn,IPR 
OI^^FUSION (5) tXV'^r 1201 .YV«-r(^.?0 J.YZ15) r ’ ‘ 

♦ ,TTov=r i5i,F \vpcr‘oi,rFTv ci^» tXMvc: r>o» 

♦ •OATO< 20 ) ,UA TT2(20 ) . T’3( 20 l, HAT’^, ( 2T U0ATT5I 20 ) 

♦ .OATr .2 ( 20 , 3 ATf.^ < ?0 I . 3 AT 0 » ( 201 <201 

♦ , OATO.? ( ->C I *0ATP3( 23 I ♦ 'A TP/* I 20) ,34 Tt?S ( 20) 

♦ ,0ATCC2 < 20) ,mATCC 3(2 3) ,n4TC04( 20) ,()ATr 05(20 » 



OAT 4 rpP Tf RNO SJPRnu^TN'* 

data =4 vEC/0., .01 ..O'*, .0 3,,0<», .0S,.0 57l,.0bo7,.0 73 ,. 07 7,.083, 
♦.091,. I , .111 ,. I26f>* .20 ..21 , .23 , .2 i. . 30/ 



OATA TTOvFC/0. ,311. ,7 50., 83 U/ 

DATA aATC/0.,0.,O.,O.,0. ,0.,0. ,0. ,‘3. ,0.,0., 

♦0. *0.*3. ,0.,0. fO.,0.,0.,0,/ 

r)AT^ ^4 TT2/0., 705., 1 ,) 36<». ,lo^^. ,1 912. , 

♦2082., 2 2 79., 2 3 62., 7 3 80., 23 51., 227b., 2 1 d^. , 20 71 . , 1^20. , 1 326 . 
♦, 1290. , I C7 0. ,9 00. ,-.7^./ 

DATA -3ATT3/0. , 815., I 1 SO. , 1970. .1 751 . ,201 3., 

•2 I 72., 2 3 52., 2'« 29., 2 99 i., 29 2 7., 2 36 J., 22 73., 2 1 60. ,201 5., 1955. 

♦ . I 300. , I 125. ,950. ,515. / 

OATA 0ATT9/0., It 06. , ) 9 22 . . I 719 . , 198Q . , 22 I 5. , 

♦ 2 395. , 29 90., 25 51 . .2 569 .,25 73 ., 2530., 29 59., 2 352. ,2200., I 568. 

♦ , 1920. , 1290. , I 115,.6^2 ./ 

OATA 'A4rT5/3.,ttP3.,19 99., I 7 86 . , 2092 . , 22 ^0. , 

♦29C(3.,2 5 32.,2 5 89. ,2 602 . , 26 06 . , 257 , 2508 . , 29 09 . , 225 5 . , 1619. 

♦, 1520., 1 390. , 1 165., 71 1 ./ 




36/ 

36/ 

36/ 

36/ 



DATA OATR P/0 ., 29. 03 ,39. I 0, 29. 1 6, 29. 2 2. ■’9. '*8, 

♦ 29.32.2 9.29, 2 9.18,2 9.39, 2-^ .69, 2 8 . 2 1 , 2 7 .6 h , 2 7. 3 0 , 2 6 . 29, 23. 0 3 , 
•22.58, 21 .65, 20. 75,1 8,5/ 

hat A 14 T F*’ /O ., 29.01,29 .10,29. I 6,29. 2 3,29, 28, 

♦2 9.30, 2 9.2 9, 29 ,1 1,2 8.9 7,2 1.6 0, 28.2 J, 2 7. 60, 27. 00, 2 6. 20, 2 3. H 7, 
♦22 .5r, ? I . 7,20.8, < 7. 7/ 

OATA hat 09/0 ., 29 .03,29.10, 29.1 6.29, 2 2, 29.20, 

♦29.20,29.10, 28 .93 ,2 rt. 79,2 S. 5 6, 23. JO, 27.60,26.16,26.20,22.97, 
♦22.55, 2 1 .6,20. 7, 18.9/ 

DATA 0ATP5/0 ., 29 .02 , 29.10,29.1 6,29.2 I, 29.20, 

♦29.20, 29.10, 28.90,26. 73,28.51,28.00,27.58,26.15,26.20,22.97, 

♦ 22.55,2 I .6,20. 7,TH./,/ 



r)4T\ cn^ QOAGN ci'IrtPUiriN^ 

3 AT A TTvrr/0, ,6.5,9.5,12.0, 15 .0/ 

0AT4 <«VFC/0. , 1.m 5,2.J0,2.10,'' .23. *. 39,2 .J5,2.S0,2.95,2.50,2.55 
♦,2.60, ?.65,2. 70, ?./5,?.rtO,2.85 ,2.90, 3. JO/ 



OATA 04 rCn2/O.0,O.l I ,0. j766, O.J 72 6 ,0.0^9) , 0 . 06 59 , 0 . 0 Y 9 ♦ , 0 .37. ’0 
*,0.06) 7, C.J609,J.C59 2, 0. 05 31 ,0 .05 7 , 0 . 0559 , P. 0*> 9Q , 0. 05 ?9 , J. i) 5 A 0 



TRA02170 
TR 40 2 180 
TP 40 2190 
TP 402200 
TP40-*?10 
TF 43 P2?3 

43 2230 
TP432290 
TP 40 ^250 
^243 ^^*60 
'P . 4 0 ^2 70 
TO 4 3 2280 

T o • 3 

*’140 2330 
70*02310 
TPA32120 
TO4O2310 
TP4 0 "• l'»0 
■"PA } > 35 ) 
TO 40 2 36 3 
TOAO ■*370 
TP 40 P180 
TPA02 390 
T P40 7605 
TO 4 3291 0 
7P 40 "*9^0 
*’3 40 2910 
TP 40 74,90 
'■0.40 ^653 
TP 40 297,0 
TP 40 "*9 70 
’P40 7410 
7P 40 2490 
TR40 ■*SP0 
TP 40 95) 0 
7R4 02523 
rR402570 
TP 4 0254 0 
yo •.0"*S50 
"P 40 25aO 
^P 40 7570 
T040 7*^80 
TC 402590 
TP 40 :*600 
R4026)0 
TP 4 0 262 0 
TPi.0 ">570 
■rp 4 J2640 
TP402650 
TR40 -*660 
TRA07670 
TP432680 
TR402690 
TR 4 02 7^70 
TO ^ 0 “*71 3 

4i) 2773 
TO 402 710 
■^P40 ^ 7<i3 
"0 4J 2 750 
TR 40 2 T60 
R 4 0 2 7 70 
TP402780 
TR 40 **790 
"040 ‘*800 
TP. 4 0 28 10 
TRAi •»870 

4 ) 2 r»o 
T r vj ^ p s, 0 
1040 ’550 
TRAP ‘9^0 
*’0 4)'* I/O 
rp4028H0 
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C 



c 



♦ ,0.05?l ,0.0504/ 

DATA DA TCn3/ 0.0,0. 1 633,0.1 I 06,0. 1 04Q ,0 .044 , 0 .0 R 53 , 0 . 0^3? , 0 . 09' ? 
♦,0.0393, C. 08 75, 0.08 53, 0.0341,0.0826, 0.0811,0.0796,0.0782,0.0769 
♦ ,0.07'44, C.0732/ 

DATA r>\T(^n4/ 0.0,0.^210,0. I 5 04, 8 ,0,1 36 1 ,0. n 00, 0.1 772,0.1 r ^6 

*,0.1720 ,0.1196,0.11 7,0.1151,0. 1L3, 3. 1109, 0. 10^*0, 0. I 0 71 ,0. 1054 
*,C.1036 ,C. 1004/ 

DATA OA TC05/ 0.0, 0.1 009,0. ’I 36, 0.203 2 ,0.1 93 8, 0. 1 355,0.1316,0. 17 79 
*,0.1 744,0.17 10,0.16 78,0. 164 7, v) .1618, 0. 1590, 0. 1 56 J , 3. 1 5 3 7 , 0 . 1 5 1 2 
*,0. 1433, C. 1443/ 



21 

22 

24 



25 

29 

31 



51 

33 



53 

35 



55 

37 



57 



r.O TP (21,21 ,21,24) , HATA 

TZ^O 

ini=20 

DC 27 I ,5 

ZVFCn 2 )»TTOVFC( TZ) 

r,C TP 29 

IZ»0 

I0A=19 

00 25 (Z*l,5 
zvrrnz) = TFTvEC(iZ) 
on ( 31,33,35, 17), lOATA 
( F AsO 

on 51 1 5 A=» , fOA 
XV-f( ( F \ )=FAVCC( IFA ) 

YVFf( 1 , I FA)snATO( (FA) 
WFCi 2, I FA )=P.\TT7( I FA) 
YvFrn , I fa) = oatti ( ifa) 
YVFC(4, I F:. )sPA TT4U FA) 
YVFC(5 , IFA)»OATTSnF.A) 

GO T'^ 19 
IFA=0 

OC -=3 I FA=1, fPA 
XVFf( I r A )=FA VFC( TFA ) 

YV^CI 1 , T FA)«PAT0( !FA ) 
YVEa2, IFA)sOATC2nF4) 
YVFC( 1 , t FA )=0A Tf. U (FA) 
YVFP(4 , ! FA )*0A T(]4 ( IF A ) 
YVFCI 5, I FA)aOATG5( ( FA) 

GO TO 39 



( FA=0 

CO T FA=1, IDA 

xv‘^ri r F A )=»-A v‘=c( ir* > 

YVf^n I . I FA)sOATO| !FA) 
YVcCI 7 , I FA )sOA T92M FA) 
YVFC( 1 , I FA )=nA TR M FA ) 
YVFf(4, ( FA)aOATP4(( F \) 
YVfC( 5, I FA )*nATR5( ( FA) 
GO "^P 39 
IFA = 0 

on 57 I FA»l, lOA 
XVFC( (FA ):=)(MVPCI (FA ) 
YVFCI I , ( FA)=PAT0( IFA ) 
YVFC(2 , I FA)snATC"'?l ( fa ) 
YVFCn , ( FA)*DA rr03( I FA ) 
YVFCI 4, I Fa U0ATC^4( T FA j 
YVFCI5, IFA)»OATC05( IFA ) 



39 IF ( (7P M .OF-06 ).LT.ZVFC( 2) )G0 TO 67 
IF ( (7P-1 ,or-C6 ) .GT.ZVrci 5) ) GH TO 67 
I F( (XP f I .0E-06).LT. xvrc( 2) )G0 tq 67 
IF KXP-l .0F-06).G^.XVrC( PA) )G0 TO 67 
IF( IOaT A.GT, 4)Cn TP 6 7 
IFA*1 

OC 61 I FA«2, lOA 

I FI (XP- 1 .OF-06)-XVFCI ( FA) ) 63, 6 3, 65 
63 IT=1 ' 

on 71 IT*7.5 

1F((/P - 1 .0‘ - 06)-ZV'*C( ; T) ) M, 7^ ,75 
73 ITPl 

on 81 I TIM, 2 
IT0=IT-2MTl 



TRA02890 
TP AO^OOO 
TPA07910 
TP A 0 2970 

TPAO '»9'>0 
T Q A 0 7 94 0 
TPA07950 
trap 2960 
TP AO ’970 
TP », 0 ->930 
T3 AO ?990 
TP AO iOOO 
TRAO 1 0 

TPAO )0?0 
TPAOlOlO 
TPAO ■'040 
TPAOlO'iO 
TRAO’060 
Tp AO 'O 70 
TP A 1 1080 
TOA.) ’ )90 

"= AO mo 

TP AO 11 10 
TRAD 1170 
TP AO ^ 130 
TR A 1 1140 
TRA J 1150 
TR AO 3153 
TRAOll 70 
TC AD 1180 
TRAO U90 
TO AO 1700 
TRAiJ J21 0 
TO A 0 12 00 
TF AD 1210 
TR A 3 1740 
TR.AO 1250 
TP AO 1260 
TP AO 1270 
TP*0 1780 
TR ^ 0 1790 
TRAO ^100 
TO AC 11(0 
’■RAO 1170 
T R • n 1810 
TO AO 1340 
TP A J 1150 
TPAO 1150 
TRAO 1370 
TR AO llHO 
TO AO 1390 
TRAO 1400 
TP* 914 10 
’■P AO 1420 
Tp AO 1410 
TRAD ’440 
TRA 01450 
TRAO 1460 
TC Al 14 70 
TR AO 1430 
Tg •n'»4«0 
TRAO 1500 
T P A 9 '1 5 ) 0 
TRAO 1570 
TRAO 3510 
TP A I’S^O 
TR A.l »S50 
TO A O ,560 

^RA*^ ^5 73 

’ f. A 3 15 3 0 
'PA.) 1590 
TRAO 1500 
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Y2 ( ITO) a (XP-XVECI IFA-l n /UVEC (IFAI-KVECI IFA-l M* 

♦( YVFCI ITC, <FA)-YVEC( I TO , (I F A- I )) ) 4*Y V EC (I T 0 » H F A- I ) ) 
81 COSTINUF 

YP*UP-2V£C( I T-l ) )/( ZVEC ( I T )-2VEC( IT-U) « 

• ( Y2( I T )-YZ ( I T- 1) ) ♦Yil IT-l ) 

D ^ T'JON 
75 CCNTIN'JF 
71 CCMIN'JF 
65 CC^'^I^UF 
bl CCNTlNtJF 



6 7 ILCP=10 

I F( 1?^. ,lT.D PETU«N 
PRINT 69.A5AR , AOAR, ZPf XP, I DATA 
69 FnPyATUx,2F7,3,5x,2Fl2,^. I3,5X,»1ISS!NG OATA TO INTER*) 
octupn 
FNO 



StPFt UT INF NOZZ 

CnMMrN,/r,FQ/AREF, AO, A I , A 2 , A 30 .A 3, A 5, A 6, L 3 , AOA R, A5 AR , I L 00 , 1 PR 
CCVmon/A I»/OA,GAl ,GA2,GA12,GA3,RH0A, TO tP 0 1 UO ♦ 0 » GRA V 
N / FLtL/PHCF , FT at , A, N 
CC^-r‘4/LCSS/PI Dl ,?I02,PIBl ,PI32.PIN 
CP'^'^CN/PLF/W A, WF ,uT , F.POQT.CST AR, PT4,TTA 
CO^MCN/uf I /RF ,GF ,CF I ,GF2 ,GF12, GF3 
CCP^^N/NCZ/'^c, PT6 ,P6 
CC,'*MCN/FES/CF . ThFgST . (SP.SFC 
ccr«r*:/r>c/p T5 ,R5 
CCVM-^/ f M /? TO ,TTO 
CC*'*-rN/rC*/P T 1 C,‘"IC,AIC. ALFA 
COMMON / TFFu/Mll,PTU,M2l,PT2l*MSl,PTSl 
CCHMf'N/ 2 , PT 12,M22 ,PT22tMS2,PTS2 
CGmmp*: / rj IF/M2, PT2,P2 
CCMMCN/CFE/M ,M3 t , P T 3 , P T 3M , M4 , A S 
R F al « 0 , ir , M n , y 12 , .M2 , M 3N , .M 3 ( , M4 , M5 » M6 
real L 3.N,ISP,M21,m22.MSI,MS2 
M6- 1.0 

CALL C AICM(M6, A5, A6 ,oFl,GF2,GF 12) 

IF( ILOn.LE.l JGC TO ll 

ILCC*3 

RETURN 

ll PTc=PT4*FIN 

P6-PT6* ( 1.»GF2*M6**2 )**(-GF3 ) n 

RFTiJOM 

END 



ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 



4 

5 

6 

7 

999 



SU^POtlT INF C ALCM(X,AN, AM,GI,G2 ,012) 

COV‘>ON/OFO/AKFF, AO, Al, A2 . A30,A3,A5, A6,L3,A0AR, A5AR, ILOOt IPR 
this SFC^ION CCmPuTFS MACH NUMBER *=R CM AREA RATIOS* 

THF PASIC FOQ.muLA IS: 

AN/AM= { Gl/(l •►G2<‘X**2 ) ) *«G1 2*X 
OaTA NCCCcO: AN = AR fa \T <Nf WN MACH number 
AMsARFA at new mach number 
G l ,G2,G12 = GAM^A PATIOS FOR: 

-AIR FUITMI GA I ,GA 2,GA12 ) 

-HOT FL0«( GF I ,GF2,GF12 ) 

X= KNOWN M\CH number AT AN(=XAN) 

X IS THE CCWPUTEO MACH NUMBERITHE SAME X IS USED FOR IN/OUT) 

( mach* I 

XANsX 

IF(XAN-l •0)‘t,5,5 
ALFA=0,99 
GO TO 6 
ALFA*1 .01 
COM INIJF 



EPS*l.E-C^ 

r,FUN=( G I /( 1. ^G2*x**z) I 
F < = X«GFUN**G12-AN/AM 
I F ( IFF . FC.3) PR INT 999,X,FX 
FOPMAT ( 6 >,2F 10.4) 

F <D'"T= ( GFUN* *G 12 ) ♦< U-2. *G2*Gl 2/ALFA •( X**2 ) / ( I . ^G2*X**2) ) 

XNF« = X-F »/FxnoT 

0<Fr=ABS(X-XNrW) 



TRA0 3610 
TRAO 3620 
^RA03630 
TRAO 3640 
^RA036S0 
TRAO 3660 
TP AO 36 70 
TRAO 36H0 
TRAO 3690 
TRAO 1700 
TRAO 3710 
TRAO i720 
TPAO 3730 
TRAO 3740 
TRA03 ;S0 
TRAO 3760 
TR Ax) 3 7 70 
TRA0 37B0 
TRAO 3790 
TP^ J3800 
TR ACHHIO 
TRAO )B20 
TR AO 38 3 0 
TR \Q ^640 
TRAO )B50 
TR \0 3d60 
TRAD 3870 
TRAO 3RR0 
TRAO 3890 
TPAO 3R00 
TR AO 391 0 
TRAO ^920 
TRA0 3930 
TRAO 3940 
TRAO 3960 
TRAO 3960 
TPAO 39 70 
TRAO )980 
TRA03990 
'^PA04000 
TPA 0«*010 
TR A04020 
TRA04030 
TRAO 4040 
TRAO 4060 
TRAO 4060 
TRA04070 
TR A040B0 
TP A04090 
TR A04100 
TRAO 4110 
TRA04120 
TRAO 4 I 30 
TR A04140 
TR AO^l 50 
TRA04160 
TRA04170 
TRA04180 
TRA04190 
TR A04200 
TRA04210 
TPAO 4220 
TRA04230 
TRA04240 
TRAOM250 
TR \04260 
TRA04270 
TR A04280 
^R AO 4290 
TRAO -*300 
TP A04310 
TRA04320 
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G 



IF(0IFF.LE,£PS)RETUR1M 

IF( I^ACH.GT. lO)GO TO 9 
IMACH- I MACH*' 1 
X* XNFW 

P=lXAN-l.C)l<», 15,15 
IFtX.GF. 1.01 X=0.7 
. I F(X,L=.C.iX=0.3 
GO rn 7 

15 IF(X.LF.l.J) X=1.0 
1 P(X.GP.lO.) X-3. 

GC Tn 7 
9 lLCn=lO 

IF( !?o .LT.l) OFTUPN 
PRIMT 10,A5AP,A0A« 

10 <=OR'*AT ( I x,2F 7. 3,3X, • SUTTA0LE SOLUTION WAS NOT FOUND AFTER 6 
IF( [P° .GE.2)CALL PR IN 
kFTIJRN 
END 

SLDR<"U'^INF RESUL 

CCyN,lN/G = n/ARFF,AO, Al, A2,A30fA3,A5,A6,L3,A0AR, A5AR , I LOO, I PR 
/\ IS /GA,GAi,GA2, GA12 ,GA3 ,RHnA, TO,PO ,U0, ^0, GR AV 
COm.-^C'M/FL EL /RhCF .FIAT, a, N 
CC‘'MrN/LCSS/ PI 01 .P1D2,P1QI.RTD2,PIN 
CC«NrM /P LR/W A, f,F ,WT , F ,ROOT ,CST AR, PTA ,TTA 
CC‘-H( n/GFI/R F.GF.CF 1,GF2,GF12,GF3 
Crv^t'N /NC2/«6 , PT7, ,P6 
CCVMr.N/ RFS/CF , TH'^UST ,l SO ,SFC 
CC'^yp»;/CFn/P T3 ,05 
IM /P TO,TTO 

ccyyoN/ccN/p t ic.^ic, mc, al fa 

CC'^yr *J/TFR0/^1 I, PTl I , ^21,PT2 I, HSl ,PT SI 
COyMC'N/N«/Mi 2, PT 12.M22 ,PT2 2,MS2,PTS2 
COyHON/DIF/N 2, Pr2,P2 
CC^M( N/CFF/N3N,*^3I, PT1,?T3M,M^ ,AS 
RFal mo l,yi2 ,M2 ,m3N,M3I ,M<.,‘^5 ,M6 

peal L3,A,lS0,y21,^^22,MSl,HS2 
pp AT^PT6/P0/ (P T6/P6 J 

CF = 2.*A6/Ar-c F/ jga*wO**2) •( PRAT ♦(! )-2»A0/AR6F 

IFfCF.GT .0.) GO TO 22 

lLCn=9 

RETtJFN 

22 C0=.5*GA*P0*^'0«*2 

rHRUST=rF*CO*ARFF*GRAV 
I SP=thouST/wf 
SFC==3600./ISP 

01 yE*^'SinN5: THRUST IN NEWTONIN), ISP IN N/(KG/SEC), SFC IN 

RETURN 

ESC 



SUPROfjT INE CHOKE 

CGVMrN/GFO/AREF, AO, Al, A2 , A30 , A 3, A5, A 6, L3, AOAR, A5AR, ILOO, IPR 
CCyyPN/A IP/G A,GAl,GA2,GA 12,GA3 ,RHOA, T0,P0,U0,M0,GRAV 
COyNON/ Ft EL/RHCF .^T AT, A.N 
COMMI N /Less/ 0 101 , PI 02, PI Bl ,P IH 2,PIN 
CCMMLN/rlF/W A, -F ,WT ,F, ROOT, C STAR, PTA,TTA 
CC'^NCN/GFI /P F.GF ,GF I , GF2 , GF 1 2 , GF3 
N/NCZ/M 6, PTe,P6 
CCFNLN/FfS/CF , thrust , ISP, SFC 
CCyyON/CH^/P T5 ,P5 
CCFMON/ IM/P TQ,TTO 
CCV‘'^N/CCN/PT1C,MIC,A1C, ALFA 
COf^MON/THP :/Hl l,PTl 1 ,M2l ,PT2 I, MSI ,PTS1 
CC.''Mo*-/n£/'- 1 2, PTl2iM22,PT22,«S2,PTS2 
COMM^u/D IF/M2, PT2,P2 
CC«MPN/CHF/'- 3N,M3 I,PT3,OT3M,M4,AS 
REAL MO IC, ''I l,y 12 ,'-2 ,N3N,M3I ,N4.M5,M6 
PEAL L 3 ,N , ISP, «2 I ,M22,MSl ,MS2 
PT5-PT^ * SORT (PIN) 



P5=PT5« ( l./GTl )^*r,F3 
IF(P5.G6.P0I pFTUR.N 



TPA04330 
TRA0<*340 
TRA0<V350 
toa043^0 
TR A 0^,3 70 
TRA0<»380 
TPAO ^390 
TR AO^AOO 
TRA04410 
TRA04420 
TRA04't30 
TP AO <.'^40 
TR AO^^SO 
TR AO 4460 
TP AO 44 70 
irgR, ITP A04480 
TO AO 44^0 



TO AO 4500 
TR AO 4510 
TP A04S20 
TR A 0 45 30 
r RAO 4 540 
TRA04550 
TR A()45^0 
TOA04570 
TR AO450O 
TRA04590 
TR A04600 
TP A 04610 
TP A04620 
TP A04630 
TRA04640 
TP A0^650 
TPAO 4660 
TPA04670 
TRAO460O 
TP A04690 
TPAO 4 700 
TR A04710 
TP AO 4 720 
"RA04730 
TP AO 4 740 
TPAO 4 750 
TRA04760 
TPAO 4 7 70 



KG/N. 



TPAO 4 7fi0 
TRA04790 
TPA0 4H00 
TRA04B10 
TR A04H20 
TPA04H30 
TPAO 4040 
TPAO405O 
TP 30 4860 
TR A04870 
Tp A04880 
TRAO409O 
TRA04900 
TRAO 4910 
TRAO 4920 
TRA04930 
TPA04940 
TO AO‘,950 
TRAO 4960 
TP A04970 



TO 304900 
TRAO 4'y90 
TRAO 5000 
TRAO 5010 
TPAO 5020 
TRAO6030 
TRA05040 
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1LC0=^ 

1F( return 

PRINT 53 .A5AP , AOAR* P5 

53 FCPHAT ( IX,2F 7. 3*3X, »N0ZZLE IS SOT CM OKED , P5- • , F 1 0, 2 ) 
IF( IPR ,GE.2}CALL PR is 
return 

PSO 



SUbRPMT INF I NL ET 

CC''«'"N/r,Fn/ARPF, AOf Al,A2fA30*A3,A5, A6fL3»A0AR*A5AR * ILOO.IPR 
A I R/G A,GAI,GA?,GA12,GA3 ,PH0A, TQ ,P0 t UO * MQ t GR A V 
CO^»^r J/ FL5L/RH0F , FT AT, A, N 
CCMMCu/t CSS/ PI0l,PI02.PIBl,PIS2»PIN 
CCM-PN /RLR/W A, WF ,«T , r ,ROOT ,CST AR, PTA ,TTA 
Cn'— o*:/Gf i/Rp ,GF ♦GPl .GF2 ,GF 12. GF3 
COFMON /NCZ/'16. PT6 .P6 
CpwM( /ags/CF, TH«?lST . I SP .SFC 
CC^'‘-r N/r fp/0T5 ,p5 

cc««f'N/ 1 M n Tc ,ttg 

CG'-'*cs/CCN/p t IC.^'IC ,AIC. Alfa 

TFPn/Mi I ,PTI I ,.M21 ,PT2 1. MSI ,PTSl 
Cn''Mr»,/NS/-l 2. P^ 12. M22 .PT2 2.MS2.PTS2 
CCM«f,S /O IF/M 2 . PT ? .P2 
CCMMn^J / r FF/M j \ I , PT3 .P T3 M.m<, , AS 
5 FAL -0 ,m IC , ''I 1. f 12 .M2 .M3N. M3I ,M6 

REAL L3.N,ISP,M2l.M22.MSl.MS2 

call CPNE 

CALL I S 

IF( IlPQ.LE.l )GC TO ll 

I LC0 = 6 I 

RFTUPS 

ll call THSCATI l.Mii ,Ai j 
I F ( ILOO.LE,1 )GC TO 13 
iLrn=o2 
rf tups 

13 CALL SSR(l,Mn,Ml2.PTll.PTl2) 
call niFFUSI 1. Al,yi2) 

IFI ILOC.GT.l ) 1L0C»63 

RETURN 

ESC 



ccc 



SLPPnilT INP CONF 

CCMMPN/GEO/AREF. AO.AI .A2 , A 30 . A 3 . A5. A 6. L3 . AOAR, A5 AR . I LOO. I Pft 

CCVMCN /A IP/GA. GAl,GA2.GA 12 .GA3 .RHOA. TO. P0» UO » MQ . GR A V 

CCM«nN/FLEL/ PHCF .FT AT, A. N 

CCMVrs/LCSS/ PI 01 .PI 02.PI81 .PIS2.PIN 

CC'^’Mf'N/RLP/M A. wF ,«T,F,R00T,CSTAR.PTA,TTA 

COMMON /GF I /R F.GF .OF I .GF2 .GF I 2. GF 3 

CCRMf/j /NCZ/M6. PTC .P6 

CCM'-ns/RES/CF. thrust .ISP.3FC 

CC««rN/CHr»/{> T5 .P5 

CCMvr>i/IM/P tq.TTC 

C^^«M^^i/ccN/P TIC. MC.A1C.ALFA 

CCMMCS/ THO^Z-l l.PTll .M21 .PT2 1. MSl.PTSl 

rCF'-rr;/NS/'’l 2. PT 12. M22 .PT22 .MS2.PTS2 

CnMM,:N/OIF/M2.PT2. P2 

CCMMn\/CHE/M3N .M3I. PT3 .?T3M.M<* ,AS 

RFAL mo , y IC . Ml I , V 12 .M2 ,M3 N,m 31 .MA.M 5 ,m6 

REAL L3.N.ISP.-2l.F'2?.MSl.MS2 

PEAL MN,Rl 

ALFA IS THt^ half ANGLE OF ThE CONE 

CP = ( .0B3 + .096/M0»*2 I *( ALFA / 10. )*«1.69 

P IP0=1 .♦'CP^A/2*MC**2 

MSsS';*- T( !.♦( PlPO-l. ) *( GA 4 .I ) / (2.*GAn 

PTlPTu-P lPO**<-l./(GA-l. n*( (GA4-1 i»MN**2/ 

/ ( I GA-l } 2 «-2 . ) } *«t GA/ (GA- I } ) 

BETi = uS I N ( '’N /MC) 

TcTAsATAN(2.«C0TAN(5FTA)*(MN*»2-l.)/ 

/ ( «0**2 * (GAtC0S(2’»«FTAI )^2.M 

1=SC- T ( ( I ( CA-l )/2 *mn**2 ) / (C 2-(GA-1)/2 ) / 

/(SlNtrt^^TA-TFTAl ) 

AUOM rl-TAslTETAJ/fAN(3ETAJ )*COS(TETA) 



TRA05050 
TRA05060 
TPA05070 
TPAO 5050 
TR AO 5090 
A05100 
TPA051 10 
TR A05120 
TRA0513O 
^PAOSl-^0 
TPA05150 
TPA05I60 
'^R A051 70 
TRA05130 
TP AO 5 I 90 
TR AO 5200 
TRA 152TO 
TPA0 3220 
TRAU 52 30 
TRAO 52^0 
rPA05250 
TP A05260 
T9 A05270 
TR A0S250 
TPA05290 
TRA05300 
TP A05310 
TPAOS320 
TRA05330 
AO 53A0 
'TRA05350 
TR AO -)3^0 
TRA05370 
TPA053fl0 
TRAO 5390 
TPA05400 
TR A05A I 0 
TRAO 5A20 
TP AO 5A30 
TRAO 5AA0 
TPA05A5U 
TO AO b-VbO 
TP AO 54 70 
TRAO 5450 
TP AU5490 
TRA05500 
TP A05510 
TPAO 5520 
TRA05530 
TRAO 5540 
TRA05550 
TRAO 5560 
TR A05570 
TP A05530 
TPA05590 
TRAO 5600 
TRA0561 0 
TRAO 5620 
^P A05630 
TPAO 5640 
TR A05650 
TRAO 5660 
TRAU5670 
TR AO 56B0 
TRA05690 
TRAO 5700 
TRA0 57 10 
TRA05720 
TRA05730 
TR A05 740 
TOA05750 
TRA0S760 
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PTlC=PTO^PTlPTO 

AlC=AO*AlAO 

return 

ENO 



SUBRPUTINP MIN 

CCMMO^/oeO/AP EF, AO. A1,A2 .A 30. A 3, A 5, A 6, L3.AOAR, A5AR . I LOO. I PR 

CC'^MON / A IR/r, A.GAUGA2.GA12.GA3 .RHOA, T 0 ,PO ♦ UO . MO . GR A V 

CCf''-rN/FtEL/ PHCF ,ETAT. A.N 

CQmmpn/LCSS/ PI 01 .Pt02.PIftl.PlB2,PIN 

COMMON /eUR/W A. UF ,WT , F ,OOOT ,CST AR . PT^ , TT4 

CCPMON/GF I/R F . GF .GFl .GF2.GF12, GF3 

common /ncZ/m 6. PT^.P6 

CC •omon/’^F $/:*=, thrust,! St>.SFC 

CnMVPN/Chr/P T5 ,Pb 

CCVMf '!/ I M /P TO .TTO 

CC^--r N/rCN/PTl C,'-* IC ,AIC, ALFA 

CC^M^-N/THfO/'-l l.PTll ,M2l ,PT21. MSI ,PTSl 

CCM‘^rN/N5 / n 2, PT 1 ?,.m22,PT22.mS2,PTS2 



CCPm< n/:) IF/’^ ?, PT? ,P2 

CQMMnN /CHC/M jN ,M3! , PT3 ,PT3M, M^ ,AS 

pCAl MO.yir.Hl 1 ,m12,M2 ,.M3N,M3I ,.m<»,M5 ,M6 

oEAL L 3 .N. rSP,M2l .M22.HSI.MS2 

CALL ThPCaK 2 ,M2l .A2 ) 

CALL NSP(2,M2l,M22,PT2l,PT22J 
M2 = .m22 

PT2^PT2? 



P 2 = PT 2 * ( I . a 2 *M 2**2 )*♦( -GA 3 ) 

CALL =tPAN 
I FdLOn.CT.l )RFTU«N 
IF IPT 3 M.LE .P T 3 )PPTURN 
I Lro= I 0 

IF ( IPR.l T.l) PFTURN 
POINT 33 .A 5 AP.A 0 AP 

.33 FQPMAT ( i>,2F7. 3, 3X. *PT3 IS LESS THAN MIN VALUE ALLOWED*) 
I F ( IPP . GE. 2 )CALL PRIN 
RETURN 
ENO 



SUBROUTINE THROATd ..MIl.Al) 

CCM'^O'j/GEO/AREF.AO, Al.A2.A30,A3.A5.A6,L3.A0AR,A5AR, iLCOtIPR 

CCMMnN/AlR/GA,GAl,UA2,GAl2.GA3,RHOA. TO, PO ♦ UO ♦ MO , GR AV 

Cnv>^j'N / FLtL/OHC*=,FTAT, A.N 

Crv“CN/LCSS/ PI 01 .PI02,Piai .PI82.P IN 

CC^'MCN /BLR/W A, *.F ,V,T, F ,Pf)OT .CST AR, PT<^,TTA 

COMMOV/Gf I/R F ,GF ,CF I .GF2.GF12, GF3 

COMMf N /NCZ/M6, PT6 ,P6 

cg*^hln/res/cf, thrust, ISP, sfc 

CCvHON/CHn/PTS ,P5 

CCH*^nN/ I M /o TO ,TT0 

CONMUN/CCN/P r ic ,M IC.AIC, ALFA 

CO'^Mn.j/ THPn/ l,PTll,M2l.PT2l.MSl,PTSl 

rCMMGN/N'S/Ml 2. PT 12.M22 ,PT2 2, MS 2, PTS2 

CGV'-oN /n if/m 2. FT2 ,P2 

CC'^MCM / CHE/M 3N , Ml { , PT3 ,P T3M , M/, ,AS 

REAL ''0, MIC, Mil, H 12,. M2 ,M3N , M3I ,MA , M5 ,M6 

real L3,N,ISP.“2l,M22,MSl,MS2 

REAL Ml I 

AN=AIC«M IC 

Gl = l.i-GA2*MIC**2 

MI l = .MlC 

CALL CALCMIM II .AN.A I ,GX, GA2.GA 12) 

IFUinO.GT.l )RETUPN 

IFIMl I .uE.l. )G0 10(17,18,19),! 

ILrnn.o 

IFI IPr'.LT.DPFTURN 
POINT I t> ,A5AP , AOAO 

15 FCRMATI IX,2F7.3,3X, ‘UNSTART CCNOITIONS*) 

IF( I0R.GE.2)CALL PR 1 N 
RETUrN 

17 PTU = PTlCMPIOl 



TRA05770 
TRA05780 
TRA05790 
TRAO 5800 
TRA05810 
TRA0S820 
TOA05S30 
TRAO 

TRAO 5850 
TRA05860 
TP A05R70 
TOAo SttBO 
TOA058M0 
TRA0S900 
TO AO 5910 
TR A05920 
TOA05930 
TRAO 5940 
TR AO ‘j950 
ToaOom^O 



TP AO 5W70 
TRA059H0 
TRA059QQ 
TRA06000 
TR A06010 
TRA06020 
TO A0O030 
TRAOuOAO 
TO AO 0050 
TRAO 6050 
TRAO 60 70 
TR A06080 
TRA06090 
TRAO 6 I 00 
TRA061 10 
TOAO 0 I 2 O 
TOA06130 
TP AO 61 AO 
TOA06150 
TRA06160 
TRA06I 70 
TO AO 6 130 
TRA06 190 
TOA06200 
TRA06210 
TR A06220 
TRA06230 
TRAO62A0 
TRA06250 
TO A06260 
TRA062 70 
TRA06280 
TRAO 5290 
TR A06300 
TPA06310 
TR*0o320 
TRA06330 
TR A0o340 
TR A06350 
TRAQ6360 
TRA06370 
TRA06380 
TRA06390 
TO A06400 
TR A06410 



TO AO 6420 
TO 106430 
TO 10 0440 
TRA ) 6450 



TO A06460 
TRAOohTO 
TRA0O480 
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RETURN 

18 »T2l = PTlC*PIOl*P102 
RPTURN 

19 FAC-I A I- Al)/ (42-41) 
PI02S=l .-FAC*( 1.-OI02J 
PTS i-PTl C*Pl Ol*P IC2S 
pptu®n 

ENO 



ccc 

ccc 



SU3RCMT IN6 N SR (I . MI 1 I 2 , PT I I , PT I 2) 

CC^Mr%/GFO/AREF, AO, A I, A2 .A30.AT,AS, A6,L3 t AOAR, A5AR, I LOO, 1 PR 
CCVMr.'^ /i iR/G a,CA I ,GA2,GA12,GA3 ,RHCA, T 0 ,P 0 , UO , 0 , GR AV 
CCV'^'.'N/FLFL/RhCF ,f TAT, A,N 
CO'-^CN/LCSS/^I ni ,?I02,PI8l,Pt92,PIN 
CC^Mr-i/RLR/WA, V,F,wT,F,ROQT ,C ST AR, PTA , TT^ 

CCmm(jN/GFI/RF ,GF ,CFl,GF2,GF12,nF3 
CC^VCN /NC2/'-6, PT6 ,P6 

CQ'-'^rN/Fes/CF, tmpijst, I sp,sfc 

CC'^'^'IN/CFC/P^S ,P5 
CCN'^'.N/ I K I /P TOfTTO 
cc^'^'^N/cCN/p tic,mic,aic,alfa 

CQ^HCN/ TKtr*/ Ml i,PTii,.-«2i,pT2i,MSl,PTSl 
COf'Mr-ZNS/vi 2, PT 12.M22 ,P^22,MS 2,PTS2 
CO^MGN/0 IF/M 2 , PT:,p? 

CnyMrN/CFE/M3N,‘^3 I , PT3 ,PT1,M, M4 ,AS 
RFAL ‘^0 ,»-lC,.Mi i,M 12, -^2 .M3 .n,m3I ,M6 

REAL L J ,N, ISP, ^21 .H22,MSl,'-S2 
REfiL I I ,MI2 

THIS SECTinN CCPUTFS MORHAl SHOCK RATIO LOSSES 
RFF. TO PT12 ,m12-SHAPIRQ VOL I , P . I 18- 1 1 9. 

RTI2UpTII*IGA l*Mii«*2/( l.>GA2*Hl l**2) I* ♦GAD 
PTI22= ( GA/GA I •y I l**2-0.5/GAl2) *♦( 0.5/GA2) 

PTI2 = P T I 21/PTl 22 

MI2^SORTUMl l**2^l./GA2)/( 2. ♦G A3»M 1 1 ♦♦ 2- I • ) ) 

RETUFM 

ENO 



CCC 

ccc 



31 

33 

35 



SLbi^r*JTIKF DIFFUS(I,AI,MI2) 

C0.“MOn/GF 0/4REF, AO, Al, A2 , A30 , A 3 ♦ AS , A 6 , L3 , AO AR , AS AR , ILQO, IPR 

CC'^MOn/A IR/GA,GAl ,GA2,GA 12,GA3 ,RHOA, TO , PO , UO tMO, GRA V 

CnyMCN/FLEL/RHCF ,FTAT,A,N 

CC'^MOM/LCSS/ PI Cl ,PI02,PI9l ,Pia 2,PIN 

fGVMPN/flLR/w A, ,wT,F, ROOT, C star. PT^,TT4 

CQMMON/GFI/RF,GF ,GF I ,GF2 ,GF I 2, GF3 

CoyyoN /Ncz/M 6, pt 6 ,p^ 

COMMON / R ES/C F, THC. US T, ISP ,SFC 

CC'^yCN/CFO/oTS ,PS 

CCy«fN/IM/PTQ,TTO 

CCFFCN/CCN/P TIC, ^IC.AIC, ALFA 

CCMMON/TFPn/ I ,PTl 1 ,M2l ,PT2l, ySl ,PTSl 

CCM'^rN/ NS/*-! 2, PT 12, M22 ,?T22 ,MS 2,P TS2 

CCy'4PN /O IF /M2 , PT2 ,P2 

CCKMfN/CFF/MSN ,M3I, PT3 ,PT3M,HA ,AS 

real 40 ,yiC, Ml I . ,^42 ,M3N,y3I ,M<»,MS ,M6 

REAL L 3 ,N, lSP,y2l ,y22,MSl,MS2 

R E AL '^12 

THIS Section computes expansion from throat of the inlet 

TC POINT 2 

AN=AI*M I 2 

Gl-1.^GA2^MI 2*<2 



M2=mI2 

CAIL CALr‘-(M2. AN,A2,Gl ,GA2 ,GA12) 
IFULTO.GT.l )RETURN 
GC T0( 31 .33, 33 ). I 
prz^PT 12*PI02 
GC TO 35 

F AC=(A2-A I )/ ( A2-AI) 

Pir.?S= I .-FaC^( U-PI02) 

PT2=PT S2 102S 

P2 = PT2* I l.niA2»«2**2»^*(-GA3) 



RETURN 

ENO 



TOA064PO 
TRA06SOO 
TR A0 6S10 
TRAO6520 
TRA06S30 
TRA0SS40 
^RA065S0 
TRA06S60 
TRA06S70 
TRA06S60 
"PAO ^S'^O 
TRA0G600 
TPA0^610 
TR AO 0620’ 
TR 4 3 66 30 
TRA066A0 

TR AO 60 SO 
TRAO 066 O 



TRA06670 
TR A066R0 
TRA066‘^ 
TP A 36^00 
TRA067I0 
^cA0^720 
TR A06730 
TRA0o7A0 
TRAO67S0 
TRAJo760 
TRA06 770 
TR A067S0 
TR \06790 
TRA06P00 
TP A06810 
TRA06BP0 
TRA0 6 3 30 
TRAC63A0 
TO A068S0 
TRA 36360 
TR A0o870 
TR4068BO 
TP A06890 
TRAOoOOO 
TR406910 
TR A06020 
TRA06930 
TRAC69^0 
TRA06950 
TRA06960 
TRA06970 
TRA369H0 
TPAO 6990 
TP AO 7000 
TRAC 7010 
TPAO 7020 
TR40 7030 
TR 40 7090 
TPA070S0 
TRAO 7060 
TR 40 70 70 
TPA3 7080 
TR AO 7090 
TRAO 7100 
TRA071 10 
TR 407120 
TRAO 7130 
TRAO 7190 
TRAO 7150 
TRAO 7160 
TR 40 71 70 
TRAO 7180 
TRAO 7190 
TPAO 7200 
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SL'KRmT INE expan 

CCVMrN/Gen/AReF,AO , a U A 2 *A30.A3,A5# A6,L3f AOAR, ASAP , ILOO» IPR 

CC>^HrN/MP/aA.GAl .GA2,GA12#GAA ,RH0A, T0,P0,U0,'-0 .GRAV 

CQ'^'^HN/FLEL/ohCF ,ETA r, A,N 

CC“vc*J/LCSS/P lOl . P ID?. P Ifl I, P !B 2.P IN 

CPmmpm / rtta/W a, , WT ,F .ROOT ,C ST AR t PT A t TTA 

CpMMfj/GF l/R F , GF ,GH I ,GF2 .GF I 2. GF3 

C0««CN/NCZ/^6, PT6.P6 

CC'-mpn/qeS/CF. Th^lST ,I SP.SFC 

Cc^^'-N/ChG/P T5 ,P5 

CCP«rM/iN{ /p ro.TTQ 

CC^' '-IN /r. C.N/P TIC, '-1C, AlC* ALFA 

CnMMO‘j/T i-RiV Ml I.PTll,M2l.PT2l,^-Sl,PTSl 

CC‘-'-^N/NS/‘-l2, PT 1 2, *PT2 2,'-S2,PTS2 

CCW'-^^N/n IF /B ?, PT2.P? 

CnvBCN /CFc/ ON , - 3 I , P'»’3 ,3T3B, Ma ,AS* 

REAL '-0 ,-lC, '-I l.>«12 . ’^2 ,B3N,.B3I ,MA,N5 ,M6 
5 = AL L 3 .N, ISP. v?i ,y22,MSl,M52 

CCC This SiJRRC'iT l*,E Cr^P^/TcS L05S=S OUE TO EXPANSION INTO COM0USTOR 
b£TA=( ( ‘-2*A2/A3 )<*2 »*( I . ^G A 2 2** 2 ) 

‘-3I=S3RT (( S3R T U , f4.>GA2*.3ET A) -1. )/( 2.»GA2 J ) 

Gl = l.t'GA2*-2**2 

PT3.'--PT2» 1 (I . «-GA2*N3 I**2 )/GL )**GA3 

RETURN 

ENC 

SUPPnjT I RE CH^CK 

Cn^PON/GFn/AREF.AO, Al. A2 ,A30,A 3fA5, A6,L3,A0AR, A5AR. I LOO, I PR 

CC'Mf'-; / A It- /GA, GAl .GA2,GA12 ,GA3 ,RHOA, TO,PO*UO,HO,GRAV 

CCNM^N/FLEL/RHCF.ETAT , A.N 

0 /LOSS/ PI ni,Pin2.PTrtl,Pie2,P!N 

CC^'-^N /BLR/r^ A. «F . «T . F .ROOT ,CST AR. PTA ,TTA 

CCvyON/GFI /R F.GF , GF I , GF2 ,GF I 2 , GF 3 

CC'^Mr>4/f.cZ/B6,PTh,P6 

CC^'-F^j/ses/Ch. THOUST ,ISP,SFC 

CC'^«CN/CFn/PT5 ,P5 

CC'-'-fn/ IM /P TO , TTO 

CGyM,N/cCS/f> T iC.y 1C , AIC,ALFA 

common/ TFO fV Ml I .PTl I ,H2l.PT2l, -MSI ,PT$l 

common /NS/*- l 2. PT 12.M22 ,PT22,M52,PTS2 

CC'^MfN/0 IF/M2. PT2.P2 

C^VMGN/CFF/M3^J ,V3 I, PT3 ,PT3M,H<, ,AS ^ 

real •-0,Mir. ‘-ll.‘-'l2,M2 .M3N.M3I *MA,M5 ,M6 
real L3.N,!SP.M?I.M22.MSI,MS2 
lFiPT3M.CF.PT3 IRfTURN 
II nn= 70 

IFi IPH . LT. 1) PFTUON 
PRINT 33,A5AR,A0AR 

33 format ( IX,2F 7 , 3, 3>. • PT3 IS MOPE THAN MAX- VALUE ALLOWED*) 

IF( IPR.GE.2)CALL PR IM 

RETURN 

ENO 

SUn.Pi ‘UT I NE HEAT 

CCMMrN/r.en/APEF .AO, A l,A2 f A30.A3. A5,A6,L3,A0AR, A5AR * I LOO* I PR 

CCVMf N /A IR/G A.GAl,GA2,GA l2,GAi *RHOA, TO ,P0 * UO • MQ • GRA V 

CQMHi\/ FLEL/RHOF .^TA T, A,N 

COvMfiN/L C5S/ P I Dl ,P 102. P I Bl *P IR 2.P IN 

CC'-mpn/olR/W A. ^F , „T , F .ROOT ,C ST Aft, PTA ,TTA 

CCf-MrN/GFI/RF,GF,GF I *GF2 , GF I 2* GF 3 

COMMfiN /NCZ/M t, PT6,P6 

COM-dN/PES/C F, THRLST , I SP ,5FC 

CQMMn.;/c^.n/p t5 ,P5 

COM'-FN/IM /P TO.TTG 

CQMMrM/rCN/P TIC. Mir , AlC, ALFA 

COMMf.N/ TF = n/Ml l.PTlI,M2l,t^T2l, MSI. PT SI 

CrvMf!/NC/vi2 . PT 12, ‘-22,Pr22,MS2,PTS2 

CCMMr»j /:) if/-^ 2 , PTJ ,P2 

CFMf^f.N /CFF/:^3N ,M M , PT3 ,PT3M,MA ,AS 

PEAL MO , M IC, »-l I,M 12 ,M2 ,Mift,M3I ,M6 

REAL L3 ,N, ISP,M2l,M22,MSl .MS2 



TRA0T210 
TRA0 7220 
TRA0 7230 
TRA072A0 
TRAO 7250 
TRA07260 
^RA0 7270 
TRAO 77FO 
TR AO 7200 
TR ^ 3 7300 
TRAO 7310 
TP AO 7320 
TRAO 7330 
TR A J 73-^0 
TRAO 7350 
TRAO 73RO 
TR AO 73 70 
TRAO 73SO 
trad 7390 
TRAO 7^00 
TPAQ UlO 
TRAO 7A20 
TRAO 7'V30 
Tf AO 7440 
TRAO 7450 
TRAO 74^0 
TRAO 7^*70 
TRAO 74dO 
TRAO 7490 
TRAO 7500 
TRAO 7510 
TRAO 7520 
TRAO 7530 
TRAO 7540 
TRA07550 
TP AO 7560 
TRAO 7570 
TRAO 75dO 
TRAO 7590 
TRAO 7600 
TRAO 7610 
TRAO 7620 
TRAO 7630 
TRAO 7040 
TRAO 7650 
TRAO 7560 
TRAO 7670 
TRAO 75dO 
TRAO 76^0 
TRAO 7 700 
TRAO 7710 
TRAO 7 720 
TR A07730 
TRAO 7 740 
TRAO '7 750 
TRAO 7760 
TRAO 77 70 
TRAO 77dO 
TRAO 7790 
TRAO 7R00 
TRAO 7fllO 
TRAO 7820 
TRAO 7830 
TRAO 7H40 
TRA07.S50 
^RAO 7860 
TRAO 7870 
TP AO 7 880 
TRAO 78^0 
TRAO 7900 
TRAO 7910 
TRAO 7920 
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REAL 

RA=29,31A 
*-4 = 0, 7 

call C ALCH(H4, A5* A3 ,GFl,GF2,GF 12) 

P-SORT ( TT4/T TO )*( 1. ♦GF •MA* ♦ 2 ) / SOR T( GA/GF*RA/RF )* ( U#*F) 

70 0EC=R*«2-4 .♦GA 

I F (OFC .L T,0. JGO TO 74 
m3N=(H-SCRT( DEC) )/( 2*GA) 

I F ( I .G"' . I )G0 TC 72 
m3i-'43n 

0EC={1 .♦(GA-1 )/2.<M31»«?)/(1.^(GF-1)/2.*MA**2) 

8 = B*SJPT (CC-C ) 

1=2 

GO TO 70 

72 *-3N=(M3N + m3i )/2, 

P IP2 = < I ♦(GF- I J /2»M4««2)**(GF/{ GF-l) ) / 
/{U(GA-1)/2^H3N**2)*«( GA/(GA-l ) )•( 1^GA*M3N*»2 )/( 1^GF*H4*»2 ) 
PT3 = PT4/FIP2 
RCtu^N 

74 ILCC=5 

I F ( IPR, LT.D RETtJON 
PRINT 75fA5AR,A0AR 

75 FGF''AT( L)i.2F 7. 3.5X, ‘MEGETIVE ARGUMENT FOR H3M 
IF ( !^R .GE.2)CALL PR IN 

RFT'JPN 

ENO 



ccc 

ccc 



52 



53 

54 

55 



SLBPOUT INF Cno VA L 

CCV'"‘'N/GFC/ARE^»AO. Al,A2 . A 30 * A 3 . A 5, A 6, L3 t AOAR, A5 AR . I LOO 

CC^MHN / a IR/GA.GA I ,GA2* GA 12»GA3 fPHOA, TO . PO * UO t MO , f,R A V 

cc^mpn/flfl/rhcf ,FTAT,A*N 

CCV'^CN/LCSS/Pini ,PIO2,PI01*PIO2»PIN 

CCV'^rN /QL4/W A, wr ,wT ,e,ROOT,CST AR*PT4,TT4 

Cqumo^/GF I/RF ,GF ,GF I ,GF2.GF12. GF3 

C 0 v^'^‘l/.\C :/‘^6 . PTotPb 

CC'-yCN/FFS/CF , THPuST . ISP.SFC 

cr''^'*j/Chn/PT 5 .pF 

CCMMO‘4/ I M /p TO ,TT0 

CCVMnN/CCN/PTlC.‘-lC.AlC. ALFA 

COMmCN/ TKR n/Mll.PTll ,M2l ,RT2l. MSI , pT si 

CCM«'^N/NS/“l 2,PT12»M22.PT2 2 »mS2*PTS2 

CCMMfN/n IF/M 2, PT2.P? 

CCMM^M/Cf-R/M 3N.M31 » PT3 ,pT3M,M4 ,AS 
real Ml I, y 12 .M2 .M3N,'-3l .M4.M5 ,M6 

P E At L 3 » N . I S P . y 2 I . M2 3 , MS I , MS 2 

THIS S'JPFnuTTNF riiCMLATES CORRECT PLACE OF NORMAL SHOCK 
ANC THFRFFORe, CORRECT VALUES CF LOSSES AT EVERY STATION 
ASL=AI 



ASR=A 2 
I AS = 0 

AS = ( ASL ♦ASP) /2. 

I A5=iAS> I 

IF{ IAS .GT.15 )GC TO 58 
call TmpcaK 3.«S1 lAS ) 

CALL NSR (3*MSl »MS2» PTSI.PTS2) 
CALL OIFFUSI 3* AS*MS2) 

CALL =XPAN 

IF( ILOO.LE.I )G0 TO 53 

ILf^= 8 l 

RRTUFN 

0EL = PT3M-PT3 

ACFL = A0S(DEL J/PT3M 

IF (ACFL .LE.O .005 JGO TO 55 

IF(OFL) 54,55.36 

ASR=AS 



GO rn 52 
PT ll = PT SI 
RT 12 =ptS 2 

Ml i^MS I 

mi 2 = *-S 2 
PCTURN 



\ 



TRA07930 
TRA07940 
TRA07950 
^RAO 7960 
TRA0 79 70 
TRAO 79R0 
TR AO 79R0 
TP A08000 
TRAOdOlO 
TRA03020 
’■RAO. -1030 
TRAO 5040 
TRAO0O5O 
TRA 0 T060 



TRAO 5070 
TR AQdORO 
TRAOttOOO 
TOAO'^lOO 
TRAOrtl 10 
TRA03120 
TRAvVdl 30 
TP AO«I40 
TPAOrtlSO 
TRAOdloO 
TRA05170 
TRAO dldO 
TRAO HI 00 
TRA08200 
TRA0«2 10 
TRAO022O 
TRA0d230 
TT AO 0240 
TRAO d?50 
TR AO 8260 
TRAO -12 70 
TRA 0 d 2 dO 
TR AOdPRO 
TRA0d3O0 

TRAO mo 

TRA0d32O 
TRA08330 
TRA0d340 
TP A0d350 
TRA03360 
TRAO 8370 
TRA0d390 
TRAO d IRQ 
TR AO 3400 
TRA0d4l0 
TRA0d4?0 
TRA0d430 
TRAO044O 
TRAO 8450 
TRAO d460 
TR AO 8470 
TRAO 8480 
TRAO d4Q0 
TRAO05OO 
TR AOo5lO 
TR A08520 
TRA0d530 
TRA08540 
TR AO 8650 
TRAO 8560 
TRAO 3570 
TPA03580 
TR A08590 
TRA08600 
TRA08610 
^R AO .16^0 



"RAO 8630 
TRAO 8640 
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56 

58 

59 



llO 

120 

130 

lAO 

150 

160 

165 

170 

180 

190 

200 



ASL=AS 
CC TO 52 
ILCO=92 

I F I IPO ,LT.I) octuRN 
PRINT 59,A5AO.A0AR 

FOP-AT I IX, 2F 7.3* »OOES NOT ^ INO COPRECT NORMAL SHOCKM 
IF I IPO •GF,2)CALL ?RIN 
Rf TORN 
6NC 



SUPR*"HT INF POT N 

COMMf N/npn/A PPF , AO, A 1, A2 ,A30,A 3, A5, A6,L3, AOAP, A5AR, I LOO, I PR 
CC'-MfN / A IQ/GA,GA1 ,GA2,GA12 ,GA3 ,RHOA, TO , PO , UO, M 0, GR A V 
FtFl./RHGF ,FTAi", A,N 
COvMfiN/LCSS/PI 01 ,PI 02,PT9l,PTB?.?IN 

/PLP/W A, hF ,«T .F.ROOT ,CST AR, PTA ,TT4 
CGVMrM/GFI/RF,GF ,GF 1,GF2,GF12,GF3 
C0M*^nN/NC2/M6. PT6.P6 
CO'^-'r'N/PFS/Ch, THRLST . I SP ,SFC 
CrvM-]r,/Chn/oT5,-^5 
CO'-'^rN/ I M /P TO , TTO 
CC»'*-*f^N/CCN/PT AlC, ALFA 

CCv?n' / Thon/ Ml l,PTi L ,M2 1 ,pT 2 I, MSI ,PT SI 
Cf^'-Mr.N/NS/Mi 2 , PT 12, M22 ,PT22 ,MS2,PTS2 
COMMON /O if/m 2, PT2 ,P2 
CCM.MrN/CH«:/M3N I , PT3 ,PT3M,MA ,AS 
RFAl MO,Mir, -U,M12 ,12 ,M3M,M3I ,MA,M5 ,M6 . 

PFAL L 3 ,N, is 3, M2 I ,.M22,MSl,MS2 
PP. INT l I C ,wA , F, „T, F ,k 00T . CSTAR, PTA 
PRINT 12C.M6 ,P^6,P6 
PRINT l3C,CF,ThOIJST,ISP,SFC 
PRINT IAC,PT5,P5,P0 
PRINT 1 5C,P0 . TC,MTO. TTO 
INT 16C,"F .GF.GA, T^A 
PRIN'T 16 = ,PT 10 ,M ir. Air.,ALFA 
PRINT 1 7C.MI 1 , PM l,.M>l ,PT2l,MSl,PTSl 
PR INT 18C,M12,PT12,M22,PT2 2,MS2,PTS2 
PRINT 19C,M2 ,PT2,F2 
PRINT 20C,'-3N, M3 I,PT3, PT3‘-,MA, AS 
FORMATI 2 >, A, «F ,^T ,F,?(r)nT,CST AP, RTA-* ,7EU*A) 
format I 2 X, »M6, PTt,P6*» ,3F 1 1 .Ai 
FaPMATI2X,»CF,F. i<P,SFC*».AFll .A) 

F3PMAT I2X,»PT5,P5.^0=»,3FIUA) 

FORMAT ( 2X, »PU. TC, c>TO, FTO = ’ ,AFl l.A) 

FOR'-'A! { 2 X, 'R F , gf , TTa* » .AFll .A) 

F0PMATI2X,'PT1C,''IC,A1C,ALFA = »,AFII.A) 

FORMAT ( 2X . »M I 1 ,PTll .M2l,PT2l,MSl,PrS 1*«,6FIUA) 

FORMAT! 2 X, 'M 12 , 12 . '"2 2 , PT 2 2 . M S2 , PTS 2* • , 6E U *A ) 

FCRMATI2X,* '^2 ,PT2, P?= ’ , 3F ll .A ) 

format ( 2X, »M3N,M3 I, PT3 .PT-iM,M<. ,AS=» , 6FU ,A » 

RFTIJCN 

FNO 



SUBRCUTINF PRlLQOdPHINI 

CC'^-ON/GFO/APEF. AO, A I , A2 , i 30 , A 3, A 5, A 6, L3 , AOAR , A5AR , I LOO, IPP 

CCmm'^N / A IR/G A,GA I ,GA2,uA 12 ,GA3 ,RHO>, TO,PO,UO,MO,GRAV 

COMMON/ FLFL/ PHOF ,FTAT, A, N 

CC««r‘|/UCSS/Pini ,PI02,PIBl ,PIB2,PIN 

CCMMfiN/bLR/w A, wF.^T ,F,R00T ,CSTAR,PFa ,TTA 

COMMrN/GFI/PF,r,F,GF 1,GF2.GF12, GF3 

C0MM^N/NCZ/M6, PTo,P6 

CCMMrN/RFS/CF, THRUST, ISP, SFC 

CCWM^,^J/r>^/PT5 ,P5 

CCMMCM/ 1 M /PTO ,TTO 

CCMMPN/rCN/P T1C,MIC, AlC, ALFA 

common/ T HR 0/ Ml l,PTll,.M2l,PT2l,MSl,PTSl 

COmmci'i/nS/mi 2 , PT12,M22,PT2 2 ,MS2,PTS2 

COVMf)N/n IF/M2, PT2 ,P2 

rOMMr N /CHF/M 3N ,M3I , PT3 ,PT3M, MA ,AS 

COMMCN/TPA/P I ,MUA ,ORAG ,RP° ,LPP , WPR, U , W0, OFLU, IFLT , TT ,TFTA, IRAM, 

, IL ,ILO. Y,YO. T0 9,Xl.Yl,X2,V2,X3,Y3,WPPB,P,T0Ra, I TR A , TETAOO, ^CB , YOB 
CCMMCN/OPG/TFTP,CCN , C F T , COWW ,C OWF , A P R, SPR , SWW, Q,XMO 



TRA08650 
TRA08660 
TRA08670 
TRA036B0 
TR A0B690 
TRA08 700 
TO A 0 37 10 
TRAO 3720 
TRAO 3730 
TRA087A0 
TRA08750 
TR AOb 7A0 



TP A08770 
TRAO -17H0 
TRAO 3790 
TRAO 3800 
TP A0M310 
TRA03620 
^RA03830 
TRAOB^AO 
TR A03350 
TRA03360 
TP AO 3H70 
TO A0drt‘>0 
TRA03B90 
TP A 0 8900 
TR A0891 0 
TP AO 3920 
TRA03930 
TP AO 39 AO 
TRAO 3950 
TR AQ3960 
TPAO 3970 
TRAU8980 
TRAO. 1990 
TSA09Q00 
TRA09010 
TRA09020 
TRA09Q30 
TRAO90Av) 
TPA09050 
TRA09060 
TRA09070 
'^RAO 9090 
TRAO 9090 
TRAO 9 1 00 
TRA09U0 
TRA09120 
TQA09130 
TR AO 91A0 
TRA09150 
TRA09160 
TRA09170 
TR A09190 
TRA09190 
TRA09200 
TRA0W210 
TRA09220 
TRA09230 
TR A092A0 
TRA09250 
TR A09260 
TP A09270 
TRA092P0 
TP A09290 
TP A09300 
TRA09310 
’R A09320 
'TO A09330 
TRA093A0 
TRAO 9 350 
TRA09360 
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fortran 



MAVAL POSTGRADUATE SCHOOL 



peal MO.J'lC.;jU,yi2,M2.M3N.M3I TPAQ9370 

REAL L3.N.^P,H2l*M22.MSl,MS2*MUA,LPP TRA093R0 

/AX T'»A093P0 

A ICAO- -» I C/AO TRA09400 

AlAO^Al/AO TRA09410 

A2AO=A2/AO TRA09420 

A3AR=A3/ARPP TPA09430 

A5AR=45/ ARFF TRAO^Uo 

A6A5=A6/A5 TOA09450 

ASAG^AS/AO TPA09460 

PlCRO^PTir/PTO TRA09470 

RIIP1C=PTII/PTIC TRA09480 

P12RU = ?T12/PTII TPA09490 

P2R12= PT2/OT 12 TR4 J9SOO 

P3P2 = ^r-^../PT2 TRa095l0 

P3KP3=PT 3M/PT3 TRAOOS20 

P4P3=PT4/?n TRA34530 

P6F9-PT6/PT4 T0A09540 

PfcP0=PT6/PT0 TR\09SS0 

«FWA=F*lCO. TRA09S60 

IF{ IPP IN.JT.2JCO TO 77 TRA09S70 

I PR ("*- 7 TRAO ^5flO 

PRINT 21C ^0A09590 

PR IN’' 2 I 2 •' , TR AO >600 

PRINT 2l4 TRAO/610 

PRINT216 TRA 09620 

PRINT 218 tpao9*j3q 

RRINT 22C TRAO‘,640 

TR(NT 222,AREF.L3.TETP TRA0O6S0 

PRINT 224 TRAO 9660 

RRINT 226 .A0AR ,AlCAOfA lA0fA2A0 *A3AR , ASAt ,A6A5 TRA09670 

RR(NT 232 TRAQ96R0 

PRINT 234.M0 *MIC.H6 TRA0060Q 

PRINT 236 TRAO‘9700 

PRINT 2 38tP{0l,PI02iPIN TPA00710 

PR INT 244 TRAO i 770 

PRINT 246 • TR^09r30 

PRINT 24E,PO,TO,PHOA ,PTO.TTO,GA- TOA0‘/740 

PR INT 75C • TP AO )760 

PRINT 252 TR^00760 

77 CC.NTINIJE TPA09770 

PRINT 2 56 ,TI ,M0, ASA0»WA,WFWA,M2,H3N| >-3It M4, ^ TR\0S7B0 

• PlCPOf P l2Pll.P3P2.P4P3tP6P0* ^ T9A00790 

• GF .TTh.CF .thrust, ISP TRAO >800 

210 format I IFI,4 0X,*CCCCCCCCCCCCCCCCCCCCCCCCCCCCC» ) TRAO 98 10 

212 FOPm-t I 4U,» CCC CCCM TRAf)9B70 

214 FORMAT I 4 IX .* CCC SOLID FUEL RAMJET CCCM TPA09830 

216 FCFmaTI 41X,'CCC CCCM ’^RAu9840 

218 format I 4 1X,’ CCCCCCCCCCCCCCCCCC CCCCCCCCCCC* ///) TRAO 9850 

220 format I 5 X . *OFp>-ctpI cal DATA:'//) TRAO9860 

222 Fjp4ATl8X,*A«fF-M"lU4, ♦M2*.3x,*l3»»,Ell.4,»M« ,5X,»TFTP«* ,F5.l//)TRA0 >5 70 
224 FORMATIICX,»AO/AFCPI, 'AIC/A0M3X, »Ai/A0M5X, ‘A 2/ AO* , TRAO 9880 

*5X,'Ai/ACFF',^X,'A5/AKEF*,4X,»A6/A5*/) TRAO 9890 

226 format I 6X, 7F 1 0.4/// J TRA09900 

232 format ( ICX, * « 0» ,dX, *MIC* , 7X, 'M6*/ ) tra09910 

234 format I 6>, 3F 10-4/// ) TRA09920 

236 FORMAT ( 5X. 'CONSTANT LOSSES:*/) TpAoqv^30 

238 format IdX,'PI0l*MF7.3,jX,'PI0 2=»,F7.3,5X,'Pn»MF7.3//) IRA 09940 

244 format I 5 X, M NI TI AL FLIGHT CONO I TI ONS : * // 1 T9A09950 

246 FCR''ATIlCX,»P0I<r,/M2)»,2X,'T0lK)',6X,'R00IKG/M3)*.lX, TRAO 9060 

• * PTOING / »^2 i* . I X, ' ttOI K ) * ,<iX, 'G A* / J Tt^ A O 90 70 

248 FCRMAr I 8 x,5c ll .3 .F10.3///1 TPA0^980 

250 FOPMaTI3X,'TImc.,7x,'MO*, 3X, 'A 5/ AO*, IX,' WA» ,4X, 'WF/WA' ,lXt TRAOPORQ 

• * M2* , 9X , 'M3N ' . 3X, *M3 I * , 3 X, ' M4' , 7X , ' T CT AL RRtS. RATIOS*) TP.AIJOOO 

252 format I 5 6X,' IC/C* ,7X , • 12 / ll * »2X, • 3/2 * . 3X, TR a I 3010 

*»4/3‘ , 2X , *6/0* ,3X,* GF* ,3X, »TTmK) ', IX, TR A I 3020 

»*CF' ,^X , 'FIN ) ' ,3X, * I SP'/ ) TRAl 30 30 

256 cr)PMATIlx,F7.2.F5.2,l2F6.3,P5.2.F7.l,F6.3,F7.l,F7.l) ^PA1 30 40 

PFTUPN TRA10050 

FNC TRAl 0060 

TRAl 00 70 

SUflPCUtI6F TRAJ TRA10030 
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CCC 



131 

132 

133 

135 

137 

140 

125 



23 



138 



COMMON /TP A/P I ,MUA*ORAG.RPR.LPo ,WPR , IT, WBi OELU , 3FLT . T ! , T ET A , I R AM . TRAIOO^O 
, IL.ILO. Y,V0.T0 8.Xl,rUX2,Y2,X3 *Y3,WPR9,R, TORB, I TRA • T ET AOO , X CB t YOB TPAlOlOO 

CC''MrN/GEr;/ARPF, AO.Al, A2 ,A30,Al,AS,46*L3,A0Afl, 45A«, ILOOt I PR TBAIOUO 

CCMMr^/A iR/3 A, GAl,GA2.G4l2,GA3 tPHOA t TO.PO , UO fHO . GR AV TPAl 0120 

COMMON /3LR/W A, v^F,»,T , F . ROO T , CST AR , P f 4 , TT4 TRA10130 

/ RES/C*" . THRUST.I SR ,SFC TRAU140 

CO-MPN/CCN/R nC.MlC.AlC.AL = A TRA10150 

CCMMCN/ORG/TFTP.CCN .CFT, Cnww.COWF ,AP R,SPR,S^W,0,XM0 TPM 0160 

RF AL >-0 L 3,«UA. rSP .MIC TRA10170 

A!TRA=(TPA TSAlOiaO 

AI TDA=AB5( AI TB A) TR4lin90 

(F( At*D A .GT. 1. )G0 TO 23 TRA1 0200 

IF (ITPA.GT.J )1TRA*2 TRA10210 

IFUTRA.IT.O )I TRA=-2 TDA1 0220 

TI=0. TRA10230 

Y0*0. TRA10240 

Y=YO TRA10250 

TFTAO=TETAOD*P!/ ISO . '^RA1 0260 

TFTAsTfTAO TRA10270 

VACUUM TIME OF IGHT( TOFVI : TRAIJPBO 

T0FV=2 . ♦li*5I N( TETA) / GP AV TRAU32<50 

ILC=4 ^RA10300 

CELr=T0FV/50./lL0 TPA10310 

I F (ITS a.LT.O JGC TC 125 . TRA1 0320 

-RITF(2,13U TRA10330 

WRITC(M131) TRA10340 

FCRMATt IHI///,40X, 'RAMJET TR AJ ECTOR Y V // .) TRMO350 

u 9 \rci 132) TRM0J60 

FHFmat ( 43X,» (0»AG COEF. )•///) TRA10370 

•^R !Tc ( 2 , 133) TRAIJ3H0 

WRITFI 3, 133) TRAI 0190 

FOP'- AT ( 4X, 'LPR • . 7X. »WP«' ,7X. *4 30* ,7X, • AO/AR* ,5X, • A5/AR» ,5X, *L3* . T»A I Q4.00 
.8X,MI0' .ex.'U* .9X,»WB» ,8X, 'TOF V*/) TPAl )410 

wPITr ( 7, 135) LPR.WPR ,A30, A0AR,ASAR,L3 tUOfU.WB.TOFV TRAI 0^2 0 

^RlTt^(3,l3 5)L?R.«PR. A3 0 . AOAR , A 5 AR , L3 .00. U.WB.TQFV TRAI 0430 

Fro*‘*AT ( 2>, 13FIC.3///) TPA10440 

wPITri2.l37) • TRA1 0450 

FCPMAT(4X,*Tl',8X,'XV,aX,*Y3‘.ax,'TETA»,6X.'M0»,4X.'P0M8X, TRAI 0460 

, 'PHrjA' ,6X,'T0' .8X, 'MUA • , 7X . ' WP P ' , 3X , 'ORAC* .2X* 'THRUST' /) TRAI J4 70 

«RITC(3,14C) TRA10A80 

cnPMAT ( 4>, 'T ETP* ,6X , 'CON* , 7X* ' COS ' , 7 X, 'COWW , 6 X i ' C OWF' , 6 X , TRAI ja 90 

, »APR',7X,*SPR' .7X,'SWw',7X. 'G* .9X,’XM0'/» ^RAIOSOO 

IL=ILO ^ TPA10510 

X1=0, TRA10520 

Yl^YO TRA10530 

X 2 =u*rns (TET A) *0EIT ^PA1 0540 

y 2='J*S I N (TET A ) *OELT TRA10550 

CALL ATM TRAI 0560 

MO=U/SCRT(GA ♦G?AV*R*TO) TRAI 0570 

TRA10580 

AIR offence CA 5E:0Q NOT LET **0 BE TCO SMALL. TO ALLOW MANUVEP INGTRAl 0590 

TRA I 060G 

IF( IRAM.GT.DGC TO 138 

IF( *'O.LT. XMO)GO TO 29 
CALL ORAGG 
IF( TLGO.GT.l )RETURN 

IF YCU taANT 0RAG==0 OR THRU5T»0RAG CASE. SPECIFY THAT HERE* 



139 

141 

127 

23 



CALL OYNA 
Y = V 3/0. 3C48 
TET0«TETA*180./P I 
IF( YB.L*" .YU) (L=r LO 
IF( IL.LT.ILO )G0 TO 28 
IL=0 

IFMTRA.lt. 0)GO TO 127 

WRI TP (2 . 139) TI ,X3,Y3,TETO,MO.PO.RHOAi TO.MUA,W»R,ORAG.THRUST 
FOKMATI ni.^E 10.3.F6.3.4E10.3.3F7. 1) 

wRITFI 3 , lAl) TPTP.CnN.CFT.Cnww.CnwF.A PR,SPR.SWW,0,XM0 
FcJP‘^AT(IX.F5.2.5X,o£10.3.2X.F5.2) 

IC(V3.L£.YOM.O to 29 

!L=fL^l 

RFTURN 



TRAI 0610 
TRAI 0620 
TRAI 0630 
TRAI 0o40 
TRAI A650 
’■'5 M 0660 
TRAI 0670 
TRAI 0680 
TR A I 0690 
TRM 0700 
TRAI 0710 
TRA1 0720 
TRAI 0730 
TR M 0 740 
'"RAl 0750 
TRAI 0 760 
TR A I : 7 70 
^RAl 0780 
TP A I 0790 
^RAioaoo 
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29 tFCTrpft*GT.O.OOI )Cn TO 30 

TnDR=TI-Of*LT 

xnftsx3 , 

ynR 3 Y 3 

<F| I P AM.NF, nr.n to 30 

TPAm^2 

cn TO 138 

30 IF (ITRA.LT.O 1 r,n 31 

PRINT l<» 6 »TETA 00 ,*^ 0 RR,X 0 fl,YnR 
146 *=OPNM ( 1 ){,///, lOx ,» T^TA= », F5.^ . 7X ,» T f Mp OF RlJRNlNG-*, 

,F6.2. I X , •';cC * r 7X, »P ANOE OF fTU' N ING» • , 10 .4, I X, • KM S 7X, 
cf BORNlNGs • ,F 10*^ , I Xt 'KMM 

<;Tnp 

31 PPpi^ 148 ,A 0 AP ,A 5 AP,T'^RH,XQfl,YnR,Tl, X 3 tY 3 ,TET’>,TFTAOO 
148 F 0 PM'.T { 2 X, 10^1 1 . 3 ) 

Tinn=?r) 

return 

END 

<;UPR»>UT I NF A 

CQP.T.., /tRA/P ! , MU A, “T rag po , l PF , wP© , U ,Wf^,'3FLU,A©LT,*f,’rTA,lRAM, 

• IL. n 0, Y,Y0, infl.xi.Yl ,X2,Y2*X3.Y3.-PRn,R,TTRM, 1 toa, T^T^OO. xrn , yob 
cn»^Mnu/r.Fn/AP©F , '.0, n ,a? .A^^OtA 3 »a5.a6, l3»aoac. a^a; . m nn,(pp 
r rwMf'*. / 4 IP/G A, GA I ,GA 3. GU.’,GA3 ,yHDA* TO, P 0 , 00, ‘^0 , G' A V 
rfal *^0 , 1 pp, I 3 , '^u ' , I <;p , Mr 

CCC AT^•^CPHCPl^ fmpmmla ©nq oc^ckSuRP, ofn S’ TY mps a r ior ,vt sc'^^TTY 

CCC UNl^K IN KG/ m 2 ,XG/m 3 , TFG. K , KG/ M.G'-G R T . ; lE I G H T ( Y I , 1 M FT, 

PO = l. 3 3 3 22 P(K*©XP{- lh. 2 -—G 6 *Y<‘« 1.091 
RNGA- I .? 2484 5 « EXP I- 7 . '♦E - 06 *Y**^ 1 . 15 ) 

IF! Y.GF. 36030)00 TO ?0 

TO = ?R« . 1 ttFXP( - 3. 70 7 ‘’ 4 ‘^- 06 »Y** 1 .0709 1 
MHA=l . 7 0 3R-05©Fxr {• l5.523<’-06^y^*0.89B4) 

CP TP 7 7 

25 TO=n/.?A 

«UA=l,4 1 24C-05 
27 C0NrT‘ii)<= 

RFTIIPN 

END 

SUPOPMT I NE 8QQS 

C^MMO^J / TRA/P I , MUA.DRAG ,RPO , LPP ,WPR.U , ^ 8, OFL!) , "*F I T,TI,rcTft, IRAN. 

, IL . Tl 0 , Y ,Y 0 , TOB, XT , Y T • XP, Y 2 , X** , Y 3 ,WPR 0 ,«, TORH, I TRA , T©TA 0 O, XC 8 , YOp 
Cr^MfM*/nEn/APcp, AO, M.A 2 .i 30 ,A 3 ,A 5 .A 6 ,L 3 .A 0 \c,A 5 \c, {LOO, 

CCWMr»j / A IP/G A, GAl ,r.A2 ,GA 12 ,GA3 TO , PO , HO , MQ . G’’ A V 

CDMMPM/BLR/W A, wF,mT,F ,R nnT,r STAR.PT4 , TT4 

CQMMn*. /pFS/FF , TH^HET , 1 SP ,<^FG 

COMMON /CPN/P TlCfMlC ,AIC.AI FA 

RFAL MO,lPP,l 3,MilA,IGP.Mir 

RHPB*! .65F03 

I SPB«P40 

VB*< A30- AS)»l 3 

wa-vp«RHrp. 

UFQ- 1 SPB*GPAV 

OEl UsuEG^ALQG (Wpc / < MPR-W 8 1 ) 

PPT(Jh^ 

©NO 

^URF^UTlNF ORAGG 

CnMMPM PA/© T , mha, OF AG ,? H '• ,LP ,WP - • U , W6, OFL J, T,TT,TPTA, {PAM, 

, I L, ILO. Y ,Y0, ©, XI , Yl, X2 , Y2 . X3 f Y3, MPP6 ,3 I TP A.tct^OO, XOrt, YPB 

POMMom /GC n/AREF,AO, A 1 , A2 , A 30 i A 3 1 4 S . A 6, L3 , A OA © . A S A© , 1 1 pO, I ©R 
CC^mom/A It /GA, GAl ,GA2, GA12 ,GA3 ,PHPA, TO , PO , U0» '^O, G© A V 
rOMMON/HlP/WA, WF ,WT F. qnm ,C STAR t©T‘ 4, TT4 
COMMPN/ PES/C ©, THKilSt , ! 5P ,^FC 
Cnyon.i/r^M/PTl C, *^IC , Al C, Al FA 

COMMON /DPG/TF TP, r''N ,r.F T,GPWW,C OWF ,APH, SPR, SWW, 0 • XMQ 
PEAl .LPP. F3 r .'^IC 

I© ( ( P AM.GF . I )f,"' T » ^0 
CALI OO ArMir FTP, “1 ) 

IFMI P'^.LF.l )GC TF c,., 

ILr'F, = 2 2 

RCT|)F*>| 



TRA10810 
’'PA 1 0R2O 
■frAI 3S30 
TRAl 0840 
TRM 0350 
FQ a I 9360 
TPAT 1870 
TO A1 OPaO 
TR ^ I )h90 
1900 
FP M 191 0 
^©11 19^0 

TOM 1070 
FPM ,)M4 0 
TRA^ )050 
FP41 )Q60 
TR •. I 0 9 70 
FRn Tiao 

’■R^ 1 0990 
^RA 1 ! 000 
FRM 101 J 
FOM 1020 
"R'\’ ’030 
^RA I 1040 
fr \ • ! O«^0 
tram 060 

TPM 1 070 
F3 Al 10RO 
TPAl 1090 
TRAl IICO 

TQM n 1 0 

FRM 1120 
TRA 1 1 130 
TRA I 1140 
FOAl 1150 
M60 
’*© M 1170 
TRM 1130 
TOA’ 1190 
1200 
^©4 1 1210 
TQA I 1220 
TRA I I 230 
TR I ; “>40 
TR.AI ! ?60 
TO ^ 1 1260 

FRA I 17 70 
tram ?30 
fr a I 1700 
fq a I I 300 
TP A M 31 0 
""PM I 3'>0 
FP Ml 330 
FRA I I 340 

TRM 1 7S0 

TR M I 360 
'’■©M 13 70 
foam 3©0 
^R M 1 390 
TRM ’ 4'^.'! 
TPAl 1410 

toa I 

TQM »4 30 
^0 A I 1 440 
TR\1 1450 
fqa I 1 460 
TPM 1470- 
’’Q M ’ 4-50 
’■SMI 40.1 
TRAl I 5'^d 
r R i I 1510 
FR M 1 520 
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cc 

c 

c 



c 



A5 call DRAGWW{ mo,GA,8WING»CWING,CDWW.)^ 

CALL CALC0{LRR.CFT» 
call CALCO{CwING»COWF) 

APRsPI *R 2 

Sp£?-2» PI *RPR *lRP 
Sk>* = b»RVi ING*Crt ING 
0 = 0. ^)*R HCA^'J 

00 AG=g* ( APO*C0N>SPR*CFT^SWW*(C0 WW^k:0WF ) ) 

RE TURN 

50 CALL CPAGPR 

petufn 

ENO 

SLPOruT inf OC aGPP 

COMMON /TOA/P I ♦ mua,ORAG ,RP»,LPO fWPR.UjWB, OELU*OELT,Tl ,TFTA, IRAN, 

, IL.ILO. r ,rO, TPB, At, Y l,X2f V2f X3 ,Y3,WpO8*R,T0Rtt, I TR a , TETAOO • X CB , YOB 
common /GEO /A REF, AO, Al,A2 , A 30 , A 3 , Ab , A 6 , L3 , AOAR , A5AR, I LOO, IPR 
CC'-MCN / A I0/GA,GAl,GA2,GA 12,GA3 ,»HOA, TO . PO , UO t MO , GRA V 
CCMM(jn /rtL® a, kF , WT , F ,3nOT ,CST AR, PTA ,TTA 

rnMMPN /oes/CF , thpust , i sp ,sfc 

CnM'-O.N/CCN/PTlC,.MlC , AlC, ALFA 

COMMrN/OPG/TETP.CCn,CFT,Cnww,COWF ,APP,SPR,S.^W,0 .XMO 
real mo . LPB, L3.MUA, ISP ,MIC 

COwF ^C. 

CDN = ( C.08 3«-O.CS6/MO**2) ♦( AL fa/ 10. )**l .69 
C0'3=< C.68 3 7-0.3l6b»M0^0.052 5*=^MO**2)*(2./PI» 

CALL CALCOI LPR.CFT) 

A pp » P 1 »po Q * *2 

SPP = 2 .*PI *0 FR*LPR 

0=0. 5*»H3 A<'U*«2 

ORAG = c*(apr *(CCN>CD8) <-SPR*CFT) 

BCTUO N 
END 

SUeorUTlNE OPAGN (TFTP, XM,CON ) 

CO?-MON/Gen/APFF,AO, Al, a 2 *A30tA 3,A5,A6,L3tAOAR,A5AR,ILOOfIPR* 

IFiXM .GF.1.5)G0 to 13 
CCN=( 0 .Cd3>0.C9F/XM**2) *1.21 

aTImo ,lT. 1. 5 J PPPGkA.M CO'nL, which CALCULATES CON.OOES NOT WORK. 
THE APPRCXI ‘-ATE VALUE HF«E FITS AT THE 80UN0ARY. 

THE APPPOX. IS EXCEPTED BECAUSE DRAG IS SMALL ANYHOW* ^ 

RETURN ^ 

13 CALL INTER(^*TETP,XM,CON) 

RETURN 

END 

SUBOOUT IKE 0»AGWW(XM,GA, BW INC, CWING# CDWW ) 

PI =^. * ATANI I . ) 

RACFG*P I /10O . 

DATA T^C/S./ ,0PEF/5./,rt/O.O9525/,C/O .0635/,T/0. 01/ 

AREF = Pl* IDREF*O.C25^)**2/'». 

8WING=8 

C«ING=C 

ALL LfcNGTHS IN METEPS. 

IFIxm.GE .l.2 5iGO TO 13 
C0WW=O. 017 

THF SL0BOUTINE DOES NOT WORK AT I M 0. LT . 1 . 2 5 ) 

the APPFOX. is EXCEPTED BECAUSE DRAG IS SMALL ANYHOW* 

RET'JPN 

13 T^-^W0*BADEG 
Xf*l=XM 
XH2=XM 

CALL PPANT(XM,GA,ANIOJ 

AN I 1= AN I C-TW 

ANI2 = ANIC^TW 

CALL pc FSIXM ,GA,P0PT) 

CALL AP PANT! XM I , GA, ANI U 
CALL PRFS (XMl ,CA ,PIPT) 
call apP ANTI XM2.GA, ANI2) 

CALL, P? ESI XM 2. GA , P2PU 



TRAU530 
TRAl 1540 
TRAl 1550 
TRAl 1560 
TPAI 1570 
TRAl 1580 
TRAl 1590 
TRAl 1600 
TRAl L610 
TRAl 1620 
TRAl 1630 
^RAl 1640 
TRAl 1650 
TRA I 1660 
TRAl L670 
TRAl I 680 
TRAl 1690 
TRAl I 700 
TRAl 1710 
TRAl I 720 
TPAI 1730 
TRA I 1740 
TPAI 1750 
TRM 1760 
TRAl I 770 
TRAl I 780 
TR M 17«0 
TRAl 18D0 
TRA I IHIO 
TRAl 1820 
TRAl 1830 
TRAl 1840 
TRAl 1850 
TRAl I860 
TRAl 1870 
TPM 1880 
TBA I 1890 
TRA I 1^00 
TRAl 1910 
TRAl 1920 
TRAl L930 
TRAl 19^,0 
TPM 1950 
TRAl I960 
TPAI 1970 
TRAl 1980 
TRAl 1990 
TRAl 20C0 
TRAl 2010 
TRAl 2020 
TRAl 2030 
TRAl 2040 
TRAl 2050 
TPAI 2060 
TRAl 2070 
TRM 2000 
TRAl 2090 
TRAl 2100 
TPAI 2110 
TRAl 2120 
TRAl 2130 
TRAl 2140 
TPAI 2150 
TPM2160 
TRA12170 
TPAI 2100 
TR A I 2190 
TBAI 2200 
TRAL 2210 
TPAI 2220 
TRAl 2230 
TRAL 2240 
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8P*8-C/2/SOR T< 

CCWw* ( 2 ./<GA *XH**2l )*( PlPT/P0PT-P2PT/P0PT)o(T*8P/AR6F) 

»6TU»N 

ENO 

SLRPnuTlNE OOANT (XMI ,GA, AN I) 

BETAsSOft T( I •♦2- 1. ) 
r,ARATsSCFT<(GA-l.J/(GA^l.)) 

AM*ATAN<GARAT*dETA> /GAR AT- AT A N( BETA ) 

RETURN 

ENO 

subroutine APRANT (XH I ,ga ,ANI ) 

I F (XMl . LE.UO iGO TO S58 
IF(ANI. LT.O.OGO to 558 
0£TA=SORT (XM I **2- U ] 

GARAT=SQRTM GA-l . J / < GA^l . ) ) 

IBETArO 

552 FsaTA^M GARAT*8ET AJ/GAi; AT-ATAN( 8ETA)- AMI 

FP=l./( I.MGARAT*hETA)**2J-1./ (U>8ETA**2} 

I8ETA=I8eTA>i 
IE ( IBE TA.CT. 121GC TO 558 
IFlFP.ro.c.O JGC TO 558 
BE TAN-P ETA-F /FP 

IF( (A8S(BETAM-8ETA|) .LE. U.E-05))G0 TO 556 
BETArABS (bETAN) 

GC TP 552 

556 x'-tsSQR T (BET A**2»*l. ) 

RFTiJPN 

558 print 56C. beta, am . F,FR .BETA 

560 FORMAT(» ’.'STCP FRON 560* • I 3, A( 2 X, E ll .4 )) 

STOP 

ENO 

SUBPQtJTINE PRES(XmI ,GA.P IPT) 

P IPT = ( I . MGA-l . )/2. ) **( -GA/(GA-l. ) ) 

RETtIPN 

END 

SUBRPtJTTNE CALCrUXL.X) 

CCHMON / TRA/P I , muA.ORAG.RPR.LPR *i<PRtU , wfl, OELU, DELT , T I , T 6T A , IRAN, 

, IL .11 0. Y .YO.TPR.Xl, Yl,X2tY2,X3,Y3 , WP R3 ,R , TOR B , I TRA , TETAQU, X CB, YOB 
CGV*-ON/GEn/AREF, AO. Al, A2 . A 30 . A 3 , A 5 . A 6, L3 t AOAR , A5AR . I LOO. I PR 
CC'-NPN/ AIR/G A.GAl ,GA2,GA12 *GA3 ,RHOA, TO ,?0 ,U0 tNO ,GRAV 
CC'^NON /HLR/W A, .r^T . c,o-)OT .C ST AR , PT<^ , TT4 
CC'^NON /WES/CF . T»^ctST , I SP . SFC 
CCm«niN/CCN/P T IC.'-IC .aic.alfa 
real *^0 ,IPR, L3 .^^UA, I SP ,N IC 
E=EXP( I • ) 

I X*J 

EPS=l.E-C7 

X=2.F-03 

R EL = PHnA»U*XL/NUA 
IF (PFL . LE.2.0 = Co)G0 TO 3 
7 FX = SOF T ( X)*ALOGIO(X*RE U -0.242 

FxoOTs I .C/SORT (x ) •( ALOG 10 < SORT (X*REL )*E) ) 

XNE.<-X-F >/FxOOT 
01 FiABS ( )i-XNE«* ) 

IF(0IF,LE.EPS)R6TURN 
X- XNFW 
I X=*IX^ I 

IF( IX.GT.13)G0 TO 9 
GC TO 7 

3 X- 1.328/ SORT (REL ) 

RETURN 
9 iLnn^20 

u crpBAT (5 Jf *00 NOT FIND A SOLUTION TO CFT; X* • , E 1 0. 3 * 5X , * 0 T F *• * 

, F 10. 3) 

STPo 
ENC ' 



TR A12250 
TRAl 2260 
'"RAl 2270 
TOAl 2280 
TRAl 2290 
TRAl 2300 
■^RAl R'^IO 
TRA I 2320 
TRAl 2330 
TRAl 23AO 
TRAl 2350 
TRA I 2360 
TRM 2370 
TRAl 2380 
"RAl 2390 
TR A1 2^00 
YR A 12410 
YR A I ?4?0 
TR A I 2430 
TRAl 2440 
TR A 1*2450 
TR A I 246 0 
TRAl >4 70 
T? A I 24H0 
TPA I 2490 
TRAl 2500 
TR A I 2510 
TRAl 2520 
’■RAl 2530 
TR A I 2540 
^RA) 2550 
TRAl 2560 
TRAl 2570 
TRAl 2580 
TRAl 2590 
TPAI 2600 
TR A I 2610 
TRAl 2620 
TPAI 2630 
TR A I >640 
Al 2650 
TRAl 2660 
TR.Al 2670 
TR A' 2680 
TR A I 2690 
TRAl 2700 
TRAl 2710 
TRAl 2720 
TRAl 2730 
TRA I 2 740 
TRAl 2750 
TRAl 2760 
yRA I 2 ^70 
TR Al 2 7RO 
^RAl 2790 
TRAl 2800 
T P A I 2 8 1 0 
TRA I 2820 
TRAl 2830 
TRAl 2840 
TRAl 2850 
TRAl 2860 
TRAl 2870 
TRA! 2880 
TP Al 2-(90 
TRAl 2900 
TRAl 2910 
TRAl ^92 3 
TRA I 2930 
YR Al 2940 
TRAl 2950 
TRAl 2960 
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SUeROUTlSE OYNA ~ TRA12970 

/TRA/P I • mija,ORAG tR PR • L PR t WPR, IJ , W8, 0 FLU 1 0 EL T , T I * TETA # I RAM, TRAl 2980 

, IL.ILO. Y.YO, TCB,Xl,Yl,X2.Y2.X3,Y3,^PRH,R,TORfl, I TR A , TF^ 100 , X C3 , YOO TRAl 2990 
CQ«MrN/CEO/ARFF,AO, AUA2 . A 30 , A 3. A5 • A 6 , L3 • AO AR , A5 AR , ILOO,IPR TRAl 3000 

CCM-ON/ A IR/GA ,GA l,GA2,G.U3,GA3 .PHOA, TO , P 0 , UO , MO , G RA V TRA13010 

CC'-MOs/t^LP/W A, WF ,wr , r.PDOT ,CST AH,PTL, TTL tr a I 3020 

rOyMf;M/PES/CF , THRUST , I SP ,SFC TRAl 30’0 

CO-mp'^/cCN/RTIc, MIC, A lC, ALFA TRAl 3040 

real •^0 , IPR, L3 .mua, I SP ,miC TraibqSO 

X3=(THRUST-0RAG)*CnSlTFTA)*nFI ?/vj pft^.2*X2-Xl TRAl 3060 

Y3={-GRAV+>((THRUST-0FAG)^SINITETA) )/ WPP) ♦0ELT**2*-2*Y2-Yl TR A I JO 70 

CCC DATA ‘jcxr LOCP: TRAl 3080 

TI = TI *>r>ELT TRAl 3090 

TETA*AT AN ( (Y3- Y2 I/I X3-X2 )) TRAl 3100 

UX=( Xi- X2I /OfLT TRAl 3110 

UY^(Y3-Y2 l/OELT TOaI 31 ?0 

U = SCPT { LiX**2 ♦UY^^Z) TRA13130 

Xl-^2 TRAl-^140 

Y l=Y2 TRAl 3150 

X2=X3 TRAIJI60 

Y2=Y3 TRAW3170 

RrTtNN TRAl 3130 

ENO TRAl 5190 

TRAn?00 

SUesnuT I NF OATTA TRAl 5210 

CC'-MrN/GEC/AOEF, AO, Al,A2,A30,A3,A5, A6,L3,A0AR, A5AR,IL00, IPR TRAl 3^20 

t . A IAO,A2A0 ,A<jA5, A3AR TRM3230 

CCVNC N/A IR/GA,GAI.GA2,GA12,GA3,RH0A, T0,P0,U0,M0,GRA V TRAl 3240 

COMMCN/FLfL/PHCF ,FTAT, A, N TRAl 3250 

CCV’ifM/LCSS/ PI 01 ,PI 02.PI 81, P IB 2, PIN TRAl 3260 

CCMMCN/0LP/W A, WF ,hT , F, R00T,CST AR, PT4,TT4 TR A I 32 70 

CCMM^^'/G^r/3F ,GF,CF I ,GF2 ,GF12,GF3 TRAl 3280 

CO'^MnN /NCZ/M6, PT6.P6 TRAl 3290 

COVMrN/PES/CF ,TI-PUST ,I SP.SFC TR Al 3300 

rcvYPN/rnr/PTS ,P5 trm3310 

COMMpM/ T M /o TO , TTO TRAl 5320 

CGYN'r>i/cCN/P TIC.MIC, AlC, Al FA • TRAl 3330 

CC'^MCN/ ThPO/ Ml I , PTl I ,M2l ,PT2 l> MSI ,PTS1 TRAl 3^40 

CC'^MnN/NS/’^l 2, PTI2,M22,PT22,MS2,PTS2 * TRAl 3350 

CC''M(]\/niF/V2, PT2.P2 TRAl 3360 

CC'-Mf'N/Cht/M 3N ,M3 I , PT3 ,PT3M ,M4 ,AS TR A 1 33 70 

CCMMnN/TPA/t> I ,.MM4,0RAG,RPC ,LPC ,WPo , U , WB, 06LU, DEL T , T T , TETA , IR AM , TRAl 3330 

,IL.IL0,Y,Y0, Tnp,xi,Yl,X2.Y2,X3,Y3 PR0,R , TORB, I TR A , T ET AOO , X 08 , YOB TRAl 3390 

CC'"MON/ORG/TcTP,Cr%,CFT,Cnww,COWF . A P R , SPR , SWW, 0 , XMO TRAl 3L00 

F.FAL “0 ,v IC, ^l I , Y12 ,M2 ,m3n, «3I ,M4,M5 ,M6 TRAl 3410 

real L3,A, ISP,M2Um22,MSI,'-S2,HUA,LPR tpai3^,?o 

CCC GECMFTOICAL OATA: TRAl 3430 

CCC -INCHES OF SU • I NC H^S ( OR I G I NA L J TRAl 3440 

CCC -kEP:4PC'/.N P.17 CR CSO P.13 TR.AI also 

CCC LOSSES TRAl 3460 

CCC PlCl*B.L ..PI 02 = SUPS0NIC 0 1 F . R E COV FR Y , TRAl 5470 

CCC PINsNO/ZLE LOSSES TRAI J4>U) 

CCC IP-.FCP OETALFC PRINTIMGS:*0 CLEAN RAM+TRAj;*l NO PRIN; TRAl 3490 

CCC *2 ALSO FRIN; = 3 ALSO LOOP UN m aCH . ; *- I TRAJ,ONLY. TRal35O0 

CCC ITOA.FCk TRAJECTOPY S lOOP CM AREAS TOAl 5510 

CCC ^ILL «ORK CN CNE P01NT; = -L WILL LOOP L PRINT ONLY SUMMARY TRAM520 

CCC ino=2 miJRN; = 3 n 022 ;*4 ChOKE; = 5 hEAT;*6 INLET(61 min, TRAI 3530 

CCC 62 throat, g3 0IFFUSI; = 7 CHPCK;=B C0RVAL181 EXPAN, 82 IAS.GT.15); TRA1 3540 

CCC «9 RFSUL ; = 20 TRAJI22 DRAGril. TRAIJ550 

CCC X'-Osmin. MACH MJMHER allowed. TPAl 3560 

CCC ItAMxC FAMJET:=1 PkOJECTlLE. TRA 1 3570 

ARFF = 19 . 3*.0 254»*2 TPAl 3580 

READ IS,*) AOAR,45AR,A lAO, A2A0 ,A3AR T^A13590 

READ (S,*l ALF A, TETAOD.TFTP TRA1 3600 

READ (5,*) I PR , I TRA , XMO, I PAM TRAl‘'6lO 

READ IS,*) P lOl,Pin2,P IN TRA 1 3620 

AC^AOAR * 4»FF TRA 1 3630 

A I =A0* A I AO TRAl .5640 

A2=A0*A2 AO TRAl 3650 

A3C=AP ^E*A3AR TRAl 3660 

13^A30 TRAl 3670 

A5=A51R*AREF TRAl 3630 
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FILH: TRAjei FORTRAN A NAVAL. POSTGRAOUATF SCHOOL 



A6*AftPF 

RPn*?. 5»C,025A 
IPF360.97»0. 02 54 
WPF=104. 7*0. 4535S 
P!=3. 141 5926 
L3*23.*.C254 
CCC Aio PLOW 
GA=l .4 

CCC ftHCA !N KG/H3,T0 IN OEG. KELVIN, PO IN Af M. ARE GIVEN FROM TRAJ 
CCC FlJfl^HTPe 
ETATrO ,9 

CCC DATA FIR ROO T*A*NAA3*»N, IN/SEC 
AsC.06 
N = C.6 

RHC^*0.0?5l 

CCC RHCF IN L3/IN3,=972 KG/M3 

flhrFrt^HnF*.4 5 3 59/(2. S4/100)**3 
CCC PLI^'^HT CCNOITlCNS: 

CCC UC !N FT/SEC;lFT.s30*46CM 
LIV.-25J0. 

UO^UC*. 3 C 4 ^ 

G? iV = 9. dC7 

CCC GRAV in ^-/SEC2 . o*PERFECT GAS CQNSTANTIFOR AIR). 

R» 29.314 

RET»M 3 N 

end 



TRAI 3690 
^PAl 3700 
TRA13710 
TPAl 1720 
TRM 3730 
TRAI 3740 
^PAl 3750 
^RM 3760 

tr A i 3 no 

TRA I 3780 
T3 A I 3790 
TPAl T3D0 
TPAL 3810 
TPAl il20 
TPAl 3630 
TRAI 3840 
TP.U 3850 
TRAI 3860 
TR A I 3H 70 
TPAl 1880 
TPAJ i 3<^0 
TRAI 3800 
TPAl 3910 
TRA I 3070 
TRAI 3930 
TPAl 3940 
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APPENDIX E: COMPUTER PROGRAM LIST OF SYMBOLS 



PROGRAM EQUATIONS UNITS MEANING 

SYMBOL SYMBOL 

GEOMETRICAL SYMBOLS 



AREF 


V 


m ^ 


AO 




m ^ 


AIC 




m2 


AJl* 




m2 


AO 2* 


''J2 


m ^ 


A30 


*30 


m ^ 


AIAJ 


A,/Aj 


- 


L3 


4 


m 


LPR 


S 


m 


RPR 


% 


m 


ALFA 


a 


deg . 



Reference area 

Area at station j 

Area behind a conical shock wave 

Area ahead of a normal shock wave 

Area behind a normal shock wave 

Initial area at station 3 

Area ratios 

Length of combustion chamber 
Length of projectile 
Radius of projectile 
Inlet cone half angle 



*When: 0=1, 2. S, the shock wave is at station 1 2 real case, respectively. 
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PROGRAM EQUATIONS UNITS MEANING 

SYMBOL SYMBOL 

ATMOSPHERIC SYMBOLS 



TO 


^0 


°K 


PO 


Pq 


kg/m^ 


RHOA 


Po 


kg/m^ 


MUA 


’PQ 


N. sec/m' 


UO 


*^0 


m/sec 


MO 


f'o 


- 


GRAV 


g 


2 

m/ sec 


GA 


'^a 


- 


GAl 


- 


- 


GA2 


- 


- 


GAl 2 


- 


- 


GA3 


- 


- 



Static temperature at altitude 
Static pressure at altitude 
Air density at altitude y^ 

Air viscosity at altitude y^ 
Projectile muzzle velocity 
Projectile initial mach number 
Gravity (9.807) 

Air heat capacities ratio (Cp/c,^) 

( Ya +1 )/2 

(y,-1)/2 

(YjY1)/[2(Yj-1)] 

Yy(Y3-l) 
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PROGRAM 

SYMBOL 



EQUATIONS 

SYMBOL 



UNITS 



MEANING 



COMBUSTION CHAMBER'S SYMBOLS 



RHOF 


Pf 


kg/m^ 


Fuel density 




ETAT 


^T 


- 


Burning efficiency 


A,N 


A,N 


- 


Burning rate parameters 


WA 




kg/sec 


Air mass flow 




WF 




kg/ sec 


Fuel mass flow 




WT 


^T 


kg/sec 


Total mass flow 




F 


F 


- 


Wf/*a 




ROOT 


r 


m/sec 


Burning rate 




CSTAR 


C* 


m/sec 


/g*R£*Tj- / r(when:r = 


RF 




m/°K 


Hot gas constant 


r_ r(j7 

_ MW(k 


GF 




- 


Hot gas heat capacities 


GFl 


- 


- 


(Yf+1 )/2 




GF2 


- 


- 






GFl 2 


- 


- 


(Yf+l)/[2(Yf-l)] 




GF3 


- 


- 


Yf/(Yf-1 ) 
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MEANING 



PROGRAM EQUATIONS UNITS 

SYMBOL SYMBOL 

THERMODYNAMIC SYMBOLS 



TJ 








TTJ 




"tj 




PJ 




p . 
J 


kg/m‘ 


PTJ 




"ij 


kg/m' 


MJ 




"j 


V 


PTJl 


, PTJ 2 


*^Tjl’PTj2 < 


1 


MJl , 


MJ2 


1 


1 


T3M 




T3M 


kg/m 


M3N, 


M3 1 


^3N’^3I 


- 



Static temperature at station j 
Total temperature at station j 
Static pressure at station j 
Total pressure at station j 
Mach number at station j 

As above with AJl , AJ2 

Maximum available 
M^ calculated from nozzle and inlet 

direction, respectively. 
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PROGRAM 



UNITS 



MEANING 



SYMBOL 

PIDl 

PID2 

PIN 



EQUATIONS 

SYMBOL 



^d' 



^NS 

% 



LOSSES SYMBOLS 

Boundary layer losses 
Subsonic diffuser recovery 
Nozzle losses 
Conical wave losses 
Normal shock losses 
Expansion losses 
Heat losses 
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PROGRAM 

SYMBOL 



CF 

Thrust 



ISP 



EQUATIONS UNITS MEANING 
SYMBOL 




RAMJET PERFORMANCE SYMBOLS 

Thrust coefficient 

N Thrust 

(or kg) 

N/kg.sec Fuel specific impulse 
(or sec) 



SFC SFC 



kg/hour/N Specific fuel consumption 
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PROGRAM 

SYMBOL 


EQUATIONS 

SYMBOL 


UNITS 


MEANING 






TRAJECTORY SYMBOLS 


Drag 


D 


N 


Drag 


WPR 


''p 


kg 


Mass of projectile 


WB 


"b 


kg 


Mass of booster 


DELU 


AU 


m/sec 


Change in initial velocity due to 
booster 


U 


U 


m/sec 


Uq + AU 


DELT 


AT 


sec 


Change in time 


TI 


t 


sec 


Time 


TOB 


^OB 


sec 


Time of burning 


TETA 


9 


deg. 


Gun elevation angle 


TETP 




deg. 


Projectile second cowl angle 






DRAG 


COEFFICIENT 


CDN 


^DN 


- 


Nose drag coefficient 


CDWW 


''DWW 


- 


Wing wave drag coefficient 


CDWF 


''DWF 


- 


Wing friction drag coefficient 


CDS 


''DS 


- 


Skin drag coefficient (1 aminar/turbul ent) 


CDB 


^DB 


- 


Base drag coefficient 
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PROGRAM 



EQUATIONS 

SYMBOL 



UNITS 



MEANING 



SYMBOL 



PI TT 

IPR 



ITRA 



I LOO 



IL, ILO 
I RAM 



XMO 



^MO 



MATHEMATICAL SYMBOLS 
3.I4I59 

Printing parameter: > 0 combustion 

results together wi th'traj ec to ry 
(on different files) : 

= 0 working results only; 

= 1 also reasons for not running 

= 2 also full reasons for not running 

= 3 also loop on mach number (CALCM) 

= -1 trajectory prints only 

Loop paramerer: 

= +1 single value for Aq/A^, A^/A^ 

= -1 loop on Aq/A^, ^ 5 /^r’ 

print summary, only. 

Check parameter: 

< 1 regular run 
> 1 doesn't have a solution 
Trajectory printing parameter 
(prints every ILO point). 

Ramjet parameter: 

= 0 ramjet in operation 
= 1 projectile without propulsion 
Stopping mach number 
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PROGRAM 


EQUATIONS 


MEANING 


SYMBOL 


SYMBOL 


SUBROUTINES 


INIT 


- 


Computes initial conditions 


BURN 


- 


Computes combustion chamber's performance 


TERMO 


- 


Thermodynamic comand subroutine 


INTER 


- 


Computes, by interpulation, thermodynaic 
conditions, or cowl drag coefficients. 


NOZZ 


- 


Computes nozzle performance 


CALCM 


- 


Computes mach number- indirectly 


RESUL 


- 


Computes ramjet performance 


CHOKE 


- 


Checks if nozzle is choked 


INLET 


- 


Command subroutine for inlet 


CONE 


- 


Computes conical wave loss 


THROAT 


- 


Computes boundary layer loss 


NSR 


- 


Computes normal shock loss 


DIFFUS 


- 


Computes subsonic diffuser performance 


MIN 




Command subroutine to compute situation 
when normal shock is at point 2 


-CORVAL 


- 


Conmand subroutine to compute situation when 
normal shock is at the correct place 


EXPAN 


- 


Computes losses due to expansion into the 
combustion chamber 


HEAT 


- 


Computes heat losses at combustion chamber 


CHECK 


- 


Check for pressure capability 


TRAJ 


- 


Command subroutine to compute trajectory 


ATM 


- 


Computes atmospheric conditions as a function 
of al titude (y) 
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PROGRAM 

SYMBOL 


EQUATIONS 

SYMBOL 


MEANING 


BOOS 


- 


Computes booster performance 


DRAGG 


- 


Computes drag 


CALCD 


- 


Computes skin drag coefficient, indirectly 


DRAGN 


- 


Cowl drag command subroutine 


DRAGWW 




Computes wing/fin v;ave drag coefficient 
(command subroutine) 


PRANT 


— 


Computes Prandtl -Meyer angle (v) 
from a given mach number 


APRANT 


- 


Computes mach number from a given Prandtl - 
Meyer angle (v) 


PRES 


- 


Pressure ratios formula 


DRAGPR 


- 


Computes drag coefficients of a projectile 
without combustion 


DYNA 


- 


Computes the dynamics of the projectile 


OAHA 


- 


Initial data 


PRILOO 


- 


Prints ramjet performance 


PRIN 


- 


Prints detailed values when does not fine 
sol ution 
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APPENDIX F: COMPUTER PROGRAM USERS GUIDE 



FI . Input Data 

AOAR, A5AR, AT AO, A2AO, A3AR 
ALFA, TETA, TETP 
IPR, ITRA, XMO, IRAM 
PIDl, PID2, PIN 
Options 

ALFA = Inlet cone half angle 
TETA = Gun elevation angle 
TETP = Second cowl angle 
IRP = 0 clean print of RAM + TRAJ 
=1, 2, 3 more details on RAM 
= -1 TRAJ only 

ITRA = 1 works on one set of data 
= -1 loop on A0£A5 
IRAM = 0 Ramjet 

= 1 Projectile without propulsion 

F2. Execution Commands (For use with IBM 370) 

F2.1 Opening Commands 

LAXXXXP 

Password 

GLOBAL TXTLIB F0RTM0D2 M0D2EEH 
F2.2 Compilation 
FORTGI TRAJET 
F2.3 Run on Terminal 

FILEDEF 02 DISK TRJ D(RECFM F BLOCK 132 PERM 
FILEDEF 03 DISK DRG D(RECFM F BLOCK 132 PERM 
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FILEDEF 05 DISK INP D 



FILEDEF 06 DISK CMB D 



EXEC 


RUN 


TRAJET 






XEDIT 


CMB 


D 






XEDIT 


TRJ 


D) 






XEDIT 


DRG 


0) 






Note: 


CMB D will 


contain 


the combustion process results 




TRJ 


D will 


contain 


the trajectory part. 




DRG 


D will 


contain 


drag coefficients. 


F2.4 


Hard Copy 






FILE 










PRINT 


CMB 


D 






PRINT 


TRJ 


D 






PRINT 


DRG 


D 
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APPENDIX G: 



Gl. PROGRAf-1 AERO^*^ : LISTING 



(*) The original program was developed by T. H. Gawain [9]. The modification 
listed here was prepared for this report to calculate the cowl drag coefficient. 



or^o ooo o o non 



FILE: ACnn 



FORTRAN A 



NAVAL POSTGRAOUMF SCHOOL 



SPECIFICATION STATCh^NTS 

implicit RFAL*4< A-H.O-ZJ . -mi *' 

0( mcms ion XC 7021 702 ) ,RPt 7 0? ), XM 7 02)* fo< 702) .QP( 70?) . 
♦GI 70? ) I 0 AOI 7 02 )i O''?! ^0 2) ,0NA (702) tOM 4( ?0** ) ,KAP ( 0 ) 

DATA <rr/*oc * / tK\t / • 'M*/,KCG/*Crt»/,<rn/tfn«/,<Tn/»Tn*/, 
4.<of / *pn • / ,<0 T/ *0T • / , <ON/ / , KpT/» 0 T* / , <M/1. 1000/ , XL/' 7, 



Nl /SOO/ ,Nt\/20 3/.\S/? )0/ 

4-kap/*x cc * » ' '■^r ri • , • n* TP* 
*'.QL/-0,073i>0/, CL2/-0 .) 7/6/.MC/ 

Sei'^CT OPTIONS 



TAB/ 20/ tXNO^P/ ‘PM* / , <T M L / T • 
, • ( ’ 0 = • , * N£ r.T • , • T ( V = • , • ) 
6/ ?/bl /, IMPPP / 0/ 



1 001 WRI T=( 6 ,1021 ) 

1021 POPMAT (•'*,• '■NTPo rnNTHDL COOP AS FOLMWS:*/ 

♦ T6,*^0 TgTMN Ex AN *.T lOM* . T4 I , • CNT c? np • / 

♦ T6,*'’"* <‘N'Tr^ P*.'' 'MCI ji p<; I ,74, ^ ^ # rvjTcp f^*/ 

♦■T6,*^n rCMP)jT<“ O^^-^PTS Tf Al ARO AYS • , T , • PNTcq C 

♦■T6,'Tn CAirULATP *tj T P I T • ^ » CNT r^*/ 

'►T6,*Tr> '^•iTpfjr • , T VI , f NfpP TP*/ 

♦•T6,*Tn plpT niJTPtiT* I , *'Nr Pt • / 

♦ Tb.’Tn OUT PRCGO^M* ,T^l , • CNTFR 0T'/| 

PFAOib , I C<il, FPOz IGAI ) <AR 

lO^l 

IFlK-P.rC.KTF) 



G*/ 



1061 

lOBl 



2001 

2021 



2.001 

2 OR ! 

2 I 

3 2BI 
^3Rl 
^601 
3721 



coonp. cgpEAT I NPUT* ) 



( F (KAP , rC-Kl Ml 
t F(K*'5 ,CC.<CG) 

I F(K*.o ,rc,xco ) 

I F (K ‘0 . rc .<TfU 
I F (K*R . PC.<Pn) 

TFit^^P , = C.K0T) 

WR ITP (6 , lOSl ) 

FOP'^AT (I • , I if4p»jT 

GO TO IJCl 
EXPLANA T ION 
WRITE ( 6 ,202 1 ) 

FORMAT!* • , • SUPPRSONir PLOW ABOUT A SLIGHTLY »*TrHPO, SLEMOER,*/ 
noiMcf) PCOY nc u fvni 'JT ION. THIS PROGRAM JTILUpS A*/ 
•'OniM'^O LC*ARTZ=n ANALYSIS OEV-LOP^^O BY T ,H , G Au A I N , * / 

ANO «'^0(FIF^ FQP CnwL ORAG calculation op a GlI'OCTt/ 
•'PojfctilF rv O.amichai.'/ 

riASSICAL L!Nr\RIZFO th=ORY SATTSFIFS TH" BUMNOAPY*/ 
rONCITIONS fNLY APPROXTmatR Y, AT THF AXIS: TH*" ropccNT*/ 

mcthOO SATIPFIES THF** VFPY Arr.URATELY AT THE BODY SJPPAf.c*/ 
^TSELP.*/ 

TFP BODY IS SCALFO SO THAT THF LAPGFST CPOSS-SrCTION IS*/ 
A CIPCLr OF UNIT RADIUS* AERnOYNAMfC rnFPF(CI=MTP APP 
ON FOONTAL 



2041 form 



TRC CONFIG JPAT Tr.>i CONSISTS OF A COM KM 
*K\r^-0 BY S CN OM ON, )P S T ' V^l Y 

OotiOnAI «I r-P E'"" (ON OF UMiFOtM IN(T ^AOI 



PR '^GW^L “OSP, 
, IN*/ 
jp , ^*‘ln AN 



AFPOOOlO 
AF»0 JOPO 
APR 3 3030 
ApROO0<»0 
Ar->ioO‘=0 

A PRO 006 J 

»(CR0 >07/) 

AP' 0 00^0 
AP' 3 00 00 
50/, AFOOjjOO 

, A^RO-^nO 

• / APR ? ) no 

Ar''0 >1 00 
AFP ^0(40 
APRO.jlSO 
AP'»0 U60 
AC'' 0 0 I 70 
4FR) >IRJ 

ACPfi '\O0 
» CO f> 7*>0 
AC-?') ) 

ACRO )220 
4CRO 1730 
\P '0 0240 
ACRO'JPPO 
APOO )2(y0 
AFRO C270 
\cn )R.BO 

\P'-» A?O0 

A = P OO300 
AP30 ) A1 0 

Ap^o n?n 
“ ACP)0330 

• coi >140 

APR 1 33S0 

AP^ 0 03 
APRO )1RQ 
APROO3N0 
AF*’0 04 00 
AC^O 04! 0 
AFRO )420 
A = RO 04^0 
AAf'O )44 1 
AC-0 >4F.l 

ACO 1 144,0 
AFT,! )47o 
AF' > 04B0 
AFOOOAOO 
A poo osoo 

ACO 10510 
BAS CO* /AFRO )5R0 
\F?1 )S30 



AFRO A543 
ACC A >560 
OOTI PN-* / AC- ) i'560 



Al TML (N -F A TRD-ifAf^r) CTNC IR '^.I'^C, F t I / j qg /q 

OP FT, P r SP' ft ( '/' LY. total L = NSTH XL 'S S t JO 0 I >/ * or o r^ro ML*/ *foo->sR') 
INTCrv.^n OF flNlFORM W I O'" M H= < L / *'»L • NOSC |CMr,TH «- j q- » / 

ccf T(nN I PNf;>w >OF PXORFSSCO 'S NA*H ANU MR*--', 3'*S®‘^CT-'/ \co0)600 

TVCLY, WHFRF M4 4NO MB AP F INTEGERS. <*ARAMCrcp Cl orMQtr-.) AC'O^blO 
WRITP(6 ,?C4l ) Aror)0b20 

T{* *,*thE F'CRCiS^ 'N RAPIIIS ovcp TMC T\IL StCTTOM.*/ AEROJbRO 

CAICIII AT I>NS \PC •-•AOE Fri2 cy/p-y STATION J OHTnijr^UT*/ AFR01640 

IS TABtJlATFn TNLY p 0® J=2*-K*NTaP ^HCPF S O . I , 2 , ? , . . . • / AFO^RbSO 

math NUM»« P mjj^;T lit W'THlM THC OANGC l.<'Xw<*/ AFRQ >660 

SOF T ( 1 , *P p**ak*»-2 ) WHE'^F ‘■Pmax rjP^nTCS T»|C MixiMiju nOY*/ AFRvA >b/0 

c 1 npc . * / K 3 1 >(,.in 

INPUT PC RAMI- 'c ARE <M, ^L, '3L,Nlt NX, |H0 , r.TA(B,CfJ 0« ’ / 

OH. Cl OB or, A*''> npTTrrni M RHAMlI'^'f R ' C IOENTT-'/ *r-' > IN'*!' 

CjCB !)p TO ,< fh^RACKRS 'M length,*/ A=oO)MJ 

A CONTpni i'll!'!! IS I fJI ATC ' AT ^=0., Al I H£ Mlf>!>'>IMF OF*/ APR'>n/20 
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FILE: AERO 



fortran a 



NAVAL OOSTORAOUATE SCHOOL 



*^ACH OF THC NL ONIFOR^ INTERVALS. ANO AT X*XL. INDEX J»l,V 
*2.3.*., (ML^-Z) OESIGNATES THESE POINTS, AT POTNT J THE AXIAL* 
*r.OnooiN ATF, the PAniM CrOROrUTF AND the body slope are*/ 

•CXPRESS=0 A$ X(JI. R(JJ ANH RP(J), PP SPFCT I VFL Y. THPSF'/ 
•At^oAYS ARC CALCULATED AFTER THE INPUT PARAMETERS HAVE*/ 

'BEEN EN Tcpfr).' / 

. THP ‘j'lcvrL and axial FCPC? COEFFICIENTS AP.F OF THE*/ 

*ROt losing form. wHERF S ^NO r OENOTF SINE AMO CQS^NC OF A,*/ 
•Thf A'iOlE AT'^ACK. THUS CN=CNA*C*S AND C A =C AQ ♦C A2 * S ** ? , * ) 
WRITF( 6. 2061 ) 

2061 FOP.’’'T(‘ This PROGRAM fOMpijIFS Th^ NUMcrICAL VALUES OF*/ 

■ ■ ’COEFFICIENTS CAC AND CA2 AS «ELL AS THE COORDINATE XAC* 

•rr ThF AFRf'-CYNAMI C C^NTfJ- , PA^AMr^E = S C.NA , CAO AND CA2 ARF*/ 
•F0U^!) ^Y INTCjO pr f DRP cSP ON 0 I NG DISTRIBUTIONS ONA, QAO*/ 

• iNO OA7 WITH £?FSP=CT to X OVER IMF. LENGTH OF THE ArnY.’/ 

*ThFSE T I ST^ I3UT ions r \n cE ODTAINED IN TABULAR AND ^.paPH-*/ 
•TCAL FT = m. th= CAN 3G CCNTROILED BY FOLLOWING THE’/ 

•PRO^-PTING I NS'’-RiJC TI DNS wH ICH APPEAR AT THE TERMINAL. TH= • / 
ipRn(TOv< STEPS "CrjTFR INPUT o AR A MF TF PS'* , '*C3 *^PUT ';F0M = T R I C ' / 
•A-RAyS" ,‘*CALCHl ATF OUTPUT" A.ND "'^LOT OUTPUT" MUST BF pcr-i/ 

• cnFvp*^ IM TM=- n.cfijc I ISTFJ SINCE EACH ST^P GENFPaTES OATA’/ 

• »;pcpcn pno THE '!=<*■ ST^p.*) 

«FITr(c'.20o5) 

2065 FQrV*AT(* *.* CHANOFS: ADDING SFCONO CONICAL TAIL.'/ 

*N‘A2.NC2 ARF fNTEC.== S. SI MI L lAR to ThF AQn^ = ;DL2 DFNOTES thE*/ 
•srrCN’D DECRFASF in radius 0VF-» THF tail SFCTinN;IN ORDER '/ 
*"0 PFCEfVP INCt^‘'ASc IN RADIJS.0 LiDL2 SHOJID RR NFGFTIVF,'/ 
•S“CCNO TAIL NAS NUT ARAN&FD FOP OUl VAL TAIL PROFILE.*/ 

*ThF NOSF pa-”*' (NA) can Kf INCLUDFD IN CAO C AL C'Jl AT I ON * / 
*(!NCS--l) 3R IT can mf ysFO ONLY TO ARRANGE THE FLOW/ 

•CONCI TI uNS( IN0SE = 0) * ) 

GO TC IOC I 

c 

C ENTFC INPUT '’APAHFTERS 

C 

2Cai WR ITE( 6 . 2101 ) XM. XL .DL »NL..NA ,NB,iNTA3 .KNOSE.KTA IL,KAPt6I.XAP(7) 
•tNC.VR2.CL2. INOSF 



!101 FOpmat ( 



oREST^T INPUT parameters ARE:*/ 



♦ • 


* . • I TCM 


1 


XM=* . 


F T.*^/ 


• 


* . * I TF V 


2 


XL=* . 


F7.3/ 


♦ * 


*.*ITCM 


3 


DL 5^* . 


Fo • 4 / 


♦ • 


' . * I TC*- 


4 


NL= * . 


I 3/ 


• 


• . . Itpm 


5 


N As* , 


13/ 


♦ • 


* . . I TCy 


6 


N K := « , 


13/ 


■f • 


* I * I T F y 


7 


. NTAB = 


*. 13/ 


♦ * 


* . * I TE M 


6 


N'JSt 


Cnncr* ,a2/ 


* 


' . * I T F M 


9 


T A !L 


CnoR = * . A2/ 


*■ • 


* , * I TP*. 


10 


I DFVT IF IC ATIDN NUMBER ■ 


♦ • 


* , • I TC M 


ll 


NC = * . 


13/ 


♦ • 


*. * Item 


12 


r)«?= * 


. I 3/ 


♦ * 


1. . irrv 


13 


D1 2s • 


,F6,4/ 


♦ • 


*. • I T = y 


14 


IMOSF 


^* , I 3// 



i2A4/ 



TO CHANGF any IT«"m ^NTFR ITEM NUMBER IN 13 FORMAT.*/ 
TO FXIT ENT^P -01*) 

2121 REA0(5 .2 Ul, ERP = 2Itol ) <AR 
214.1 coP'*aT( 13) 

I E (<AR . rc.oo I ) or TO 2201 

I FIKAR .FC.0U2) GO TO 2 26 I 
I F(kA‘ . rc .03 3 ) GO TO 2A2 1 
I F (KAP . rc. JO^) GO 2 3BL 
IFIKAP .EC.0J5) GU TO 2^^/, I 
IF (KA'’ . FC.DD6 ) GO tq 24B I 
IFIKiS ,=C.U07) CO rr] 252 1 
IP (KAR . J^C.OO rf ) GO tq 256 1 
IFIKAD .FC.00<M go to 2621 
I =(<AP .EC.JIO) GO ’’0 26ttl 
IF(KAP .tC.Ol I ) GO TO 2 60 3 
I . EC. Jl 2) GM ’•0 Z6H5 

! F iXAR . rC.Dl 3 ) GO TO 266 T 
I F{<A= . FC. )l H) G3 TM 2689 
I F(K AR . EC.-3 I ) GO ro lOUl 



AFR00730 
/AER00740 
AEP007S0 
AER00760 
A ED 00770 
AFP00780 
AFR00790 
AEPO JflOO 
APDDOblO 
AFR0U820 
AEPOORDO 
/AFRQ0840 
AER0C850 
AFD00860 
AEPD0670 
AERDO08O 
AFP00890 

aeddqdoo 

AERO 3910 
APPDQ970 
AFK 0G930 
AERO 0940 
AFR0G950 
AER 30960 
AFR00970 
AEE00980 
AFP0G990 
AFRO 1000 
AEPO lOlO 
AFRO 1020 
AFPO 1030 
AERO 1040 
A = RO 1050 
AERO 1060 
AERO 1070 
afro lOdO 
AFRO 1090 
AERO 1 100 
AFRO 1 1 10 
AFCOU20 

AERO mo 

AFRO I 140 
AERO I 150 
AERO 1160 
AERO 1170 
\FPO 1100 
AFRO 1190 
AFDQ 1200 
AFR01210 
AFRO 1220 
AER01230 
AFPO 12 VO 
AERO 1250 
AFR 01260 
AHPO 1270 
AFO01280 
APRO 1290 
AFRO 1300 
AFR01310 
AeR0l320 
A=RO 1330 
APR3 1340 
AERO 1350 
APD01360 
AER01370 
AERO 13R0 
APP01390 
AER ) 1400 
AFRO 1410 
AFRO 1420 
ACRO 1430 
AER01440 
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F!L6t AERO 



FORTRAN A 



NAVAt P0STGRA0UAT6 SCH03C 



2161 

2181 



INPUT ERROR MESSAGE ! 

mRITF(6.2181 ) 

FGP‘'AT(« •,» INPUT ERROR. REENTER ITFM ^UMBER.*] 
GO T'O 2121 

ENTER MACH NUMBER XM 



2201 w® ITF|6,2221 ) 

2221 FOSMATI* »,• enter aM IN OFCIMAL FORMAT. 

RFAT( 5, 22‘»l. rRR=2l6l ) XM 
22^1 FCFMAT(F7.^I 
GO TO 2051 

FNT®R length XL 

2261 -PIT® ( 6 , 2231 ) ' 

2281 FQR-AT(' •.'FNTFR XL IN DECIMAL FORMAT, 
OFaT(5.230l,cPR = 2l6l ) XL 
2301 FCP^AT lFt.3J 
GO TO 20 El 



* ) 



ENTER tail taper DL 
, 2341 i 

2341 FCR''M(* '.'^^NTPR 01, IF APPU CABLE, IN DECIMAL FORMAT.*) 
RF aH( 5, 23ftl» F0R52161 ) OL 
2361 FCPMAT(Ft.A) 

GC TO 2051 



EME® BQCV LFNGTH INTEGER NL 



2381 woTTFi o, 2401 ) 

2401 F0RMAT(» '.'CNTFR ml in 13 FORMAT.*) 
PEA0l5,2-42l,rPR*2l61) NL 
format (13) 

GO TO 2u81 

ENTER NOSE LENGTH INTEGER NA 



2441 WRIT® ( 6, 2461 ) 

2461 ®PPmaT{* '.'cnTFQ NA IN 13 FORMAT.*) . 

REA015, 242lf FRR*2l6l ) NA ^ 

GO TO 2061 

ENTER MICSECTICN LENGTH INTEGER N8 
2481 WRI T=( 6. 2501 ) 

2501 FOCmAT ( * *,*®NTER NB IN 13 FORMAT.*) 

RFAO(5,242lf FFP=2l6l ) NB 
GC TO 2061 

ENTRO tabular INTERVAL NTA8 
2521 wR ITP( 6 , 2541 ) 

2541 FORMAT! • '.’ENTFR MTA3 IN 13 FORMAT.*) 

PEAn(5* 2421 , FROs 2161 i NTAB 
GO TO 2061 

C FNTRR COCE for CONICAL OR OGIVAL NOSF 

C 

2581 FPPMAT? » ^ * ?* FN TFR nose code "C N" OR "ON**.*) 

R ®A0(5 , 2601. ER R3216 I ) <N0SE 
2601 F0PMAT(A2) 

GO TO 20 81 

C ENTEF CUCE for C^'NICAL or ogival TAIL 

C 

2641 F?PMAT*( • FNT®R TAIL CODE "C ’’•* OP ”0T" IF APPLICABLE.*) 
P£A0( 5,26bl, EPR»2l6l ) K.TAIL 



AER01450 
AEP01460 
AEPQ 1470 
AERO 1480 
AEPQ 1460 
APPO 1500 
AFPO 1510 
AFQO 1520 
AFQO 1530 
AEPQ1540 
AFPO 1550 
AEPO 1560 
AEPO 1570 
AFPO 1580 
afro 1500 
AEPO 1600 
AEPO 1610 
AEPO 1620 
AFPO 1630 
AEPO 1640 
AB^^O 1650 
AF^O 1660 

AEPO 1670 
AFQO 1680 
AFPO 1690 
APPQ I 700 
AFQO I 710 
AFROl 720 
AFP01730 
AFPO I 740 
AEROl 750 
AFQO I /60 

AFPO I 770 
AEPO 17R0 
AFQO I 790 
AEP01800 
AFQ 3 1810 
APQQ L820 
AEP01830 
\EP0 1340 
AFPO 1350 
AFPQ1860 
A®PO 1370 
AFQO IflRO 
AFQO 1800 
AFPO 1900 
AFQO 1910 
AFP01920 
A5P ) 1930 
AFPO 1940 
AF®0 1950 
AFQO I960 
AFRO 1970 
AEQO 1980 
AFRO 1990 
AFP02000 
AFR02010 
AER02020 
AEPQ2030 
AFPO 2040 
AFP 02050 
AFQO 2060 
A ®R 0 20 70 
AE»0 2080 
AFQO 3090 
AFPQ2100 
AFR021 10 
AFO02120 
AEP02130 
AEPO 3140 
A®Q02150 
AER02160 
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PILE: A£RO 



PORTRAN A ^4AVAL POSTGRADUATE SCHODU 



2661 



2601 

2701 



2721 



2683 

2703 



2685 

2705 



2741 



P0P»<AT(A2) ' 




AER02170 


GO TO 2081 




AER02180 

AFR02190 


PNTFR rOENT! F ICATICN NUMBER 




AEP02200 
AERO 2210 


WRITPI6.2701 1 




AER02220 


FnRMaT(» t,.cMXpq OPTIONAL 10 NUMBER 


IN A8 FORMAT . * 1 


A=R02?30 


R 5 , 2 721. PPR* 21611 i^A P( 6 1 .K AP( 7) 




AFR02240 


FCf '^AT ( 2 A4) 




AFR02250 


GO TO 2081 




AERQ 2260 
AE<?02270 


PNTE6 first TAIL LENGTH INTEGER NC 




- AERJ2280 
AER02290 


.^RITC(ft,27Q3 1 




AFR02300 


F'"P*-«ATI* »,»=NTER NG'' IN 13 PORMAT.M 




APO02310 


REAOl 5. 2 421. ERR* 21 61 INC 




AER02320 


GC TO 2081 




AFR02330 

AER02340 


= NTEP SECOND MIOSeCTTON LENGTH INT*^GER N82 


AF002350 

APO02360 


WRITFI o. 2705 1 




AER02370 


FORMAT { • ','PNTER NR2 IN 13 FORMAT.* 


1 


AFR02380 


REA0(5.242l.FRH*2l6l )Nfl2 . 




AFR02390 


GO TP 2081 




AFRO 2400 
AFR02410 


ENTER TAIL SECCNO TAPER 0L2 




A=R02420 

AEP02430 


WRIT = (6 .2707 1 




AF902440 


FCPMaTI • •,* PNTFR 0L2.IP APPLI CABLE, 


IN OECIMAt FORMAT.* 1 


AER02450 


REAP! 5, 2 361. FRPa2l6llOL2 




AFR02460 


GC TO 2081 




AER02470 

AFR024d0 


ENTER iriCSE OPTION 




AER02490 






AFR02500 


WRITFI 6 , 2709 1 




AFP02510 


FORMAT!* *.*PNTER INOSF IN 13 FORMAT 


*,/• • 


AER 32520 


IF I NCI UOFO IN CAO .*^R ITE 


1*./* * 


AFR02530 


.•IF ONLY TO ARRANGE PLOt»rf. 


WRITE OM 


AFO02S40 


R F AD! 5, 2 42 1. ERR* 2 161 H NOSE 




AER02550 


GO ’'0 2081 




AER02560 

AFC02570 


CALCULATE GEOMETRICAL ARRAYS 


\ 


AFRO250O 

AER02590 


SET UP INPUT AND CALCULATE X(J) 




AFR02600 

AFR02610 


IE(N82 .EC.Ol nC=NL-NA-N8 




AER02620 


NC2 = Nl -NA-NB-NC-NB2 




A=R02630 


DFLX = X L/FLfjAT( NL ) 




AER07640 


AL=3ELX*FLCAT( NA) 




AER02650 


SL^O^'LX^FLOATINB) 




AER 02660 


CL=PElX*RLnAT(NC) 




AER02670 


CL2=0FLX^FL3AT(NC2) 




AF902630 


JAL=NA>1 




AER02690 


JARaJAL ♦ 1 




AER02700 


J0L = NA 




AER 02710 


J8R-JEL> 1 




AERQ 2720 


JCL=NA^Ne*NC ♦! 




AFR02730 


JCR^JCL^l 




AER02740 


J»2L = N 4*NR^-NC4-NB2>1 




AER02750 


J-B2R = JB2L*l 




AER02 760 


X( l)=0. 




AER02 770 


p ( 11=0. 




AFR02780 


NLPI=NL ♦ I 




AFP02790 


NLP2=NL>2 




AER02830 


X (NL^2 1 - XL 




AERO 2810 


00 2761 J=2.NLPl 




ACR02H20 


OrLJ=FLnAT( Jl-1.5 




AERO 2830 


XU 1 =CELX *0ELJ 




AF502840 


COMINUF 




AFP 02.850 


IFIK'PSP .cO.KONl GO TO 3001 




AFR02860 


‘ 




APOU2870 


CalCLATE RU) and RPUl FOR CONICAL 


NOSE 


AER 02880 
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C 

2 801 



2 921 
2 9<*i 
2 861 



3001 



3024 

3026 

3021 
3 04 1 
3C61 



PP(l)sRPA 
00 2921 J*2. JAL 
RU » = XU) /AL 
PP( J ) -POA 

continue 

IFINA. lT.nl ) GC TO 3081 
R (r:L‘^2 ) * 1. 

KP (NLR ’ ) sPRA 
GC TO 3241 

calculate rij) anc rpu) for ogival nose 

00 302 I J=l, JAL 

(F((Ql^QA-(AL-XU))**2).GE.0.0)f*0 TO 3026 
WRITP{ 3. 1024 )J.PA,ALfXtJ),R(J) ,RP(J) 
FC-MAT 1 • • ,* 30 24 » ,2X. 13*2X,5(E10. 3.2X)L 
GG Tf; 17 77 

R< J ) s SORT (RA*PA-{ AL-X( J) )♦• 2)-RA> 1. 
3 P(J)=UL-X(J))/{Ra-I.^R{J)) 

CONTINUE 

IFINA.lT.nl) go to 308L 
RILLP2)=l. 

RP(ULP2 ) = .Q 
GC TO 3241 



CALCULATE P(J) ANO RPIJ) FOR CYLINDRICAL MIOSECTION 



3081 |F|Nn.-0.0) GO TO 3L21 
30 3101 J=JAR,J8L 
F U ) = I . 

KP( J ) =.0 
3101 CONTINUF 

IF ((NA «-Ne) .LT.NUGO TO 3121 
RINLP2 )= 1. 

RFI.NLR3 )*.0 
GO TP 32^1 

3121 IFUTAIL .cO.xnT) GO TO 3161 



CALCULATF RIJ) ANO RPIJ) FOR CONICAL TAIL 



QPTs-Ol/CL 
RXL»AL ♦•81 «*CL 
DC 3141 J*J8P.JCL 

RIJ) = l.^RPT*l X(JJ-RXL>CL) 

RP I J ) *FPT 
3141 CONTl^’Hf 

IF I (NAt>NE^*NC ). LT.NDGO TO 3123 

R|NLP2)*L.-">L 

RPIMP2)=RPT 

GO TO 3241 

CALCULATE PIJ) ANO RP(J) FOR SECOND HIOSECTIQN 

3123 IFIN32.EC.CI GO TO 3241 

00 3125 J = JCF , JB2L 
PU) = l.-OL 

®P( J ) 5.0 
3125 CONTINUE 

1 F ( (NA fNP^NC ♦•N82 ) .LT.NL) GO TO 3127 
R INLP2 ) 5 l.-OL 

»F(NLP2)s.O 
GO TO 3241 

CALCUL ATE .aIJ ) AND RPIJ) FOR SECOND CONICAL TAIL 

3127 RPT25-"*L2/CL 2 

DC 3129 J*J923 .NIP2 

P I J ) 5 1-OL ♦? PT2 *iXIJ )-XL>CL2 ) 

RP ( J ) aPPT2 



AER028*30 
AE»02900 
AFP02910 
AER02920 
AFR02930 
At^R02940 
AEPQ ^950 
AE® 02960 
AEP02970 
AER02990 
ACOQ2990 
AE®0 3000 
AFR03010 
AER O 3020 
AERO J030 
AERO 3040 
AERO 3050 
AFRO 3060 
AFRO J070 
AER0 3080 
AE20 3093 
AF»03100 
AERO 31 10 
AERJ 3120 
AFP03130 
4ER0 3140 
AEP^O 3150 
AER03160 
A*=oo 3170 
afro 3180 
AP®0 3190 
AERO 3200 
AER 03210 
^FP0 3220 
AFR0 3230 
AERO 3240 
AFR03250 
AFRO 3260 
AER03270 
AERO 3280 
AE® 03290 
AFP03300 
At^R033l0 
AFRCT3320 
AERO 3330 
AE® 0 3 340 
AERO 3350 
AER03360 
AE®0 3370 
AERQ3390 
AER03390 
AER03400 
AFRO 3410 
AFP03420 
AFO03430 
AER03440 
AFRO 3450 
AERO 3460 
AER03470 
AERO 3480 
afro 3490 
AEOQ3500 
AERO 3510 
AERO 3520 
afro 3530 
AER03540 
afro 3550 
AER03560 
AFRO 35 70 
AFRO350O 



AF703590 
AFRO 3600 



123 



ooo ooo 



FILE 



AERO 



FOR TP AM A MAVAL POSTGRAOUATS SCHOOL 



3129 CONTtNUE 

GO TO 32AI 

C CALCULATE R(J) ANC RPUJ FOP OGIVAL TAIL 

C 

3161 IF(OL.GT.O..AND.OL*L£*U J GO TO 3201 
«PITr(6 *3161 ) 

3ldl pC»maT(‘ 'REVISE INPUT. OL .G T.O . . A NO.DL *L£ .1 
GC TP root 

3201 PC = (rL*CL*0L»0L)/(2.»0L) 

00 32^1 J=J 6 P,NL 02 

P(J) * l.-RC^Sgr T(RC*PC-( X(J)-XL<L )»«2 ) 

Pt>( J ) =-(X( J)-XL*-CU/(RC-l.>R(J)) 

3221 C□^JTINU«= 



COf^PLSTICN 'lESSAGE 
32^1 ^PIT»=( 6 * 32 ol ) 

3261 tOP^ATI* • ,* GECHFtrICAL ARRAYS COMPLETED*) 
GO ro lOCl 



CALCULATE OUTPUT 



3231 

3284 

3286 



C 



3343 

3345 



3346 

3348 

C 

C 



8FTA=S0RT(XM*XM-1.) 

XI (11=0. 

T = x(2) ♦•2-(9FTA*P (2) 

IF(T.GE. 0.0)00 TO 3236 
v<RITF( 2 . 3264 )X(2)tP(2).8ETA,T 
FORMAT! » * *• 32 84* *4{El0.3f2X)) 

GO TO 3777 



T=S 0 RT(T) 
G( 1) =0. 
OAOI l)=0. 
0A2( I ) =0 . 
QNA{ 1) =0. 
CMA( IJsO. 
CAG=0. 
CA2=0. 

C-NA = 0. 
CMA=0. 

CAOl =C. 
CA03=C. 



00 3341 Js2fNLP2 
SllMC = C. 

SH«PsC. 

SUMPaC. 

SX = 0. 

TX=0. 

TR = 0. 

IF(P( J ) .NE.0.0)G0 to 3345 
WRITF( 2 * 3343 )J»X( J) *P( J) *RP( J) 

FORMAT ( • ' ,* 3343* ,2X, I 3. 31 2X,£l0.3n 
GO TO 3777 

XI ( J) sX(J )-8ET A«R (J ) 

HaX(J)-X( J-l) 

HAaXi (Jl-xi (J-ll 

T=xu )** 2 -(eeTA*p(jnR *2 

IF(T.GE.C.O)GQ TO 3348 

^PITM2»3 346 )X(J),R(J).RP(J1.8ETA*T,J 
FORMAT ( * * ,* 3346* *5( ElU. 3,2X), 13) 

GC TO 3777 
TxSORT(T) 

IFU.EQ.2 I GO TO 3321 



OP 3iCl I =2 f Jl 

T? = T 

T = 1X U )-X ({ I) ) **2-(3ETA*P (J) )**Z 
IFfT.LT.C.C) *0 TP 3302 
T^SQRTdl 



AER03610 
AERO 3620 
AEO03630 
AER03640 
AEP03650 
AFRO 3660 
AEP03670 
AFRO 3660 
AF'’03690 
AEP03700 
AFR037I0 
AER03720 
AFR03730 
APR03740 

AER037S0 
AER03760 
APRO 3770 
AE'^0 3 780 
AE903790 
AF30 3800 
AFRO 38 , 0 
AEP.O 3420 
AER03330 
AcRO 3840 
AER03650 
AFR03360 
AFO03870 
APRO 3680 
AFP03890 
AFRO 3900 
AEP039I 0 
^£*>0 3920 
AERO 3930 
AFR03940 
APPO 3950 
AFPO 3940 
\'^R03970 
AF^O 3980 
A PR 0 3990 
AFPQ4300 
APR04010 
A£R 04020 
AEP04030 
APP04040 
AER04050 
AFR04060 
APP 04070 
APR04030 
AER04090 
APR04100 
ACR 041 10 
AFR04120 
AFR04130 
AFR04140 
AEROJrlSO 
AFP04160 
AER04170 
AER04180 
AER04190 
AER04200 
AFR04210 
AFR04220 
AE'>04230 
APR04240 
APR04250 
AER04260 
AFPQ 4270 
AER04230 
APR0-*290 
AFR04300 
AER04310 
APP04320 
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(F nx( J )-xn I) >T ) ,NE .0,0 )GO TO 33(J7 

i J.I . X( JK A! a ) ,X( I-l ) ,R ( I )• R{J J , 3ETA*T, TP 

3?0 5 3305* .ZX.ZC ( 3*2XJ , 9(E10. 3*2X ) ) 

3 30 7 4=<<(JI-nl((-l)*.TP)/(X(J)-<{(()>Tl 

IF(A,L E.C.O) jO to 3302 
AsaLQG (A I 

P- TJ /3( J ] 

CsL^oo(J)eA 

0*U X<J»-Xl( I-l) )»TP-(XU)"Xni) )»T)/< 2, ♦RI J 
E»(X< J )-xn l^l )) »6/R(J i-0-3E TA*3FTA*A/2. 

FsF*-r3=TA-PFT \« A+^P ( J )*9 
= S'JmC ♦•C^Fp ( I ) 

<»jrtF = <^UMP *F *r,p ( ( i 

SU*n (- I ) 

SX-S X ♦'A«'cp( I j 

Tx = T X I ) ' 

iPsTS ♦•BFTA*^FT A* A*GP( { ) 

3301 CONTINUE 

3321 continue 

Tpa T 

T*0. 

IF( ( X( J) -Xl( J J ),NE,0,0)GO TO 3325 
AaU 

GC TO 3327 

332 5 As( X( Jl-xn j-U>TPj/(x( J)-XI(J)) 

{ F(A,GT ,0.0)00 TO 3329 
3327 R»l, 

C*l, 

0*1. 

p*u 

F*l. 

GC TO 3331 
3329 A*Af.OG(A) 

B=To/P( J) 
r*a-^s F U) ®A 

0*( X( J )-<n J-L) ) •TP/( 2.*R ( J i**2) 

F*( X ( J )-X l( J- I n ®B/PU i-D-6FTA®BE'^A*A/2. 

FsF^-PET A»PETA*A^PP( j)®B 
IF(C.NE.C.0)G0 TO 3333 
GC TC 3 3 31 

3333 IF(P.NE.C.0)GG TO 3337 . 

3 331 WRITP( 2.3 335 1J,X(J),XI (J-l),XI ( J) , TP . « U ) .RP( J ) , d£TA , A J b • C . 0, E • F 
3335 FCPmaTI* *.*3335* ,2X.(3.2X.6(E10.3,2X)/* S 7< E iq. 3, 2X) I 
GG ^0 3777 

3337 SU*'B-SU'-B4'b*G( J-U 

fp( j )*(PP ( J)-SumC1/C 

GP( J )s( 1. -SUHF-<;u.1o/ft ( J ) I /F 

G( J ) *G( J- I ) ♦HA®G? (J i 

SX= 5 X + A»F P( J) 

Tx* T X ♦P<= ; p ( J j 

T? s T 5 ♦ 5 "’^ A® BETA* A*GP( J) 

SP*P P ( J ( I .-SX) 

TPa 1 ,-TP-FP ( J ) ®TX 

GAO* I I . J) **2 )• ( I ,-SX )**2®( XM*SX)*®2 
GA’s-CX'^^Sx I, ♦^PPJ J) **2)*a,-SX) **2® 

♦ < 3 ETA*oETA-c.PU)»® 2 )®(TX**P)/ 2 .-ni.®TF)®® 2 )/ 2 - 

Cf.'As ( ( [.♦'Ri* J.J 1* < I .-SX )♦ XH»X'1*SX)*Tx 

0A0( J )* 2. ®P U ) »P PU )*G AO 

QA2(J)*2.«^ ( J i *0P ( J ) *GA2 

QKA ( J )^2. U i ®GNA 

QMA ( J ) = (0 ( J )®PO( J )^X( J) ) *ONA(J ) 

I F ( TNnSE.GT .OJGO TO 3339 
( F ( J.r,T . JAUGO TO 3 3 39 
CAOUCAOl <•( 0 A 3 i Ji ♦QAO( J -1 n*H/2 . 

3339 C'-'^;TlNMF 

CA03*CA03«-(CAOUn'OAO( J-t )) *H/2. 

CA2 = CA2MQA?(jJ«>aA2U-l))»H/2. 

C^'A»CNA4-( wN A( j )^»wNA u-l ) ) *H/2. 

C''A = C^ A^l :jMA( J ( J-l ) ) ®H/2, 

3341 CCNTPI'.IF 

CA0*CA03-rAfn 



AER04330 
AEP04340 
AFR04350 
AEP 04360 
AER04370 
AER04380 
AFP 04390 
AEP04400 
AFO 04410 
Ar'^04420 
AFQO 44 30 
AEP04440 
AER 0-.450 
aFR0^<.60 
AF9 04470 
AE«044o0 
AFP0440Q 
ACP 04500 
AFR 045 I 0 
A = P04520 
ABP04530 
AFO04540 
AEP0^5S0 
ACR04560 
AEQ 04570 
AER04580 
4 CO 04590 
4ERQ4600 
AFOO4610 
AF004620 
AFP04630 
AFOQ 4640 
AFP 0‘.650 
ACP04660 
AERQ4670 
AER04680 
AERO 4690 
AFR04 700 
A = P 04710 
AERQ4720 
AER04 730 
AF.P0 4 740 
AFR04750 
AER04760 
AEO04770 
A=onH780 
AFO04790 
AFQO '.800 
AFP 04810 
AEQ04820 
AFP 04830 
AEQQ48 VO 
AEQ04fe50 
AFP04860 
AEQ 04870 
AFP048flO 
AFR 04390 
A CD 0 4900 
ACP04910 
AHP 0^920 
AEP04930 
AE= 0 4940 
ACP049S0 
AER04960 
AFR04970 
AER 04980 
AEQ04990 
AFP 0 5000 
AEP 05010 
AFRO5O20 
AeQ05030 
AERO 5040 
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FILE: AERO 



FORTRAN A NAVAU POSTGRADUATE SCHOOL 



IFCCNA.NE.O.OJGQ TO 3355 \ 

WRIT’F(?, 2353 )CNA 
3353 FC'J^ATt* ' 3353* ,E10*3) 

GO TP 3777 
3355 XAC=C'^A/CNA 

w« ITFC 6 , 3361 J 

3361 FPRM4T ( • *,* CALCULATION completed* ) 
GC TO lOCl 
C 

c tabulate output 



338 1 «RITP( 2, 3A01 ) 

3^01 f OP-*AT ( • i* a6X** FORC E COE'=c I C I ENTS p qR A SLENDER, POINTED*/ 

>18X, •i^nCV Of ccvclJTICN in supersonic FLOy^*) 

WPIT^I 2, 342U X-. XL ,0L ,0L2 ,NL , NA, N8 . NC ,Nfi2 *NC2 , NT A8, I NOSE, 

♦ KNCSF. KT*A IL. KAP(o), KAP(7) ,CNA,CAO»CA?,XAC 

3^21 fpcmaM * C* .* XM-* ,F7.^, Tld, * xL^* ,F7. ^ , T36, *0L3* .F6.-V.T54, 

> *C! ?=* . F e.4/ 

• .'M , 13 . T 13, *NA^' , I 3.T36, *N3s* , I 3 , T5 4 , * NC « * , I 3 / 

, 'NA2 = * ,I 3, TIB,* NC2** . I 3 . T3 6 , * NTA 0 = ' , I3.TS4, * (NOSE** , I 2/ 

♦ • ',*Nn<:F COOr s' , A2 , T? T, ' TAI L COD E» • , A 2 , T S^. , 

•I'^=',2A‘./* '.'CNAs* ,F6.4, Tlb,*CAO=* ,F6.V, T36,*CA2 = *,F7.4, 

> TS4, *XAC=* ,F8.3) 

>^RITP t 2 f 3441 ) 

3A41 FC'=-^ATC*C*,2X,'J*,15X,*X*,15X,'R*,L4X.*RP*) 

OC 3461 Js2,NLP2,NTA8 

^FI"p(2.356 l) J, X(J ),PtJ) ,C PU ) 

3661 CONTIMUF 

wsirec 2, 35>*l ) 

3541 FCR'^AT I » C* ,2X, ' J* ,15X, *X* ,13X, *QNA* , 13X, *QAO* , 13X, *QA2* ) 

DC 3581 Js2.MP2.NTA8 

WR I TF 12,3 56 1 ) J. XU », ONA(J) .OAOtJ ) ,QA2IJ) 

3561 FrR‘3An» ',I3,VC16.4) 

3581 CPNTINUF 

GC ^0 lOCl 



PLOT OUTPUT 

WRITF OUTPUT FOP PLOT INTO 



file FTOlFOOl 



3601 WR ITPI I ,3621 ) XM, XL .OL .NL.NA ,NB,NTAB , KNOSE, KTA IL , 

♦ KAP.CNA ,CAU,CA2,XAC v 

3621 FOP*^AT(* '.3£l4.5,4l2X,t3),2I2X.A2)/* *,8A4/* *,4E14.5) 

00 3661 J = 1,NLP2 

WP I TP (1,3 64-1) XIJ ),R( J) ,QNA IJ) .OAOU) ,0A2IJ) 

3641 FpC8AT(* *,5E14.5) 

3661 continue 

C 

♦iR ITP( I , 2681 ) 

3681 FOP >^*r ( • *,*/♦*) 

wPITTCo , 3701 ) 

3701 FORMAT!' • , * TO OBTAIN PLOTS. FIRST ENTER QT TO QUIT PROGRAM,*/ 

*•* ','thfn issue the following command;'/ 

*-• *,8X,*FQP. RRINTPP plots* ,T36 ,' enter «CHARTS PRINTER***/ 

'.8X,*cCR plotter graphs* ,T36,‘ ENT ER •*CHARTS PLOTTER"*} 

GO TO lOCl 



QUIT PROGRAM 



3721 CCNTI'IUF 

WRITE! 6, 3741 ) 

3741 EORMATI* ', 'EXECUTION TERMINATED*) 
3777 STT3 
END 



A6R05050 
AER05060 
AFR05070 
AER05080 
AFR050D0 
AER05100 
AER05UO 
AER05120 
AER05130 
AFR05140 
AFR05150 
AEDO 5160 
AFR05170 
AEP 0 5180 
AER05190 
AEP05200 
AERO 5210 
AFP05220 
AFRQ5230 
AER05240 
AER05250 
APR05260 
APR 0 52 70 
AFR05280 
AER05290 
AFP05300 
AER05310 
AER05320 
AER05330 
AFR0S340 
AER05350 
APR05360 
AER0 5370 
AER053B0 
AFR05390 
AER05400 
AER05410 
AER05420 
AFP05430 
AFPO 5440 
AER05450 
AFR05460 
AEP05470 
AFR054B0 
Apoo 5490 
AFP05500 
AER05510 
AERO 5520 
AFR05530 
AERO 5540 
afro 5550 
AFRO 556v) 
AFR0557C 
AER055B0 
AER05590 
AER05600 
AFR05610 
AFO05620 
AFP05630 
AER05640 
AFO05650 
AEP05660 
APR J5670 
AER05680 
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G2. PROGRAM COWL 



ADJUSTMENT OF AERO FOR DO-LOOP ROUTINE 
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ooooooo ooo ooo 



FUEt CCWC 



FORTRAN A 



NAVAL POSTGRADUATE SCHOOL 



XM. 1.95 

GAsl.A 

nC 33 J*l *22 , 

X*<aXM<- 0 ,C 5 

IF( J.^ 0 . 1 )XM = 1.45 

CALL nr agN(XH.CAO> 

: call ORACWWl XH ♦GA,CDWW» 

IF( J.FO.nXMsl .95 
33 CCMINJF 
STOP 
ENO 

SUHPPtlTlNF HP AGN( XN,CAO» 

I»^Pl TCI T ocm.* 4 < 4 -h,Q-Z» , INTEGER<»^( I-N) 

DIMENSION X( 7021 tP( 70 ? ) • RP< 7 C 2 > • Xll 7 02 ) • F PI 702 > • QAD ( 702 » • <A PI 8 » 

calculate neoMFTptc^L arrays 

DATA FOR DL 2 : 

- 0 .l 2 C 77 a 6 . 50 FG. ;- 0 . 1770 * 9. 5 OrO . ;- 0 . 2 250*1 2 nFG. ;-o. 28 AO« 
*15 DEG 

OATA XL/ 17 . 50 /tNL/ 500 /.NA/? 03 /*NB/ 200 /,NC/ 6 /,Na?/ 61 / 

♦ , 0 L/- 0.0 7360 /,nL 2 /-O.l 7700 / 

A* XL/ 15. 7 6 /.NA/ 2 l 2 /*Mn/ 209 /,NC/ 79 /»NB 2 / 0 / tDL/- 0.2 7 fl 0 / 

SET UP TNptiT and CALCULATF X(J) 

2741 IF(fm 2 .FC. 0 )MC*Nl -MA-NH 
NC 2 *'''L-NA-Nrt-NC-NB 2 
DeiX=XL/FLnAT(NL) 

ALaO'^LXPFLPATI NA ) 

8 L* 0 FLX*FinAT(Kfl) 

Cl*OFl X«FLflAT( NCJ 
CL 2 *nrLx*FLD AT I nC 2 ) 

JAt=NA> I 
JAPa |*,L> I 
JRL*NA 
J 0 P*JAL^ I 
JCL*N^ ♦we*NC 
JCP*JCl 

J 02 ) *N A ♦A 0 *NC*'N 82 +l 
Jfl 2 P = JD 2 lM 
XI l)* 0 . 

P ( l)*7. 

NLPUML ♦'I 
NLP 2 *‘U >2 
XINL''? J*XL 

DC 27 ol J* 2 *NLPl , 

OFL J*FL 0 AT(J)-U 5 / 

XI J )*CELX*DELJ 
2751 CCMTNUF 

CALCUl ATF RIJ) ANO RPIJ» FOR OGIVAL NOSE 

3001 flAa(U«.AL **2 »/ 2 . 

00 3021 J*l*JAL 

tF|(RAtRA-(AL-XIJn** 2 ».G". 0 . 0 >GQ TO 3026 
wRtTPI ?* 30 ? 4 ) J*RA tAL.XIJ ) ,RI J) ,RPI J» 

3024 FOPHATI* • , * 3024 * , 2 X • I 3 • 2 X .5 I C 10 . 2 X M 
STOP 

3026 RI J)*SCRTIRA*PA-(AL-X(JI )** 2 )-«A 4 - 1 . 

PPIJ)»IAL-XIJI)/IRA- 1 .>R(J) > 

3021 CONTTNUP ^ 

3041 IFINA.IT.NU GO ^0 300 1 
3061 RINLP 2 |*U 
PPINLP ? )*.0 
GG TO 3241 

CAICULATF PU) AND pO(J» *HR CYLINDRICAL MlDS'*rTiON ♦ 

3001 IFINP.'O.OJ GO TC 3121 



COWOOOlO 
cnwoDO’o 
CPW3C030 
C0900040 
r PW00050 
rortO 0060 
C0Wf)u070 
rnwYioRO 
rn'^0 IDf’O 
)0100 
cnwo'Jl 1 0 
CDW00120 
T 1^0 
Cn'<.'nGT '»D 

CDWO J’ 50 
COWO 1 laO 
rnw(oi7D 
CPWO )IS0 
row ) )iQo 
cnw-^ )20D 
)^]o 

CPWJ 32PO 
CP40 )?.30 
rnuaD2!*0 
cnAO )?5o 
rnwo )T53 

CON-* }?7D 

rn.<0!)?«o 

cow TO ’90 
cnv 1 rr 
CDWODTIO 
cnwo >320 
CNW00330 
cn‘^0 )’40 
C0«*3'13S0 

ro ,0 ^’<60 

1370 
CPW ) 03 AQ 
CDWO )‘390 
COv^D >400 
COW0 3410 
cnwoo4’o 
COWOOOlO 

rnwo )'^40 

rnwo 3^»5D 

rn -^0 0460 
cnw.)04 70 
r nwo 34HD 

cnuoo^90 
cnwo J500 
cnwo 1510 

cnwo 3520 

rnwo 05 10 
cnw00540 
cnwo 1550 

rflAl ) 5^0 
Cn^OC570 
rowa 35 AO 
rnwo )5'io 
rn^O 1600 
cnwo 1610 
cnwo )620 
rnwoo630 
cnwo 1640 
cnw00550 
cnwo 067 0 
rnwo )6 7 o 
rnwo T6A0 
cnwo iMo 
r nw 1 j / 10 
cnwo : 71 0 
rnwi 1 
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ooo r»oo ooo 



FILE: 

3101 

3121 

3l<»l 

3123 

3125 

3127 

312«> 

C 

32<>l 

C 

C 

3281 

328V 

3266 

C 

3343 
3 34 5 



COWL FORTRAN A NAVAt, POSTGRADUATE SCHOOL 



RP(J)*.0 




C0W00730 


CONTINUE 




C0W00740 


IFMNA^NEKLT.NDGO TO 3121 




CQW00750 


Rl^LP2 1 =*U 




C0W00760 


i^PCNLP? ) =.0 




Cnuoo770 


GO TC 3241 




C0»J00730 


CCNTIN'JF 




cnwoo79o 

C0*^00800 


calculate R(J) ANO RP(J) FOR CONICAL TAIL 




cnwooaio 

cnw 00820 


»PT*-0L/CL 




C0W00830 


BXL*AL>3L*CL 




CGW00840 


DO 3141 J=J6R.JCL 




cnwo0850 


R( i) »l.^RPT*( X(J )-RXL+CL) 




C0W0C060 


>?P( J ) =RPT 




CQW008 70 


CONTINUE 




cnwoos50o 


IF ( (MAfNg^NC ) .LT.M )C0 TO 3123 ’ 




COW00890 


R(NLP2) =l.-0L 




cnwoQ9oo 


R P (NLP 2 ) -ffPT 




cowo 'mo 


GG TC 3241 




COWO092O 

cnw00930 


calculate r(j) anc rp(j) for second hiosfction 




COWOO940 

C0W00950 


IF(NB2.£C,0) GO TO 3241 . 




cowo 0960 


00 3125 J«JCR,JB2L 




COW00970 


R(J ) 3 l.-OL 




CGW0G930 


RP( J ) «.0 




cowo 0990 


CONTINUE 




COWO 1000 


IC({MA4.N04.fx:<>NB2).LT.NL)GO TO 3127 




cowo lOlO 


R (NLP2 ) = l.-OL 




cowo 1020 


R P(NLP2 )s.O 




COwO 1030 


GO TO 3241 




cowo 1040 
COW) 1050 


CAlCUIATE R(J) ANO RP(J) FQR SECOND CONICAL TAIL 




cowo 1060 

COWO 1070 


RPT2*-0L2/CL2 




cnwoioao 


00 3129 J*JR2R.NLP2 




cowo 1090 


R (J ) = l-OL i*R PT2* ( X (J )-XL«*CL2 ) 




cnwo i iQO 


PP(J)-PPT2 




CQWOIIIO 


CONTINUE 




cowo 1120 
cowo 1130 


ccntinup 




rowo 1140 


CALCULATE OUTPUT 




cowo 1150 
cowo 1160 


8ETA = S0RT( XM*XM-U) 




r.nwo 1 1 70 


XI ( 1 j = 0. 




cowo 1 180 


T*X(2)*»2-(8ETA*p(2J )**2 




cowo 1190 


IF(T.GE.C.O)GO TO 3286 




cowo 1200 


WR(TF{ ? , 22 84 ) X(2) »R( 2) .BETA* T 




cowo 1210 


FnoMATl * 3284*,4(E10.3#2X)) 




cowo 1220 


STC® 




cowo 1230 


T=SUFT( T) 




C0W01240 


QAC( 1) -a. 




COWO 1250 


CAOl *C. 




COWO 1260 


CAC3aC. 




cnwo 1270 


CA0=0. 




COWO 1230 
C0W01290 


00 334 1 J»2,NLP2 




cnwo 1300 


SU*<C = C. 




cowoi 310 


SUMFsC. 




C0W01320 


Sli'^R *C. 




C0W01330 


Sx=0. 




COWO 1340 


TX=0. 




COWO 1350 


Tff = 0. 




cnwonso 


IF(R(J I.AE.O.O)GO to 3345 




cnwo 1370 


WRITF(2.3343 )J#X(J),R(J),RPU) 




cowo 1330 


FORMAT ( • S' 334 3* ,2 X, I 3. 3( 2X,E 10, 3) ) 




COWO 1390 


STCP 




cnwo 1 400 


x( ( J J *X(J )-PRTA*a(J ) 




cnwo 1410 


H=x ( J )-x( J- 1 ) 




cowo 1420 


HA*x I (J)-XI { j-l) 




C0W01430 


T*X { J 1**2- ( 8ETA*» (J J ) •♦a 




COWO 1440 
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puet cQviu 



FCaTRAN A NAVAL POSTGRAOUATP SCHOOL 



3346 

3348 

C 

c 



3302 

3305 



3307 



3301 

3321 



3325 

3327 

3329 



333 1 
3 335 

333 7 



3339 

3341 

3344 



IF(T.Ge.C.O»GO TO 3348 

W«ITF(2 ,3346 )X(J)tP( J) ,RP(J) ,6FTA,r, J 
format ( • • 3346* ,5(Pia.3,2Xi, 13) 

STOP 

T»SORT(T) 

IF(J.EQ.2) GO TO 3321 



Jl = J-l 

00 33CI !-2,Jl 
TP*T 

T = (X( J )-X(( I ))*«2-(B€TA*R(J) 1**2 
IF(T.LT .C.JI GG TO 3302 

T=SOmT(T) ^ 

I F (( X( J )-Xl( I I I .NF.0*0)GO Tn 3307 

wo (TP(2 , 3305 )J , ( , XU) ,XU ( ) ,X( l-l ),R ( n, R(J),^R( J) ,a€TA,r,TP 
FORMAT ( • • 3305* ,2X,2 ( { 3 ,2X),9(E10. 3»2X ) ) 

STOP 

A=(X( J J-X<( I-l)>TP)/(X(J)*Xl ( I)4-T) 

IF( A.L E.C.Ci) GO rrj 3302 
a = ALOG(A ) 

P=(TP-7)/P( J) 

Cs-j^^PP ( J ) «A 
SU'<=SUMC^C*Fpn ) 

SX = SX>A*FP( ( ) 

CCMTINUF 

CONTIMJE 

TP*T 



T = 0. 

IF((X( J)-XU J) ).NF.O.O)GO TO 3325 
A=U 

GO 3327 

As( X( JI-XK J-l )»T?)/(X( J)-XUJ)) 

IF(A.GT.C«0)GO to 3329 
B=l. 

C* 1. 

GO TP 3331 
A=ALOG(AI 
fl = TP/R(J) 

C=StRP( J)*A 

IFCC.NC.C.OIGO TO 3337 

WRIT6( 2 , 3335 ) J,X( J) , XI ( J-l) , XI ( J) ,TP »P { J ) ♦«?( J ) , BET A, A, B, C 
FORMAT ( • »,* 3335* ,2X,I3,2X,6(F 10. 3,2 X)/* * ,4( 6 lO . 3 , 2X ) )' 

STC® 

FP( J) = (RPCJ )-SUMC)/C 
SXaSX ♦A^F P( J ) 

GAO* l.-( I .^PP( J)»*2 )♦( l.-SX )*»2^( XM»SX)»<«2 
0A0( J ) = 2.*R ( J)»OPU)*GAO 

iF(j ,GT.JAL )G0 TO 3339 
i f AOl *C AOl»(QAO( J |>0A0( J-m*H/2. 

CONT INUF 

CAi)3*CA03«' « QAO(J ) HJ AO ( J- 1 ) ) ♦H/ 2 . 

CCNTINUF 



CA0aCA03*CA0l 

WRITF( 2 , 3344 JXM, XL, Nl, NA.NR, NC ,NB2, 0L2 ,C AO 
FOP*^AT ( • *,F5.2.2 x,F7,4,S( 2X.I 3),2(2X,F7.4)) 
RETURN 
END 



SUBROUTINE ORAGUW(XM,GA,COWW ) 

P(s4.*A TANU . ) 

RADFG»o I/IBO . 

0A7A TWO/5./ ,OREF/5./,8/0.09525/,C/0.0635/,T/0.01/ 
AR6FsP( •(ORFF*O.C254j*»2/4. 

CC ALL LENGTHS IN METERS. 

Tw^TkO*R AOFG 

XMLsXM 

XM2*XM 

CALL PRAM(XH,GA,ANI0) 

AN I IsAM I C-TW 
\NI2*AN I C*-TW 
CALL P^cSIXM .GA.POPT) 



C0W01450 
COWO 1460 
COW01470 
COWO 1480 
CnW0I490 
COWO 1500 
C0W01510 
COWO 1520 
COWO 1530 
COWO I 540 
COWO 1550 
COWO 1560 
C0WO157O 
CnwOl580 
COWO 1590 
CO^O 1600 
CTwO 1610 
COWO 1620 
Cnwo 1630 
COWO 1540 
COWO 1650 

Cnwo 1660 

COWO 1670 
C0W01680 
cnwo 1690 
COWO 1700 
C0W01710 
cnwo L 720 
COWO I 730 

cnwo 1 740 
cnw0l750 
cnwo I 760 
cnwo I 770 

COWO 1780 
cnwo 1790 
COWO 1800 
COWO 1810 

cnwoid 20 

COWO 1330 
COWO 1340 
C0W01850 
COWO I860 
CCWO 13 70 

cnwo 1800 
cnwo 1890 
cnwo 1900 
cnwo 1910 
COWO 1920 
C0W0L930 
C0W01940 
cnwo 1950 
C0W019S0 

cnwo 19 70 
cnwo 1930 
COWO 1990 

cnwo 2000 

C0WO2010 
cnwo 2020 

COWO 2030 
C0rt02G40 
COW02050 
COW02060 

cnwo 20 70 

CnW02030 

CPW02090. 

cnwo 2 ioo 

C0W021 10 

Cnw 02 l 20 
cnwo 2130 

C0rt02l40 

CnwQ2l50 

cnwo2i60 
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FILE* COWL 



FORTRAN A NAVAL POSTGRADUATE SCHOOL 



C 



CALL APOANT( XNl,GA*ANtl) 
call PRFS(XN1 . GA.PIPT) 
CALL APOANT( X:-*2.GA,ANI2) 
call ®ReSlXM2.GA,P2PT) 



BP = a-C /? / T( 

CDww* <2 ./(CA *X ^^••2) ) *( PIP'T/POP T-P2PT /POPT )*(T*BP /AREF) 
«RITF( 2 t33^» x‘4,XMl,XM2 ,POPT»Pl PT, P2 P T , Tw ♦ 8«rtP, C , COWW 
334 FORMATI* 334» ,6( 2x,Fe.5)/* »♦' • * 5 ( 2X • F 8, 5 ) / 

♦ • 75 I * ♦•)) 

ocTifC'j 

ENO 



SLSRCJTINF 3RANT ( I fGA. ANt J 
aFTAsSORT(XM I**2-l. ) 

GARA^'sSORTh oA-l. ) / (OA^l . ) ) 

AN I* AT AN (GAR A T *86 TA ) /GAR A T- A TA NO ET A ) 

rftupn 

ENC 

SLHROUTI^E APRANT(XMI .GA.ANI ) 

8ETA*SOR T( XN ) 

GAPATaSQRT(( GA-1 . ) / ( GA*l . ) ) 

I 8FTa=0 

552 F = ATAN( gap it *8 ct a )/ GAR AT- A TAN( SETA) -AN! 

FP^l. / ( I .♦•(GAPaT*06Taj •«2)-l./ (l.OETA**2) 

IBFTA= I BETA^ I 

tP ( TA .GT. L2 )GO TO 558 

IFCFP.FQ.O.OJGO TO 558 

a^TANxRETA-F /FP 

IF((ABS(eETAN-8ETA)).LE.(Ue-05))Ga TO 556 
BETAsiCT AN 
GO TO 552 

556 XMIsSOR T(8Er A»*2^Ui 
acT.)p^ 

558 wRITP(2t560) I B ET A , AN I , F • FP , A 

560 FGRMATi' • * • 560« t 13 t4( 2X , <=U .4 ) ) 

STCP 

END 



SUBROUTINE »RES(XMI ,GA,P IPT) 

P( PTs( I .XGA-l .)/2. •XNI**2)**( -GA/(GA-l.) ) 

return 

END 



COW02170 
COWO210O 
COW02IRO 
C0WO22OO 
CQW02210 
C0W02220 
C0W0223O 
COW022AO 
COW02250 
COW02260 
COW ) ^’70 
CCW)22B0 

C0W022O0 
C0W02300 
C0W02310 
CGW02320 
COW02330 
COW 02 340 
Cnwo 2350 
C0V02360 
COiJ0 2 3 ^0 
CCW023H0 
C0W02390 

cnw 02400 
cnwo 24 10 

CDwO 2420 
COwO 2430 
CDWO 2440 
C0W02450 
C0WO2460 
cnwo 2470 
CDH02480 
C0WU2490 
C0W02500 
C0W02510 
CnwQ2520 
C0W02530 
C0W02540 
Cnw02550 
C0W02560 
cnwu 2 5 70 
C0W02580 
C0W02590 
C0W02600 
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G3. PLOT ROUTINES FOR USE WITH AERO [9] 
-PREPLOT P (PLOT ON PRINTER) 
-PREPLOT G (PLOT ON PLOTTER) 
-CHARTS (CONTROL) 
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PILES P96PLOTP FORTRAN A NAVAL POSTGRADUATE SCHOOL 



//AMTC3250 JOP (3250,0258 )**AMICHAI X29 J5 • ,C LASS^A 
// EX6C FRTXCIGP 
//POPT.SYSIN CO ♦ 

c 

C SPECIFICATION STATEMENTS 



IMPLICI T RFAL»4( A-H,n-X) , lNTEGEft*4( I -N ) 

01 MENS ION XI 502) »Ri 502 I # ONA ( 50 2) ♦ OAO ( 502 ) iQA2( 502) *KAP(3) 



REflO(5, I COO) XM* xLf OL.NL .NA, NB ,NTA<3» KNOSE*KTAIL* 

■►KAP.CVA, CA0tCA2tXAC 

1000 FCBMAT ( I x,3E 2X, 13 ) .2( 2X. A2)/l X.aA4»/lX.^El4.5) 

NLP2 = Nl ^-2 
CO 1040 J = UNLP2 

RPA0(5,10 20) X(J),AU ) .QNA( J) ,0A0 (J). QA2U) 

1020 PC-’M AT( IX .5E14.5 ) 

1040 CCNTIMIJE 

ITFl 6. 1060 ) 

1060 FDR'^AT < • ' SU'^'^APY GF RESULTS*///) 

WP ITF ( 6 , 1080 ) xm,xu* 3L»NL.NA,N8»NTA8,KN0SE*KTAIL*XAP(6)*KAP(7)* 
>CNA .r AO .CA2t XAC 

1080 FOPmaT ( » C* XM*» ,F7.4,Tia,* XL* » .F 7.4 *T36. *0L** *F6.4 iT 54,»NL** t 13/ 
«- *0* . *NA = • , 13 » T 18, *N8 = * , I 3,T36t *NTA8« * » 13 .154. ’NOSE CODE** . A2/ 
^*0*.**^AIL COOEs* .A2.T3 6, *10**, 2A4 / * 0 * , *C N A* * ,F6,4 ,T18, *CAQ* • , 
^F6.4.T36,*CA2* •.P7.*,,T54, 'XAC* *,F8. 3 ) 

I TP ( 6 , 1 100 ) XM»kAP(6),XAP(7) 

IlOO PGPMAT ( * I* ,36X,*RA01US R VS X XMs*,F7.4,» I0«*,2A4///) 

CALL PLOTPlX.R .NLP2.4) 

WRITF(6. U20) 

1120 FORMAT! * 1* .3oX,*NCRMAL FQPCE COEFFICIENT QNA VS X*/) 

ITF I 6, 1140 ) CNA ,X AC . < *,P ( 6 ) , K API 7 ) 

1140 FORMAT! »C* .36X.*r\A* * ,F7. 4. • XAC»*.F9.4,* I0**.2A4//) 

call PIGTPIX ,0NA .NLP2.4) 

*RITP( 6, ll60 ) 

1160 FC«'*AT( • 1*,36X. ‘AXIAL FG®r. E COEFFICIENT OAO VS X*/) 
wR ITr ( t>, 1180 ) CA0,;<AP(6) ,KAP (7 ) 

1180 FORMAT! *C* ,3oX,*rA0=*,F7. 4, < I0»*,2A4//) 

CALL f^LOTPlX ,OAO,NLP2,4) 

WR ITFI 6 , 1200 ) 

1200 FORMAT I • 1* ,36X,*AXI AL FORCE COEFFICIENT QA2 VS X*/) 

WRTTF! 6 , 1223 ) C A2 ,K A p ( 6 ) ,K AP ( 7 I 

1220 FORMAT! *C* ,36X.*CA2**.F7, 4, • I0**,2A4//) ' 

call PL0TP(X.0A2.NLPP,4) 

STOP 

ENO 

//GO.SYSIN DO • 



PREOOOlO 
PRP00020 
PR E 00030 
PRF00040 
PRE00050 
PRF00060 
PRPQ0070 
PRE00080 
PPFO 0090 
PREOOlOO 
PRPOOIIO 
POP00120 
PR^OOISO 
PRF 30140 
PP=00150 
PR=00160 
PRF00170 
PRP0C180 
PRpQ 0190 
PRPO J200 
PPf 00210 
PRP00220 
ORFO 3230 
PPE03240 
PRF00250 
PRE00260 
POC002 70 
PRF00280 
PRFQ02PQ 
PP‘=Q0300 
PRC00310 
PRF00320 
PRPQ0330 
PR«=0 0340 
PRFO 3350 
PR=00360 
PRB0Q370 
PP = 0 3300 
PRF00390 
PRF00400 
PRE00410 
PRPQ0M20 
PRC00430 
POF00440 
PRF0C450 
PRE00460 
PR PO 94 70 
ORF00430 
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FILE: PREPLOTG FH^^TRAN A NAVAL POSTGRADUATE SCHOOL 



//AMTC3250 JOB ( 32 50 • 0258 J , • AH ICHAI X29 3l»* ,CLASS«A 
// SXFC FPTXCIGP 
//FORT.SVSIN CO * 

SPFCIF TCATION STAT6HFNTS 

implicit R*^AL»A( A-H,0-X» , INTEGFR*A( I-N» 

C 

DIMENSION X( 50 2) t'^I 502) . ONA t 502 ) t OAO < 502 ) tOA2( 502) .KAP(8) 

C 

PP/J0(5, I CCO) XM, XI ,0L.NL .NA, N8 ,NT A8, XNOSE.KTAIL. 

♦ KAP.CNA ,CAOf CA2, XAC 

1000 F CRH AT (lX,3El<».5.^(2X,13)f2( 2X , A2 ) / I X, 8AA/ IX, AE 1 A. 5 I 
NLP2 = NL<‘2 
DO lOAO J*1,‘!LP2 

SFAOI 1020 ) X (J I ,R (J ),QNA( J ),CAO (J ), QA2(J) 

1020 RCRMaK U,5EU.5) 
lOAO COMINUF 

WRITC( 6,1060 ) 

1060 FCR*»AT (• l» .» summary CP P6SULTS*///) 

wR ITF( 6 , lOsO ) XM,XL,0L,NL,NA,NB,NTA8,KN0Se,KTAIL,KAP(6»,<AP(7), 

♦ CNA ,CAO ,CA2, XAC 

L080 FQbmaT ( »C* ,* XMs* T13, • XL=^ » ,F7.<, ,T36. 'OL^* , F 6 , T5A , • NL* » , 13/ 

♦ *0* , »NA=» , 13 ,T 18, »N6*» , I 3,T36, *NTA3a » , 1 3, T5^, ’NOSE CODE** ,A2/ 

♦ * 0 * . *T1 IL rOOEa* , ^2 ,T36, » to* * , 2A4/»0 * , *CNA** ,F6.A ,T18, •CAO**, 

♦ F6.<,,T3 6, *CA 2= *,F A,T5<*, • XAC= • ,F8.3 ) 

CALL ^LnTG(X,P ,NLP2,l, l,D,KAP. 32. MAOIUS R • , 3, 0 . , 0 . , 0. , 0. , 8 • 5 • 6 . ) 
CALL PL OTr,{ X ,CNA.NL P2, 1 , I ,0,KAP, 32, * ONA* ,3,0.,0.,0.,0.,3.5,6,J 
CALL RLOTGIX ,OAO,NLR2, I, l,0,XAP, 32, * OAO* , 3,0.,0. ,0. ,0. ,3.5,6. I 
CALL PLOTCIX ,0A2, NLP?, 1,1,0, KAP, 3 2, ♦QA2» ,3,0., 0.,0.,0., 8.5,6. ) 
CALL PLOT (0.0,0. U,V9V) 

STCP 

FNO 

/♦ 

//GO.SYSIN 00 * 



PRFOOOlO 
PRF00020 
PRE00030 
PRFOOOAO 
PPPQ0050 
PRF00060 
PRCJ0070 
ORCQ0080 
PRP00090 
RRPO 0100 
PRFOOllO 
PRP00120 
PPP00130 
PPPOOl^O 
PRF00150 
PRP0C160 
^R=00170 
PPPO 3180 
PRF0G190 

Po=00200 

PRF00210 

PR=00??0 

RRR00230 

PRF002^0 

PRPQ0250 

PRF00260 

PRP0J270 

PRF002B0 

OPP00290 

PRE0U300 

PR=00310 

PRP00320 

PRP00330 

PPPOOIAO 

PPP00350 
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//AMIC3250 JCe ( 3250f 0258), » AMICHAl >2936 • ,CL ASS*A 
// £XeC FRTXCLGP 
//FCPT.SYSIN CC • 

/• 

//GC.SYSIK CC • 



//AMC3250 JC8 (3250 ♦ 0258) ,♦ AM ICHAl X2936 * ,CL ASS« A 
// exEC FftTxCLGP 

CCC: THIS IS the FORTRAN IV H tXTENOEO PRCCSnURE THAT 
••• COMPILE, LiNK-eOIT AND EXECUTE mlJH PLOTTER 

• ♦ « 

PRCC IMSL^SF 
PjM^IFEAAft ,SEGICN=6A0K 



PREOOOlO 

PRE00020 

PRE00030 

PKt00340 

PRE00350 



IS A 



XXFRTXCLG 
XXFCRT EXEC 



PREOOOlO 
PRE00020 
OOOIOOOO 
0002000 I 
OOOJOOOO 
OOOhOO I M 
OO050O 



6 


X X s><term 


QO 


CUM'^Y 


000 /OOUO 


7 


X X < V S L T 1 


00 


UN I T = SYSCA, SFACE = (CYL . (1 ,2) J ,OC8*0 L KS I ZE»3 Ao5 


OOooOO 1 A 


8 


XX SYSUZ 


00 


UNI T^SYSCA. SPACE = ( CYL . (2 .2 ) ) .0 CB*0 L K S I Z E-20 A3 


OOUSOO 1 


9 


X XSTSLIN 


00 


OSNsCCej . SPACE » (30^0 ♦ ( AO.AO) ) ,UNIT«SYSCA, 


♦•uuiGuooa 




XX 




CISPa(>CCtPASS) . 

CC8={ ELKS IZE^30aC, LRtCL*<10,RECFMsFaS,BUFNa*U 


<»ouiioo()a 




X X 




O0L20OOO 


10 


/ /FCR T.SYSIN 


OC * .ccE^eLKSiztaeo 


PKEOOOaO 


ll 


XXLKEC EXEC 


oGM^IEWL ,CCNC*<9,LT.FCPT) 


00130017 


12 


X X ST SFR INT 


OQ 


SYSCUT 


GO lAOOOo 


13 


XX STSLMOO 


00 


SPACH*(CYL , <A,4, 1) ) ,DSN-=CLUO(X) ,01 SP*( .PASS ) , 


♦-O0 IbOOO i 




X X 




UMI T=SVSCA.CCS*aUFNGsl 


UGioOuOO 


lA 


XXSTSLTI 


00 


USfi^SSYSLri ,SPACE=( 102A, ( 120, 120) ) ,UNITa$YSCA, ^ 


♦•O017OJOO 




X X 




CCa*dUFNC=* 1 


JO 1 oOOOO 


15 


X XSYSL IN 


00 


OSN = CCaj,CISP^(CLa,OELETE ) 


JOiOOJOO 


16 


XX 


00 

uO 


C5N=SYSl .PRCCL la (VMAPP) ,0ISP»SFR 


002000U2 


17 


X X 


C0\aM6=SYS (N 


«)U210U03 


18 


XX5YSLIB 


00 


UNI Ts 3 35(1,0 ISP^SHR ,VLL*SE!<=MVS 002 ,CSN= SYS 1 .PP.FURTLIB 


0022JO 16 


19 


XX 


00 


LNI T=3 350 ,0 (SP=SHR , VC L- SER^MV S003 , C SN*S YS 1 .M PS LI 6 


t)023O'J l*t 


20 


XX 


CC 


UNI T- 3 3 50,0 (SP^SHR , VC L =S tR-MV S 003 ,0SN*SYS3 . (ySU-UMSL 


OGZhOOiS 


21 


XX 


00 


UNI T= 3 3 50,0 ISP »SFR ,VCL=S 6 R = MVS003 ,CSN= SYS 1 .V TEC PLOT 


U025OO 


22 


xxcc EXEC 


PG«’=*.LKEU. SYSLMCD.CCNCs^I ( 5 . L T , LKE U , < 5 , LT , F CR T ) ) 


l)G2ftO() 1 A 


23 


XXFT05FO01 


CO 


CONAME^SYSIN 


1)02 70000 


2A 


XXF TCtFCOl 


00 


SYSCLl =• 


002b0OiK> 


25 


Xtf TC7F001 


uo 


SYSCUT=a 


0029OU00 


26 


XXVfCTRl 


00 


aSN = cevECTHl ,0 1SF=(,HASSJ ,SPACE*(THK,< 1, 1) ) ,LNIT=SYSOA 


003UO009 


27 


XXVECTR2 


00 


CSN=?CCvrCTR2 ,01 SF= ( .PASS) , SP AC ( C Y L ♦ ( 5 , 5 ) ) .UNIT-SYSOA 


003 iOi)00 


28 


XX F LC TLGG 


oc 


SYsnuT^^ 


U032OOO9 


29 


XXF T15F001 


DO 


Q0NAME=PLCTP ARM 


0033000S 


30 


//CC ,SYSIN 


QC 


• .CCe^HLKS IZE*aC 


PKcOOjSO 


31 


XXFLCT EXEC 


PGV= ( F VMAPP ,CONO-( A,LI ) 

LNlT* 3350 ,OlSP=SHP.,VCL»SERaMVS003,CSN*SYSl .VIECPLOT ^ 


003hOQG9 


32 


XXSTEPLia 


00 


003500 Lh 


33 


XXFLCTLOG 


oc 


syscut** 


00360009 


3A 


X X STSVEC TR 


00 


SYSCUT*< A,,5555) \ 

CISP= (CLO.OCLE Tfc J ,QSN=CCVECTRl 


00370009 


35 


XXVECTRl 


00 


00360009 


36 


XX V6CTR2 


00 


CISP«(CLC,OELETE ) ,CSN*CCVECTR2 


0C390009 


37 


XXVECTTAPE 


00 


DUMMY 


00AOO009 
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FUE: charts exec a ^iAVAL POSTGRADUATE SCHOOU 



CCCNTPCl ERROR I'' 

CIF .CL EO .PRINTER CGQTn -PRINTER 
CIF .Cl FO .PLCTTcr cgoto -PLOTTER 
CTVPF REPEAT. SPECIFY PRINTER OR PLOTTER. 

CEXIT 

-PRINTER CCONTINUF 
CN&'^E = PRBPLOTP 
CGCTQ -00 

-plotter CCHNTINUE 
CNAME - PPFPLOTG 
-CC CCONTINU'^ 

COPYFILE C^AME FORTRAN A FILE FTOlFOOl A PLOT FORTRAN A (REPLACE 
EXEC SUBPIT PLCT FORTRAN 



\ 
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G4. Results from AERQ/COWL 



G4.1 The symbols used in programs AERO and COWL are defined in AERO 
and are presented in Figure G4.1. The values of a, , c-j , b2» C2 are 
normalized with respect to r-j . NA, NB, NC, NB2, NC2 are the appropriate 
numbers of points used in the program. The cowl angles a-j , 02 were 
selected as 20°, 9.5°, respectively. 

To create flow at the first cowl in Figure G4.1 which is the same as 
for a ramjet inlet, an extension to the body was used. The extension con- 
sists of the cone of length a and cylinder of length b-j in Figure G4.1. 

The cone angle, 6-| , is 8 degrees. The value of b-j was varied until the 
pressure at the first cowl was equal to ambient; a value of b-j equal to 
7r^ gives this condition. 

G4.2 The normalized values of the various variables which were 
selected are as follow: 

a = 7.11 , b-| = 7. , = 0.2 , b2 = 2.13 , C2 * 1.06 units 

- DL = 0.0736, - DL2 = 0.177 units. 

The appropriate numbers of points are: 

NA = 203, = 200, NC =6, NB2 = 61 , NC2 = 30 

The dimensional values are: 

c^ = 0.40, b2 = 4.22, C2 = 2.1 inch 
r^ = 1.98, T2 = 2.13, r3 = 2.48 inch 

G4.3 It was found that the cowl drag coefficient is sensitive to the 
magnitude of projected cowl area. For example, another combination of 
these variables: 
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a = 7.11, b.| = 7., = 0.38, = 2.16, C 2 = 1.46 

NA = 196, NB = 193, NC = 11, NB 2 = 60, NC 2 = 40 
-DL = 0.138(20°), -DL 2 = 0.243 (9.5°) 

Gives higher values of the cowl drag coefficient; 



Mo 1 


3.0 


2.3 


G4.2 


0.0732 


0.0953 


G4.3 


0.1214 


0.1562 



Therefore, attention was made to select an inlet shape as smooth as 
possible. 

G4.4 AERO can produce also graphical results. A typical example is 
illustrated in Figure G4.2. 

G4. 5 Computer Program Users Guide 

Use opening cotimands and compilation similar to that defined in 
F2.1, F2.2. 

The routine to run AERO on the terminal is defined in program 

itself. 

The routine to run COWL on the terminal is as follows: 

FILEDEF 02 DISK OUT D (RECFM F BL0CK 132 PERM 

EXEC RUN COWL 

PRINT OUT 0 I 

[ optional 

XEDIT OUT 0 ’ 
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Figure G4.1 Geometry for Calculation of Cowl-Drag- 
Coefficient (Programs AERO & COWL) Showing 
Definition of Symbols 
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Figure G4.2 Typical Results from AERO 
XL-18.88, NL=600. RA.= 188. NB-185. NC»20. NB2=56. NC2-51 



APPENDIX H: RESULTS 



HI. SUMMARY 






Ac/A 
5' r 



A^/Aq 



A2/A0 



As/A^ 



0 . 25 - 0.40 0 . 25 - 0.40 0.47 0.887 0.426 9.5 15 

25 

35 

40 

45 

65 

80 

0 . 25 - 0.40 0 . 25 - 0.40 0.50 0.827 0.426 9.5 7 

25 

45 

65 

80 



0 . 25 - 0.40 


0 . 25 - 0.40 


0.42 


0.827 


0.426 


9.5 


45 


* 




0.58 










0 . 25 - 0.40 


0 . 25 - 0.40 


0.47 


0.75 


0.426 


9.5 


45 








0.827 














0.91 








0 . 25 - 0.40 


0 . 25 - 0.40 


0.47 


0.827 


0.27 


9.5 


45 










0.32 














0.47 






0 . 25 - 0.40 


0 . 25 - 0.40 


0.47 


0.827 


0.426 


6.5 


45 



11.5 
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